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BILL TAUFMAN, a scientist with the Research & 
Development Division, studies edible plants grow- 
ing in pressurized glass jars under simulated 
high-altitude conditions. Tests thus far indicate 
that the pressure for optimum growth is of the 
order of 8 to 10 psia, some plants growing con- 
siderably faster at this pressure than at the earth 
sea-level pressure of 14.7 psia. 
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zero gravity on plant growth. Although zero gravity can- 
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down growth, showed only one-half the development of 
normally grown plants. 

Every facet of flight technology, from advanced 
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computer-based trajectory studies for planetary recon- 
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and Development program. Distinguished advances have 
already been made in space guidance concepts, in plasma 
and nuclear propulsion systems, in radiation physics, in 
new materials and processing techniques, and in proto- 
type development of hardware (as an example: hydraulic 
systems that operate reliably up to 1500F.) 
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for a Fleet Ballistic Missile, the first Polaris test vehicle was $ 
cessfully fired. Firings continued throughout 1958. By Septem™iyn4,) 
the test vehicles were close to the final configuration of the Pola ice y,. 
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Sometimes forgotten during the thundering ascent 
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DVANCES in ramjet combustion or, more broadly, ram- 

jet propulsion, are likely to have a considerable impact 
on our way of life in coming years. Today, ramjets are in use 
only in guided missiles, small helicopters and a few re- 
search airplanes. But the quest for speed in terrestrial trans- 
port has not ended with the recently introduced commercial 
jet transports or the Mach 2 military planes. Already the 
aircraft companies are looking forward to supersonic trans- 
ports (1).1. Predictions of hypersonic airbreathing aircraft 
appear frequently, e.g. (2,3), and ramjet propulsion will un- 
doubtedly be used in such aircraft. 

When one considers the time required for passengers to get 
to and from airports, or the time required to assemble and 
load and then unload and deploy military personnel and ma- 
terial, there appears to be no great advantage in reducing trans- 
port time below an hour or two. Such transport times—for 
intercontinental distances—are within the reach of airbreath- 
ing aircraft, but do call for speeds of Mach 3 and higher. A 
Mach 7 transport could make a nonstop flight of 6000 nautical 
miles in 2 hr. The ramjet engine would be at its best in this 
flight speed range. A slight extension of present technology 
indicates that the ramjet will reach its greatest overall ef- 
ficiency somewhere in the range Mach 6 to Mach 8 (2,4). 
This means that direct operating costs for hypersonic aircraft 
using ramjets for high speed phases of flight will probably be 
no greater and may be smaller than those of the subsonic craft 
of today (3). Nuclear ramjets are currently under study and 
should have a place in a second generation of supersonic 
terrestrial transports. However, as long as our petroleum re- 
serves hold up, the kerosene-fueled ramjet probably will be 
prevalent, and beyond that, it is likely that hydrogen or some 
other chemical fuel will be burned with air to supply the 
energy for supersonic and hypersonic ramjets. 


Scope and Background > 


Interpreted literally, the title of this paper might be ex- 

ected to restrict the following review and discussion to recent 

ork directly applicable to current ramjet engines. However, 
a broader viewpoint is taken for the following reasons. 

On the one hand, there have been no significant changes in 
the ramjet engineer’s understanding and application of com- 
bustion principles in engine design in the last few years, but 
rather a general refinement and drawing together of loose 
ends in the scientist’s picture of how combustible mixtures are 
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author is project supervisor for hydrocarbon combustor development and research on ad- 
vanced ramjet engines at Applied Physics Laboratory. After receiving his B.Ch.E. and M.S.E de- 
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was related to fundamental combustion research and particularly to flame velocity and stability. 
This research resulted in numerous papers and his Ph.D. thesis in 1953 for Case Institute of Tech- 
nology in chemical engineering. In 1954 Dr. Dugger joined the staff of the International Minerals 
and Chemical Corp. in Florida, where he worked until coming to APL in 1957. 


prepared, how heat and mass are transferred in a flame, and 
how—more empirically here—flames are stabilized in high 
velocity (subsonic) streams. The basic concepts and engi- 
neering practice are much the same as they were a few years 
ago, and the actual design and production of a particular ram- 
jet engine to fill a particular mission is still a matter of engi- 
neering judgment and compromise. Furthermore, design and 
combustion performance details on recently developed ramjet 
engines, such as those for the current and advanced versions 
of the surface-to-air missile Talos (developed by the Applied 
Physics Laboratory of the Johns Hopkins University for the 
U.S. Navy and produced by Bendix as prime contractor) and 
Bomare (produced by Boeing for the U. 8. Air Force, using 
Marquardt engines), fall within the realm of classified in- 
formation, so that they cannot be specifically disclosed here. 

On the other hand, many papers have appeared in recent 
years which are not related to current engines but do foretell 
future ramjet applications. Of particular interest are reports 
of successful experiments with stabilized detonation waves 
and with external combustion and theoretical papers which 
show the potential benefits of supersonic combustion for hy- 
personic airbreathing vehicles. In order to determine whether 
such unconventional combustion schemes promise advances 
in ramjet propulsion, the combustion process must be con- 
sidered in context as part of a propulsion system. The scope 
of the present paper is therefore extended to include a review 
and discussion of experiments and analyses related to uncon- 
ventional ramjet systems which are likely to find application 
in airbreathing vehicles within the next decade. 

For those readers who have not followed the ramjet litera- 
ture closely, some of the survey papers and documents which 
have appeared during recent years should be mentioned 
briefly as providing a background for subsequent sections of 
the present paper. 

Avery (5) succinctly reviewed ramjet history and progress 
up to 1955 and described the basic operating principles, rela- 
tive performance and advantages of ramjets to Mach 4, and 
current and foreseen applications of ramjets as of that date. 
He updated this picture and extended his considerations into 
the hypersonic regime in classified papers in 1958 and 1959 
(unclassified titles listed as (6 and 7)). Hawthorne (8) wrote 
a well-referenced review of the aerodynamic problems of air- 
craft engines, which includes the aerodynamic aspects of com- 
bustion, in 1957. Morris (9) briefly outined the 1958 ramjet 
state-of-art. Emmons (10) reviewed the nature and problems 
of combustion as an aeronautical science. Brief general re- 
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and future prospects appeared in Asrronautics, April 1959 
(11). Weber (4) discussed several unconventional ramjet 
concepts, including supersonic combustion, external burning, 
jet lift, fuel-rich operation, oxygen recombination and nuclear 
power. 

Leduc (12), Daum (13), and Roy (14) have reviewed the 
notable French accomplishments in jet propulsion. Par- 
ticularly interesting are their demonstrations of the use of 
ramjets in manned aircraft. Leduc’s first ramjet powered 
airplane was completely designed in 1939! Its combustion 
chamber was huge, and its overall efficiency was quite low, 
but its little-acclaimed first flight in 1949 is likely to take on 
greater historial importance in future years when ramjets 
finally achieve their deserved place in powerplants for manned 
aircraft. The Daum paper covers in some detail the major 
design features and accomplishments of the Nord Griffon air- 
planes, which have so ably demonstrated the advantages of a 
combined turbojet-ramjet (or turboramjet) powerplant for 
supersonic aircraft. The combined engine consists of a turbo- 
jet in parallel with a ramjet, with a common inlet and with air 
flow division phased gradually from pure turbojet at takeoff 
to pure ramjet at Mach 2. The great advantage of paralleling 
the two engines between common inlet and exhaust lies in the 
way their air flow and thrust coefficient vs. Mach number 
characteristics complement each other at supersonic speeds; 
the ramjet’s air flow and thrust coefficient increase as the 
turbojet’s decrease, so that the combined thrust coefficient re- 
mains nearly constant at supersonic speeds. An important 
corollary advantage of this characteristic is that less geometry 
variation in inlet and exhaust is required for optimum per- 
formance of the combination than would be required for either 
alone. Thus an engine comprising a small, relatively inex- 
pensive and easily serviced turbojet inside a lightweight ram- 
jet can produce more thrust per pound with less sensitivity to 
flight speed and to angle of attack or yaw than can the best 
and most complicated turbojets. The speed limitation of a 
turboramjet will depend upon the thermal protection system 
provided for the “swallowed” turbojet; Mach 4 appears to be 
a very reasonable goal, and a Mach 5 turboramjet would 
appear to be more practical than the Mach 5 turbojets with 
afterburners which have been proposed (1). Roy (14) con- 
siders conventional ramjets to be practical to Mach 6 and is a 
proponent of the detonative ramjet. 

Probert (15) and Jamison (16,17) presented fine reviews of 
British ramjet progress, ramjet applications in missiles, such 
as the Bloodhound, and expected future performance trends 
and applications. Jamison strongly supports the turboramjet 
for supersonic (around Mach 8) transport aircraft. 

Among the books on combustion and jet propulsion, the 
following works may be cited. A 259-page report entitled 
“Basic Considerations in the Combustion of Hydrocarbon 
Fuels with Air” by the NACA Lewis staff (18), quite thor- 
oughly covers that subject through 1955. Zucrow revised his 
text on jet propulsion with a two-volume series published in 
1958 (19), and includes chapters on ramjets (20); both of these 
are good as introductory textbook references on ramjet engine 
principles, but neither presents much information on ramjet 
combustion nor seriously considers speeds beyond Mach 4. 
A Russian book on ramjets appeared in 1958 (21); although 
it quotes heavily from American work, it is an impressive 
volume, and its publication indicates their strong interest and 
technological capabilities in the field. A book on ramjet 
technology by the staff of the Applied Physics Laboratory and 

various contributors under the Bumblebee program is to be 
published in 1960 (22). The books resulting from the bi- 
annual Combustion Symposia, e.g., (23,24), contain many 
valuable papers on combustion, flames and flame stabilization 
which bear on ramjet combustion problems, as do the proceed- 
ings of the AGARD Combustion and Propulsion Panel Meet- 
ings, e.g., (25,26). Other valuable texts on combustion which 
have appeared in recent years or will appear shortly include 


Basic Information on and Interpretation of 
Combustion Processes 


In the jet engines of today, heat addition is accomplished in 
the combustor by a number of overlapping physical and 
chemical steps, which might broadly be called atomization 
(and initial distribution) of fuel, vaporization of fuel, mixing 
of fuel and air, stabilization of flame through provision of re- 
circulation zones from which hot, burned gas can pilot the 
main flame, flame spreading, mixing of hot gas with cooling 
air (if any), and reassociation of dissociatied species as the gas 
cools by expansion in a supersonic nozzle. As combustor inlet 
conditions change due to flight schedule, or combustor h:rd- 
ware is changed to best suit the requirements of a partic ilar 
engine for a particular vehicle’s mission or missions, the rela- 
tive importance of the various steps listed above may s)iift. 
Breitwieser (11) gives a good overall exposition of the nature, 
interactions and changing roles of these processes in rarijet 
combustors, but a few specific references and comments are in 
order here. 

The literature on atomization and evaporation and com})us- 
tion of droplets and sprays is growing rapidly (see papers in 
(23) and reference lists therein). Weiss and Worsham’s 31) 
results on atomization in high velocity streams show that for 
simple tubes or pressure-atomizing nozzles the relative velocity 
between fuel and air has an overwhelming effect on drop size, 
which masks lesser effects of nozzle size, injection pressure 
and air stream density. Most fuel injectors are oriented 
contrastream for this reason. Air-atomizing nozzles are more 
effective than simple tubes or pressure types, but at relative 
velocities approaching 1000 fps any nozzle will produce a 
mean drop diameter of less than 20 uw (petroleum wax at 300 
F). Agoston et al. (32) and others (33) show that data on the 
evaporation of drops are usually correlated by semi-empirical 
equations, of the type introduced by Fréssling, as follows. 

For evaporation with relative gas motion in the absence of 
combustion or large temperature differences, diffusion controls 
the rate 

dm dD? 

dt 
Foster and Ingebo (34) found that when temperature dif- 
ference between drop and air was large, spray evaporation 
ate varied with AT®-%, 

For the evaporation of burning droplets (or droplet burning 
rate), heat transfer controls the rate, and the above equation 
is altered simply by the replacement of Sc’/* by Pr'/*, Bolt 
and Saad (35) measured burning rates of freely falling drops 
of hydrocarbon at 1500 F and found that most of the heat 
transfer required for evaporation had to come from conduc- 
tion and convection rather than radiation. Since both 
evaporation ahead of the burning region and droplet burning 
occur in ramjets, and these are mixed up with fuel-air mixing 
by turbulent transport as well as by molecular (diffusional) 
transport, the picture is indeed a complicated one. 

The physical picture of the flame stabilization process has 
been delineated rather clearly in recent years. There is 
general agreement that recirculation zones are established; 
in “baffle” combustors, these zones are in the wakes behind 
the baffles, whereas in “can” combustors, they fall behind the 
metal which separates the air (or mixture) admission holes. 
The wake behind a baffle is somewhat wider than the baffle 
and several times as long. Hot, nearly completely burned 
gases (average wake temperature is about 90 per cent of 
adiabatic flame temperature (36)), recirculate in regular 
eddies, which are parallel to the burned gas at the wake 
boundary and opposite in direction in the wake core. Turbu- 
lent transfer of heat and mass occurs between fresh gas and 
wake at the boundary. The critical parameter for flame 
stabilization can usually be explained in terms of reaction 
kinetics of approximately second order. See (36) in light of 
(37), or (38), 


= —C(1+ 0.3 Reé/Sc'/?) 
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The ignition delay concept of flame stabilization proposed 
by Marble and Zukoski (36) may be expressed quite simply: 
If the fresh gas, moving at velocity V2 past the transfer region 
at the boundary, does not reach its effective ignition tempera- 
ture by the time it reaches the end of the recirculation zone of 
length L, then the main body of the flame will blow out 


¢Vo/L > 1, flame stable 


5 TV2/L < 1, flame blows out 


The delay time 7 required to reach the effective ignition tem- 
perature is a chemical factor which is a function only of 
chemical properties of the mixture, i.e., varies with fuel type, 
oxygen concentration and fuel-air equivalence ratio; has 
Arrhenius-type temperature dependence and varies inversely 
with pressure (37,39). The ignition delay concept can be 
used to predict lean screech limits (large amplitude, stable 
transverse oscillations) in combustion ducts (36); it is found 
that ra/h is constant in a given system, where a is acoustic 
velocity and h is transverse duct dimension. The wake length 
and the relationship between V2 and approach velocity V, are 
charicteristic of a given flameholder-duct system. Wright 
40) finds that for flat-plate flameholders 


where A is duet height and d is flameholder width. This 
formula predicts that maximum blowoff velocity will be 
found for d/h = C;/C2 = 0.35 for flat plates or 0.56 for cylin- 
ders. 

Weiss, Rohrer and Longwell (38) present evidence that 
reaction rate is the controlling factor, based on an excellent 
correlation of lean blowout equivalence ratios for various 
fuels and flameholders with those for spherical reactor blow- 
out at a loading factor equal to that for lean blowout of iso- 
octane with the same flameholders. (They did not find as 
good correlation with an ignition delay criterion using Mullins’ 
data (37) or with a laminar flame velocity criterion based on 
Spalding’s (41) proposal that Vio ~ dpU,*.) The spherical 
reactor experiment was invented by Longwell (42) following 
papers by Longwell et al. (48) and Avery and Hart (44) de- 
scribing flame stabilization and combustor heat release rates 
for cases limited by chemical kinetics of near homogeneous 
reacting mixtures. Most of the results indicate overall or 
“global”’ reaction kinetics of approximately 1.8 order, so 
that for a given fuel and temperature the loading factor, in 
terms of air weight flow rate divided by reactor volume and by 
pressure to the 1.8 power, is a function of fuel-air ratio. For 
hydrocarbons, a “global activation energy” of 42 to 45 keal/g- 
mole is indicated (45). Other interesting results (46) from 
spherical reactors support the premise that reaction rate 
controls, but suggest that the pressure exponent may vary 
with fuel type or with feed injection momentum (mixing 
effects are not completely eliminated). Williams (47) points 
out that Khramtsov’s (48) turbulent flame results show that 
the feed rate per unit volume of turbulent flame zone is 
also proportional to p!-8. 

Weiss et al. (38) also studied the effect of heat loss from the 
recirculation zone to the wall of a rectangular can-type com- 
bustor on blowout loading factor for the combustor. In one 
case an independently variable, annular cooling air flow sur- 
rounded a metal can; in the other, the can walls were made of 
insulating firebrick. In the metal can, heat losses averaged 20 
per cent of the heat liberation rate; in the latter, only 6 per 
cent. Such a decrease in heat loss (14 per cent difference) 
would be expected to reduce blowout loading factor for a given 
fuel-air ratio by a factor of 4, but the observed factor was 
only 2. The authors concluded that in a jet engine applica- 
tion it would be better to make use of available space to in- 
crease recirculation zone volume rather than to insulate a 
smaller recirculation zone volume. 

Hottel et al. (49) have studied rod-stabilized flames in 
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volatile fuel sprays with the flameholder wetted to varying 
degrees by fuel. For wetted flameholders, part of the fuel in 
the recirculation zone comes from vaporization of the fuel on 
the flameholder, so that heat transfer between the recircula- 
tion zone and the wetted flameholder y lays an important part 
in the establishment of stability limits. 

Various studies of the propagation of confined (stabilized 
inside ducts) and unconfined (above Bunsen burners) turbu- 
lent flames through prevaporized, premixed gases have shed 
considerable light on the part of laminar flame velocity (most 
turbulent flame velocity correlations involve U,), intensity of 
approach stream turbulence (50,51), shear regions between 
burned and unburned gases and flame generated turbulence 
(51,52,53), pressure (48) and temperature effects (related to 
temperature effects on laminar flame velocity and approach 
stream Reynolds number (18)), and acoustic disturbances 
(53). Blackshear’s tests (53) indicated that significant in- 
creases in turbulent flame velocity could be obtained with 
the introduction of transverse velocity gradients by means of 
“benign” acoustic disturbances (92 db, serious damage to 
hardware occurs at 200 db). This corresponds to the frequent 
observation that combustion efficiency goes up when a ramjet 
“sereeches.”’ Blackshear suggested that improvements in en- 
gine combustion efficiency might result from putting whistles 
on the flameholders or making the combustor wall wavy. 

Basic information on heat and mass transfer in flames and 
on transport properties at high temperatures is still of im- 
portance in ramjet combustion because turbulent flame 
velocity—or flame spreading rate in a combustor—is related 
to laminar flame velocity (or at least to the way in which fresh 
mixture is burned at a local interface or in a pocket or an 
eddy) and to the more basic thermal and molecular diffusion 
rates as well as the momentum transfer rate. Fristrom and 
Westenberg (54) have done some excellent work in unraveling 
the mysteries of laminar flame front structure. From recent 
measurements on a low pressure methane-oxygen flame, 
Fristrom et al. (55) hope to get at the high temperature reac- 
tion kinetics. Such experimental information on reaction 
kinetics at flame temperatures would permit a much more 
accurate assessment of the problem of dissociation in combus- 
tion chambers and degree of recombination to be expected 
in supersonic nozzles, which is of major importance to ramjet 
(or rocket) thrust obtained with the high combustion tempera- 
tures associated with hypersonic flight. However, the reac- 
tion rates could not be determined from flame front structure 
experiments without knowledge of transport properties at 
high temperatures; the work of Walker and Westenberg (56) 
makes a valuable contribution to this field. Transport proper- 
ties may also be important in detonations (57). 


Advantages and Disadvantages of Higher 7 
Temperatures 

The variable which has the greatest impact on both the 
ramjet heat addition process and the overall ramjet design is 
temperature. The inlet air temperature to a conventional 
ramjet combustor is some large fraction of the inlet total tem- 
perature, which is a function of ambient temperature, flight 
Mach number and effective specific heat ratio 


7, — 1 


The combustor inlet air temperature which might be encoun-_ 
tered by a typical engine on a typical flight program is illus- 
trated by curve A in Fig. 1. The external skin temperature 
for the same aircraft, for areas 1 ft or more from leading edges, _ 
might follow curve B. For a ramjet burning kerosene, the _ 
temperature of the combustion products might fall between 
curves C and C’, depending on whether the engine follows a 
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Fig. 1 Temperatures which may be typical in ramjet engines at 

supersonic and hypersonic speeds. (Altitude assumed to in- 

crease with Mach number to maintain constant free stream dy- 

namic pressure of 350 pfsa) 


maximum thrust fuel schedule (near stoichiometric fuel-air 
ratio) or a maximum range fuel schedule (optimized for 
minimum specific fuel consumption; fuel-air ratio increases 
with speed until stoichiometric is reached near Mach 7).? Ap- 
proximate melting points of a few structural materials are 
shown for reference. 

Let us consider first the trends in conventional ramjet com- 
bustor design brought about by the increasing inlet air tem- 
perature with flight speed. The ramjets (and turbojets) in 
use today operate at speeds below Mach 3 and encounter com- 
bustor inlet air temperatures in the range 200 to 800 F. In 
order to obtain high combustion efficiency, the fuel must be 
injected some distance upstream of the flameholding point 
to allow time for vaporization and mixing, and a sheltered 
region must be provided for flameholding, as discussed in the 
previous section. Such flameholders have appreciable drag 
(see for example, Wright’s results (40), which are quite con- 
servative with respect to full-scale engine practice). Flame 
spreading rates are moderate, so that a tailpipe length of 2 
to 4 ft must be provided between flameholder and exhaust 
nozzle. 

Engines now under development for speeds of Mach 3 to 5 
are entering the spontaneous ignition regime. The flame- 
holder can be moved closer to the fuel injector, because 
vaporization rates are considerably higher. Less blockage or 
burner drag is required to achieve the same flame stability, 
and less tailpipe length is required for efficient combustion, 
because both blowoff velocity and flame spreading rate in- 
crease as power functions of inlet temperature. The maxi- 
mum space heating rate for a lean hydrocarbon-air mixture in a 
well-stirred reactor increases by an order of magnitude when 
the inlet temperature is doubled (45). 

Some possibilities for the virtual elimination of burner 
hardware drag in higher speed ramjets are seen in aerody- 
namic flameholding and flame spreading techniques and re- 

? Both the stoichiometric curve C and the minimum SFC curve 
C’ are estimated for an aircraft following a moderate, constant 
dynamic pressure trajectory of q~ = 350 psfa, with a diffuser 
kinetic energy efficiency of 0.94, giving combustion chamber _ 
sures of the order of 1 to 10 atm. 7 
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cessed geometry experiments. Various experiments have 
shown that flames in small ducts can be stabilized by opposed 
jets (58) or transverse fluid sheets (59), even at mild inlet 
temperatures, and it has been shown that stabilization stil] 
occurs through the formation of recirculation zones. Lertin 
and Salmon (60) described “aerodynamic combustion” 
studies culminated by flying-test-bed tests in a full-scale turbo- 
jet afterburner (45 per cent thrust augmentation on 6600-lb 
thrust turbojet). Their basic concept is to speed u)) the 
flame spreading process by using transverse air jets within 
the recirculation zone behind a gutter-type flamehc der; 
the air jets force the rapid turbulent spreading of avtive 
centers from the recirculation zone into the fresh mixture. 
(In the full-scale tests, they also injected a small fraction « f the 
total fuel behind the gutter.) Their preliminary tests in a 7-in. 
pipe had shown that tailpipe length required to produce 85 
per cent combustion efficiency could be cut in half by ble: ding 
2} per cent of the inlet air through the jet flame spre der, 
Such aerodynamic flameholders and/or flame spreaders may 
find use in higher speed ramjets, because stabilization and 
rapid flame spreading will require only very small (if any) 
piloting or recirculation zones. Fluid jets or sheets could ssue 
from the burner walls, from an otherwise required str:am- 
lined centerbody, such as a diffuser plug, or from the vicinity 
of smaller than usual flameholders. Or, the flame could be 
stabilized by the recessed geometry technique (61), in which 
a recirculation zone is established in a cavity in the engine 
wall (e.g., at the base of the ramjet diffuser centerbody) rather 
than by placing an otherwise unnecessary drag-producing 
obstruction in the main stream. 

However, the simplification of the conventional ramjet 
combustor hardware requirements from the performance 
standpoint will be compromised by two major problems as 
flight speeds enter the hypersonic regime, say, speeds above 
Mach 5. One problem is the problem of containing the high 
temperature gases without seriously compromising engine 
weight or performance, and the other is the loss in temperature 
rise attainable in the combustion chamber. In current ram- 
jets, the tailpipe and exhaust nozzle are often protected by 
passing a film of the relatively cool combustor inlet air between 
the tailpipe and an interior metal liner or “cooling shroud.” 
At hypersonic flight speeds this will no longer be feasible, be- 
cause the temperature of the combustor inlet air will be too 
high. Some small gains can be made by internal insulation 
of the tailpipe and nozzle with ceramics (62), or by the use of 
“refractory metals,” such as molybdenum or columbium al- 
loys, but it is obvious that elaborate, forced cooling systems 
would soon be required. Further discussion of the structural 
problems associated with hypersonic conventional ramjets is 
beyond the scope of this paper, but unconventional combus- 
tion systems which will relieve this problem are discussed in 
later sections. 

The effect of flight speed on the temperature rise due to 
combustion is seen by comparing either curve C or C’ with 
curve A in Fig. 1. As the flight speed and inlet air tempera- 
ture increase, less temperature rise is obtained for a given 
energy addition, because more and more of the energy is being 
absorbed by dissociation of the combustion products. That 
is, the static enthalpy of the ‘“‘burned”’ gas represents a larger 
fraction of the stagnation enthalpy, so that there is a smaller 
conversion of potential chemical energy to kinetic energy 1 
the combustor. As noted previously, part of the dissociation 
energy will be recovered in the exit nozzle. Moreover, evenif 
there were no further change in the chemical composition of 
the burned gas leaving the combustor, part of the energy ab- 
sorbed by dissociation would still be useful, because for the 
partially dissociated gas both the specific heat ratio and the 
molecular weight are lower, and a reduction in either of these 
increases air specific impulse and exit stream thrust obtained 
with a given burned gas temperature. Nevertheless, the loss 
in engine thrust due to dissociation will be felt more and more 
as flight speed increases in the hypersonic range. 
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Ways to Circumvent the Effect of T,,, on AT, 


The reduction in combustor temperature rise with increased 
flight speed can be mitigated by using fuels which form stable 
products of combustion, for example, fuels containing metals 
which form stable oxides not subject to dissociation. A further 
gain will be made if the stable products appear in a condensed 
phase, that is, as small droplets or fine powderlike particles, 
so that an exothermic phase change releases more heat in 
the combustor; again, oxides of metals provide examples. 
Another type of tailored fuel would be one which would disso- 
ciate or react to give a burned gas of lower molecular weight 
and hence higher specific thrust for a given temperature rise. 

On the subject of fuels which form metallic oxides, Olson 
and Setze (63) discuss the advantages and problems of the use 
of high energy fuels (borohydrides or organoborohydrides) 
and fuels containing magnesium or aluminum. The high 
enerzy fuels have higher heats of combustion and higher flame 
temperatures than hydrocarbon fuels. This means that for 
a given thrust level, a lower fuel-air ratio is required, so that 
for given fuel weight, the aircraft using high energy fuel will 
have a considerably greater range. The familiar Breguet 
range equation may be written (for moderate speeds) 


L 1 


er" > 
where 
h. = heat of combustion of the fuel 
n = overall engine efficiency 
L/D = lift to drag ratio 
Wy, = weight of fuel consumed in cruise 
Wy = gross weight at beginning of cruise 


It is seen that h, first enters directly, but, even more important 
for very long ranges, the lower fuel consumption rate per- 
mitted with higher h, keeps the gross weight from increasing 
as quickly; that is, it puts a restraining leash on “longer range 
means more fuel means more structure means more fuel, 
ete.’ An example for aircraft carrying 45,000-lb useful pay- 
loads is cited: For 6000-mile range the plane using kerosene 
(JP-4) would weigh 450,000 Ib compared to 212,000 Ib for a 
similar but smaller plane using a high energy fuel with a 40 
per cent higher he. 

Aside from the heat of combustion advantage, the borohy- 
dride fuels are attractive due to their very high flame veloci- 
ties and very wide stability limits. The maximum laminar 
flame velocity of pentaborane is 22 fps, compared to 9 fps for 
hydrogen measured by the same method at the same condi- 
tions (64), and less than 2 fps for propane. As noted in an 
earlier section, flame spreading rate is related to laminar flame 
velocity Us, and space heating rate (loading factor) is roughly 
correlated with U7; it becomes apparent then that much 
smaller recirculation zones are required for flame stabilization 
with the borohydride fuel, and that good combustion efficien- 
cies could be achieved in much shorter combustion chambers. 
Berl and Renich (65) review the available unclassified infor- 
mation on the combustion of boron hydrides through July 1958. 
Breisacher et al. (66) report studies cn the structure of laminar 
flames of borohydrides and borohydride-hydrocarbon mix- 
tures. In the flames of the borohydride-hydrocarbon mixed 
fuels, the borohydride is always oxidized first; the hydrocar- 
bon may only be cracked and partially oxidized if the borohy- 
dride flame is lean, or it may burn in a second luminous reac- 
tion zone following a hotter (nearer stoichiometric) borohy- 
dride reaction zone. The wide flammability limits'and low 
ignition temperatures of these fuels—as well as of metal 
alkyls, such as aluminum trimethyl, boron triethyl, etc. (67, 
68), and mixed metal hydrides, such as aluminum borohydride 
(6S)—suggest that they could also be used as igniter or pilot- 
ing fuels to widen the operating limits of jet engines using a 
hydrocarbon fuel as the primary energy source. (External 
combustion of such fuels in supersonic air streams is discussed 
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later.) However, these same combustion properties, plus 
toxicity, are disadvantages from the safety standpoint. 

It should be pointed out that the heat of condensation of 
B.O; usually represents a large part of the heat of combustion 
advantage of a borohydride fuel over a hydrocarbon fuel. 
This means that if temperatures are too high to obtain the 
condensed phase in the combustor, one must again rely on the 
lower temperatures in the supersonic nozzle to bring about the 
condensation before the vapor leaves the region in which it 
can influence the pressure field of the engine. If condensation 
does not occur within the pressure field of the engine, as much 
as two thirds of the high energy fuel’s heat advantage may be 
lost, depending on its boron and hydrogen contents (63). Of 
course, condensation offers problems too; if BxO3 condenses on 
relatively cool engine hardware, it forms unwanted deposits. 
If particle size is too large, the heat of condensation may not be 
effectively transferred to the gas. Olson and Setze (63) de- 
scribed studies of particle growth and deposition. 

The metals magnesium and aluminum do not have high 
heats of combustion, but due to their high densities (com- 
pared to liquids) and high stoichiometric fuel-air ratios, they 
do look attractive on the basis of short range, high thrust 
applications, because they can add a lot of heat per pound of 
inlet air. For this reason, the burning properties of the metals, 
and particularly the properties of slurries of magnesium in 
JP-4 fuel, have been studied (63). A 50 per cent slurry of fine 
Mg particles (1u) in JP-4 gave a combustion efficiency one 
third higher (0.8* vs. 0.6) than JP-4 alone in a 19-in. long tail- 
pipe. 

As noted above, at high flight speeds and high combustion 
temperatures, the products of “combustion” of hydrocarbon 
fuels will be highly dissociated in the combustion chamber 
and must recombine in the exit nozzle if all of the potential 
chemical energy is to be converted to kinetic energy; this also 
applies to condensation of metallic oxides in the case of high 
energy fuels. Little is known about high temperature kinetics 
and recombination rates, even at fixed conditions; in super- 
sonic nozzles the problem is compounded by rapid changes in 
pressure, temperature and velocity, which in turn depend 
both on the changing properties (y and molecular weight) of 
the gas mixture and on the design of the nozzle. Unfor- 
tunately, pressure decreases much more rapidly than tem- 
perature, and velocity is increasing at the same time. When 
more is learned about the kinetics, it may be possible to tailor 
the nozzle (perhaps by adding more slowly diverging sections 
at points where favorable changes would occur) to obtain 
higher efficiency. Meanwhile, some very interesting work is 
noted in which the situation has been turned around—a di- 
verging reactor, or supersonic nozzle, may be used to deter- 
mine the kinetics of a simple reaction—work such as Wegener 
(69) has done for 2 NO, = N2Ox. 


Use of Unconventional Engines for Hypersonic 
Speeds 


An interesting proposal for a ramjet capable of flight speeds 
from supersonic all the way to Mach 18 or 20 has been ad- 
vanced by Breitwieser and Morris (see (4)). Their proposed 
“fuel-rich ramjet,”’ is not limited to stoichiometric fuel flows, 
but is deliberately operated at higher and higher equivalence 
ratios as flight speed increases. Although fuel specific impulse 
falls rapidly as equivalence ratio is increased, it remains higher 
than rocket specific impulses, and an adequate accelerative 
thrust margin is provided to speeds of the order of Mach 20. 
Advantages cited for such a system are: 

1 Rich operation means lower burned gas temperatures 
and less dissociation. 

2 Varying fuel rate can be used to match inlet and exit 
flows in a simple fixed geometry system. 

3 The large fuel flows can be used for cooling the engine 
and possibly the airframe as well. 

Another approach to higher propulsive efficiencies at very 
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high speeds is to diffuse the inlet air only to some moderate 
supersonic velocity before adding the heat (2,4,70). Such an 
engine might be called a supersonic combustion ramjet or 
SCRJ (Fig. 2). Heat addition to a supersonic stream results 
in a decrease in Mach number. The ideal supersonic engine 
performance would be obtained by adding just enough heat to 
reduce the Mach number to unity and then expanding the gas 
supersonically without ever taking a strong shock loss. In 
other words, the desired heat addition level would determine 
the desired combustor inlet Mach number. For any given 
heat addition level at any given flight speed, the highest static 
temperature reached in the supersonic combustion engine 
would be 10 to 12 per cent lower than that obtained with the 
same heat addition in subsonic flow. This 10 to 12 per cent 
reduction in static temperature is very important insofar as 
reducing dissociation is concerned. The higher total pressure 
loss due to fuel-air mixing and heat addition in supersonic 
rather than subsonic flow would be compensated (to some de- 
gree dependent on M ,.) by a lower total pressure loss in the 
diffuser. Analyses suggest that the gain due to lesser dissocia- 
tion (and at high /,, lesser overall total pressure loss) may 
make this engine superior to a conventional engine at Mach 
numbers higher than 7 or 8 (2,4,70). Thrust augmentation of 
ramjet engines by afterburning in the supersonic nozzle may 
prove to be another valuable application of supersonic com- 
bustion, in that it would provide a means for increasing ramjet 
thrust without affecting the diffuser performance (4). 

A modification to the latter approach is the detonative or 


M2 ~ 2.5 to3 


Mg>2! 


EXTERNAL = INTERNAL HEAT 
DIFFUSION DIFFUSION ADDITION 
Fig. 2 The supersonic combustion ramjet (SCRJ) as discussed in 
(70) and (2). The reduction in static temperatures in diffuser 
and combustion chamber, as compared to temperatures which 
would be experienced in a conventional ramjet (CRJ), means less 
dissociation in the burned gas 


EXPANSION 


Fig. 3 The standing wave (detonative) ramjet (SWRJ), with an 
oblique detonation wave stabilized by a streamlined centerbody, 
as proposed in (71) 
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standing wave engine, which might be designated SWRJ. 
See Fig. 3, after (71). Fuel would be added to a supersonic 
stream ahead of a strong shock wave held stationary by the 
duct geometry with the heat release occurring behind the 
wave. Roy (72) proposed such an engine in 1946. Nicholls 
et al. (73) and Gross et al. (57,74,75) have demonstrated 
stabilized detonation waves in hydrogen-air streams at Mach 
numbers in the range 2 to 5. Most of Gross’ fascinating ex- 
periments were done with hydrogen addition upstream of a 
Mach-reflected normal shock. The latter was obtained by 
placing opposing wedges on the top and bottom walls of a 
two-dimensional, slightly divergent test section of a nominal 
Mach 3 tunnel (Fig. 4a). The geometry was chosen so that 
oblique waves from the wedges generated a normal wav:: at 
the center of the tunnel. Hydrogen, introduced through a 
small streamlined injector just upstream of the nozzle thrvat, 
diffused through the air as it was swept toward the test sec- 
tion. At air stagnation temperatures of approximately 1860 R 
and higher, the resulting mixture ignited behind the shock 
system and caused the plane wave to move forward and be- 
come larger in vertical dimension. Even though the ‘uel 
concentration profile was far from flat ahead of the shock, the 
shock itself remained rather flat. The extent of the loca‘ion 
and size change was directly related to the fuel injection ‘ate 
and inversely related to temperature. This was interpr:ted 
as suggesting that Q/c,7i, the dimensionless 
second parameter,’ where Q is heat added and 7, is static 
temperature before the shock, plays an important role in 
characterizing detonation waves. 

An ignition “hysteresis” effect was found (57,75): After ig- 
nition at high air stagnation temperature, the air stagnation 
temperature could be lowered to as little as 860 R (static 
temperature at Mach 3, 310 R, or —150 F) before ignition 
ceased, and the normal wave regressed to its ‘“no-burning” 
position. The extent of the hysteresis effect, which was also 
related to Q/cp7,, was purported to imply that transport 
mechanisms within the wave are important once the combus- 
tion has been initiated, thus supporting the theoretical work of 
Hirschfelder et al. (76), but denying the assumption of zero 
transport properties in the von Neumann-Zeldovich type of 
analysis which separates shock wave and chemical reaction. 
Supporting evidence that no ordinary “ignition delay” time 
was involved was obtained by superimposing Schlieren and 
sodium emission (from injected NaCl) photographs, which 
showed that combustion began in less than 5 microsec. By 
using a dual fuel injection system, exothermal waves of hy- 
drogen, methane and air were also stabilized. Pure methane- 
air waves were never successfully stabilized, but it was be- 
lieved that careful experimentation would reveal combinations 
of conditions conducive to their stabilization. Oblique deto- 
nation waves were also produced with hydrogen by using a 
single wedge. Ignition hysteresis was observed in all cases. 

Nicholls et al. (73) did their experiments with a small, 
underexpanded, free jet, which produced a Mach reflected 
shock in the expanding jet at a local Mach number of 3.5 to 
6.2 (Fig. 4b). Their experiments, chronologically ahead of 
Gross, indicated an ignition delay of the order of 25 microsec 
with a stagnation temperature of ~2150 R, an M, of ~5, 
and roughly double the stagnation pressure Gross used. Con- 
trary to Gross, Nicholls contends that the separation of shock 
and combustion, and the corresponding ignition delay are real; 
perhaps at some later date an explanation will be found which 
ties together the apparently incompatible results from these 
different experimental techniques and experimental condi- 
tions. This same group also did an analysis (71) for a detona- 
tive engine as sketched in Fig. 3. They estimated the stagna- 
tion pressure loss for the mixing of the fuel with the supersonic 
air stream ahead of the oblique detonation wave. The esti- 
mated mixing loss was reasonable, 10 to 20 per cent, but it was 
based on the assumption that no oblique shocks were caused 
by the fuel injectors. This problem of getting the fuel into the 


supersonic stream without shock losses—which applies also 
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to the SCRJ previously discussed—is likely to provide the ae 


greatest obstacle in the course of demonstrating the calculated 
potential performance of detonative or supersonic combustion 
ramjets. 

Various estimates have been made for the performance 
potentials of these engines compared to conventional ramjets 
(2,4,14,57,70,71). Fig. 5 compares overall engine efficiencies 
estimated for SCRJ and SWRJ with those for a conventional 
ramjet (CRJ) and a turbojet. Ramjets are superior to turbo- 
jets for flight speeds above Mach 3 or 3.5. The conventional 
ramjet probably will reach its maximum efficiency at some 
speed between Mach 6 and Mach 8, depending upon the ef- 
ficiencies of the various engine components. The supersonic 
combustion ramjet may surpass the conventional ramjet for 
specds greater than the CRJ’s maximum efficiency speed, and 
the detonative or standing wave ramjet will follow slightly 
behind the SCRJ due to its shock loss. The curves for Fig. 
5, except for the “optimistic” CRJ curve, are taken from (2); 
the CRJ curve is based on a higher diffuser kinetic energy 
efficiency (0.94 compared to 0.90) and higher combustor 
tot! pressure recovery (0.96 compared to 0.86). All curves 
are based on the use of kerosene fuel on an approximately opti- 
mui fuel schedule for an aircraft following a moderate, 
constant dynamic pressure trajectory (qo = 350 psfa). 
Dynamic equilibrium is assumed in all cases for the burned 
gas expansion—an optimistic assumption—but the favorable 
effect of centrifugal lift at high speed, which would serve to 
increase range, is not included in efficiency estimates. 


External Burning Experiments Point to 
Still Higher Speeds 


The previous discussion pointed out ways to extend the per- 
formance barrier due to combustor temperature rise, but an 
equally important problem is that of the large structural cool- 
ing requirement imposed by the staggering temperatures 
encountered in the diffuser, combustor and exit nozzle at 
flight speeds above Mach 6. Either supersonic combustion 
ramjets or standing wave ramjets would have a smaller dif- 
fuser cooling requirement than a conventional ramjet due to 
the lower heat transfer rate per unit area with supersonic flow 
(70). The combustor and nozzle areas for the standing wave 
ramjet may be much smaller than those for a conventional 
ramjet (14,71) which may further reduce overall cooling re- 
quirement. An SCRJ would have a lower heat transfer rate 
per unit area in the heat addition region than a subsonic com- 
bustion engine (70), but more area would be required. (Lower 
temperature means lower reaction rate; coupled with higher 
velocity this means at least an order of magnitude increase in 
combustion chamber length, which in turn means more area 
for heat transfer.) It is difficult to estimate, therefore, 
whether the combustor cooling requirement would be greater 
or smaller. 

A better way to push back the structural temperature 
limitation on ramjet engine flight speed might be to turn the 
engine inside out, so that all of the compression-heat addi- 
tion-expansion cycle is carried out on external surfaces, which 
are free to radiate to the atmosphere. Thus the all-external 
engine (ERJ) is simply a supersonic wing beneath which 
mass and heat are added by injecting fuel into the supersonic 
air stream. The addition of mass and heat beneath a wing 
causes the supersonic air stream reaching the heat addition 
region to be deflected away from the wing, and increases lift 
in much the same way that an increased angle of attack would 
increase lift. If the heat or mass addition is concentrated 
under the rear part of the wing (Fig. 6, after Weber (4)), where 
the increase in surface pressure will serve to reduce drag, then 
the maximum lift to drag ratio may be increased by a factor 
of 2 or more (77,78). If the pressure rise under the rearward- 
facing surface is great enough, the “wing” will have a net 
thrust, and it becomes a “propulsive wing” or “engine.” 
Luidens and Flaherty (79) estimate that a simple wedge air- 
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tion waves, as reported in (75) and (73) 
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Fig. 6 Sketch of a ‘‘propulsive wing,’’ or external ramjet (ERJ), 
after (79) 


foil of small thickness to chord ratio and large rear to fore 
angle ratio could become more efficient than a conventional 
ramjet engine plus wings at speeds greater than about Mach 8. 
Since the major part of the burned gas expansion force goes 
into lift rather than thrust, lift to drag ratio and thrust specific 
impulse cannot be separated for an overall engine efficiency 
comparison, such as Fig. 5. The comparison of Luidens and 
Flaherty is made on the basis of a range parameter L/Q, 
where Q = wyh,-. For constant velocity cruise with thrust 
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equal to drag, this parameter, for a given h,, is proportional to 
I ,L/D. 

Dorsch et al. (78,80,81) have conducted some interesting 
experiments on external burning. Their results for the com- 
bustion of aluminum borohydride under a two-dimensional 
supersonic wing in Mach 2.5 and 3.0 air streams (78) appear 
to agree reasonably well with the “conservative” curve of 
Luidens and Flaherty (79) for a propulsive wing with con- 
centrated heat addition. Dorsch et al. (78) show that their 
results are also in reasonable agreement with the predictions 
of a linearized theory and a graphical method. This agree- 
ment between experiment and theory at low Mach numbers 
might be optimistically interpreted to mean that the predic- 
tion (79) of external ramjet superiority to conventional ramjet 
plus wing for /., > 8 may be realized, and that the same 
propulsive wing or ERJ will be superior in fuel economy to 
chemical rockets to 7,, > 12. Dorsch et al. have burned, or 
attempted to burn, several other pyrophoric fuels in a Mach 
2.5 air stream (81). Gross has also demonstrated burning in 
a supersonic boundary layer with hydrogen injection (75). 


Concluding Remarks 


Combustion science and ramjet engineering technology 
have reached a stage such that efficient ramjet engines can be 
designed and produced for vehicles with flight speeds up to the 
threshold of the hypersonic regime. The higher temperatures 
associated with hypersonic flight speeds will have a con- 
siderable impact on both combustor temperature rise and 
overall engine design, but various approaches are envisioned 
for overcoming or bypassing these obstacles. For conven- 
tional engines with subsonic combustion, it may be possible 
to tailor a fuel to mitigate the deleterious effect of low com- 
bustor temperature rise on thrust. The thrust obtained from 
any hypersonic ramjet engine will depend in large measure 
on the effectiveness of recombination reactions and/or favora- 
ble phase changes in the supersonic flow in the exhaust nozzle. 
With more information on the kinetics of these changes (which 
appears to be forthcoming), it may be possible to tailor the 
nozzle to obtain higher thrust. 

Supersonic (SCRJ), detonative (SWRJ), or external com- 
bustion (ERJ) engines may be employed for flight at Mach 
numbers above about 8. The crossover points obtained from 
performance estimates are of course strongly dependent on the 
assumptions used to make the estimates, but the trends indi- 
cated by comparisons such as Fig. 5 are probably realistic, 
provided that the promise of early experiments with the vari- 
ous combustion systems at low Mach numbers holds up to 
hypersonic Mach numbers. In any event, these recent com- 
bustion experiments and aerothermodynamic analyses have 
suggested several ways in which ramjet propulsion might be 
used to power hypersonic airbreathing vehicles in coming 


years. 

Nomenclature 

a = acoustic velocity, fps 

C = constant 

Cp = specific heat at constant pressure, Btu/lb-R 

D = drop diameter 

d = flameholder width (characteristic dimension), ft 

h = duct height, ft 

= heat of combustion for fuel, ft-lb/lb 

I; = fuel specific impulse, gross thrust divided by fuel rate, 
(2000/1, = SFC) vie 
L, = gth of recirculation zone, ft 

L/D = lift to drag ratio 
M,., = flight (free stream) Mach number 

Pr = Prandtl number i 

p = absolute pressure, atm or psfa 
Q = heat added, Btu/Ib_ 

Q = heat addition rate, Btu/sec 


dynamic pressure of free stream, yM 2/2, psfa 


= range, ft 
Re = Reynolds number 
Sc = Schmidt number 
Tac = static temperature of air at inlet to combustion, R 
AT, = total temperature rise in combustor, R 
T:.. = total (stagnation) temperature of free stream, R 
T, = static temperature before shock, R 
T.. = static temperature of free stream, R 5 dehy, 
t = time 
U, = laminar flame velocity, fps 
V = velocity, fps a | 
V,: = reference velocity upstream of flameholder, fps 
V. = velocity of fresh gas mixture past wake, fps ™ 
Wy, = weight of fuel consumed during cruise, lb 
W, = gross weight of aircraft at beginning of cruise, lb 
wy, = fuel flow rate, lb/sec 
q = effective ratio of specific heats during compression, |' 
ne = overall engine efficiency, ft-lb work/ft-lb input 
T = ignition delay 
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D. B. SPALDING! 


Imperial College 
London, England 


High Pressures 


Conventional theories of particle combustion become inaccurate when the gas pressure ries, 
and are completely useless above the critical pressure. The paper discusses the transient effects, 
neglected in conyentional theories, which dominate the burning process at high pressures. A 
quantitative theory is presented, valid for these conditions, of the time variations of flame radius, 


burning rate and combustion efficiency. 


N CONSIDERING the influence of droplet vaporization reaction kinetic considerations to be eliminated. The «p- 

on the performance of gasturbine, ramjet or rocket com- propriate boundary condition for a phase interface is «!so 
bustion chambers, it is customary to calculate the burning stated. 
time of a droplet and to compare this time with the residence Then the quasi-steady theory is reviewed. It is shown ti:at 
time of the droplet in the ee Various formulas are this is not even a first-order approximation to the transicnt 
available for the burning time (1,2,3)?; they have in common problem. A new theory is therefore derived which, when the 
the assumption that the process is quasi-steady, i.e., that the property variations are neglected, gives an explicit formula 
transient processes in the gas can be neglected.* The pre- for the burning time which is valid for particles at high pres- 
dicted burning times have been found to agree well, in the sures, when the gas state is not too near stoichiometric. 


main, with the results of laboratory-scale experiments on 
single droplets (2,3,6 to 10), most of which have been carried 


out at atmospheric pressure. q 
Despite the good service which the quasi-steady formulas JA ee 


have given, their validity and utility are doubtful at the in 
high pressures prevailing in current jet engines. There are eneral conservation and flux laws 


pressure of the injected fuel, the intent heat of vaporisation steady state. Other treatments of a similar nature may be 

decreases to zero; the quasi-steady formulas correspondingly found in (13 to 23). 

predict a vaporization rate which rises to infinity. Secondly, The conservation and flux laws of the chemical element 

spective of changes late > mass of gas 

and sn poctave of changes to the lat heat, the m ‘= of ga a@ in a single-phase fluid mixture are expressible through the 

in the region influenced by the droplet increases steadily with equation 


pressure; its heat capacity therefore also increases, rendering 
illegitimate the neglect of the transient term in the energy 
equation. Kumagai and Isoda (11) have pointed out that 
the observed ratio of flame diameter to droplet diameter, even 7 
at atmospheric pressure, suggests the influence of transient where 
processes. = 
The transient burning of a particle has been considered G 
by Summerfield (30) in connection with his theory of the 
burning of composite solid propellants. The present paper 


C(Vna) — ne,jDip (Vm;)] = 
j 


= vector mass velocity of mixture, lbm/ft? hr 
Ne,j = mass of chemical element @ per unit mass of chemic: cal 
compound j 


was stimulated by, and may be regarded as an extension of, sa i die 
Summerfield’s approach. 4 
= mass of chemical element @ per unit mass of local mix- 
ture 
Contents of the Present Paper m; = mass of chemical compound 7 per unit mass of local 
mixture 
The present paper attempts to predict the burning rates D; = local diffusion coefficient of chemical element 7 in mix- 
of fuel particles in a way which is valid at high pressures, ture, ft?/hr 


= local mixture density, lbm/ft? 


whether or not the phase boundary has completely disap- 
= time, hr 


p 
peared. t 
First we state the differential equations which express the >> = summation over all chemical compounds which are 
conservation and transport of matter and of energy when all j present 


the diffusivities are equal. The latter assumption permits : : eet 
The first law of thermodynamics and the Fourier heat con- 


Received Jan. 1959. duction law may be expressed, for a fluid in which thermal 
! Professor of Heat Transfer, Department of Mechanical Engi- radiation and the kinetic energy of the mean motion are 
i negligible, through the equation 
Numbers in parentheses indicate References at end of paper. 
3 Agafanova, — h and Paleev (4) claim that the originator a dh . Op 
of the theory was Varshavskii (5); the present author has not G(Vh) —V [>oA;Dio(Vm;)] —-VIKVT)]=-p—+-— PI 
had an opportunity to see the latter publication. j oat 
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h = p hym; 


= local specific enthalpy of mixture, Btu/lbm 


k = local thermal conductivity of mixture, Btu/ft hr F 
7 = local mixture temperature, F 
p = local pressure, Btu/ft* 


Simplified conservation and flux laws = 

If we assume that all the diffusion coefficients are equal at 
any particular point in the fluid, we may drop the suffix 
j to D;; then Equation [1] simplifies to 


G-(Vne) — V[Dp(Vna)] = 


One such differential equation can be written for each chemi- 
cal element. 

If we assume, in addition, that Dp = k/c, where c is equal 
to the local mixture specific heat at constant pressure, and if 
the variation of pressure with time can be neglected, Equation 
[2| reduces to 


— p(Ong/Ot) [3] 


G- (Vh) — V(DpAh) = —p(dh/dt) 


A mixing theorem 


quations [8 and 4] differ only in their respec tive ndent 
variables. If the boundary conditions of the two equations 
appropriate to a particular flow are also similar, the solutions 
of the equations are bound to be similar: nz and h must be 
linearly related. 

In particular, in a “two-stream” situation in which the 
whole of the fluid under consideration results from the mixing 
of two streams in a region with adiabatic impervious walls, 
it is easy to show that the linear relation between nq and h is 
precisely the same as that which holds for an adiabatic steady- 
flow mixing chamber. For concreteness, we will regard the 
two streams as consisting of pure fuel and of pure oxidant. 

The last conclusion permits us to ascribe to each mixture 
particle a “fuel fraction” f which is related to ng and to h by 

Na Na.o h ho [5] 


Na,F — Na,o hp — hp 


where suffixes O and F refer to the oxidant and fuel streams. 
The quantity f/(1 — f) is equal to the ratio of the mass rates 
of fuel and oxidant into a steady-flow adiabatic mixing 
chamber if the outgoing stream is to have the properties n_ and 
h. 

Equation [5] may be substituted in Equations [3] or [4] 
to yield 


G-(Vf) — = —p(af/at) 


This equation will be used extensively later. — 
The Mass Transfer Boundary Condition 


Equations [3, 4 and 6] hold for a region of uniform phase. 
as is shown in (1) 


At the boundaries of such a region we have, 


re 
where 
m"” = component of G normal to phase boundary, directed 
into fluid phase under consideration, lbm/ft*hr 
y = distance normal to phase boundary, in same direction, 
ft 
P = na,h, f or any fluid property obeying a similar equation 
S = state of the fluid under consideration immediately ad- 
jacent the phase boundary 
T = value of the property in the “transferred substance,”’ 
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h; = local partial enthalpy of chemical compound j, — 


i.e., in a fictitious stream which is entirely equivalent 
in its effect on the fluid under consideration to all the 
material and heat fluxes across the phase boundary 


_ When the quantities Ps and P; are uniform and constant 
for all parts of the phase boundary, it is convenient to follow 
(1) by defining a new variable b through 


b =P/(Ps — Pr) [8] 


The differential equation and boundary condition then take 
the form 


G-(Vb) — V[Dp(Vb)] = —p(db/at) [9] 


> 
and 
= (Dp)s(0b/dy)s [10] 


Equations [9 and 10] will be used extensively. We there- 
fore note already two particularly important forms of Equa- 
tion [8] valid for the case when the transferred substance is 
fuel: 

1 If the property P is the fuel fraction f, clearly fp = 1. 
There results 


b = f/(fs — 1) 


2 If we choose the zero of enthalpy so that unburned fuel, 
products and inerts have zero enthalpy in the gas phase at 
state S, and if, as is true for burning volatile fuels, no un- 
burned oxidant is present at the surface, it is possible to re- 
place hs — hz by Q, this quantity being the heat transfer from 
the gas to the surface per unit mass transfer (Btu/lbm). 
Q is, of course, often equal to the latent heat of vaporization. 
Equation [8] then takes the form res * @ 


=h/Q [12 


Quasi-Steady Theory 


The System Considered 


Fig. 1 shows a spherical particle, for example a fuel drop- 
let, in a large gas reservoir. The gas is at rest except for a 
radial motion occasioned by the mass transfer outward from 
the particle. Often there will be a spherical flame surround- 
ing the particle (at the radius where f = fetoicn, the stoichio- 
metric fuel-oxidant ratio), but we need not consider this at 
the moment. The assumption of spherical symmetry is 
known to be well-justified when the droplet is very small. 

It will be supposed that the rate of change of droplet radius 
is so slow (later the meaning of “slow” will be analyzed quan- 
titatively), that the 0/d¢ terms can be omitted from the 
equations. The process is therefore treated as quasi-steady. 
The first application of this assumption to droplet vaporization 
appears to have been made by Fuchs (12). 


4 


Fig. 1 Conditions adjacent to burning fuel droplet 
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and in addition we have 


_ r=rg: b = bg (gas in contact with the surface) 15] 
= ©: b = bg (gas in the reservoir) [16] 


Writing the partial differential equation of conservation [9] 
for this case in spherical coordinates, we obtain the ordinary 


differential equation 
ad db , ts2db If the gas can be taken as in thermodynamic equilibrium 
ao. ( - ;) aie at all points, Dp is a known function of 6, which may have 
, ne some such form as that of Fig. 2; the maximum value occurs 
oe as in the neighborhood of the flame where the temperature is 
ig highest. It is known that departures from equilibrium are 


r = radius 
radius of particle surface 
mass transfer rate at the surface 


always small if the flame stays alight (14, 24). 


Solution for Burning Rate 


Equation [13], with boundary conditions [14, 15 and 16], 
is easily solved. The burning rate is then obtained from 


mrs Dp db 
(Dp)s ’s (Dp)s (b — bs + 1) 


[17] 


which may be evaluated by quadrature; the integral repre- 
sents the area shaded in Fig. 2. 

It is often convenient to express [17] in terms of a mean 
value of Dp, namely (Dp)m, and the “transfer number” B, 
defined by 


B=be — bs 
Then Equation [17] becomes 
th"rs/(Dp)m = In (1 + B) [19] 
The mean value of Dp must therefore be chosen so that 


(Dp)m = Jos (b-bs 


As bs 
ak 

Se Evaluation of burning rate 
Fig. 2 Determination of burnin te when Dp d d 
a Equation [19] can be evaluated, in any particular case, 


A) “ 4 ae > : by insertion of the appropriate values for bg and bs, in which 
can be a dimensionless measure of either composition or 


[20] 


enthalpy. The difference B commonly has values of the 

order of 5 for liquid fuel combustion. Reference (25) con- 

: tains discussion of the expressions for B which are useful in 

particular circumstances. In certain cases the exact deter- 

mination of bs presents difficulties. Reference (26) describes 

2) a method of doing this graphically on an enthalpy-composition 
: 1 - chart; this matter will therefore not be discussed further here. 


- Solution for Burning Time 


TIME 


>. 


We now note that the transfer of mass from the particle 
to the gas phase causes a diminution in the particle radius 
in accordance with 


Fig. 3a Appearance of burning droplet: Quasi-steady theory 


—d m 
[21] 


dt Pp 


where p, is the particle density, which is assumed much 
_ greater than the gas density. 

If it can be assumed that B and (Dp) m are independent of 
time, it is an easy matter to combine Equations [19 and 
- 21), to integrate, and so to obtain the time for complete 


_ disappearance of the particle t. We find Cr AL 
ty = + B) 


where rs, is the particle radius at time zero. 
The variation of flame appearance with time may also be 
deduced from the equations. It is illustrated in Fig. 3a for 


comparison with the alternative predictions of the transient 
Fig. 3b Fuel fraction profiles: Quasi-steady theory eer (Fig. 5b). 
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Burning time of a liquid close to the critical pressure 


Equation [22] shows that the burning time of a particle of 
given size and density decreases to zero as B increases to in- 
finity. Now B can be written (hg — hs)/Q; and in many 
cases Q is equal to the latent heat of vaporization of the 
fuel. As the critical pressure is approached, the latent heat 
tends to zero. So the burning time tends to zero also. Com- 
forting though this conclusion is, it must be subjected to ex- 
amination before acceptance. 


Influence of the Neglected Transient Term 


\lthough the burning-rate law (Eq. [19]) has been fre- 
quently confirmed experimentally (2, 3, 6 to 10), some of 
the implications of the quasi-steady theory are entirely in- 
correct. This fact is illustrated by Fig. 3b, which shows 
the variation of the fuel fraction f which follows from the 
solution for b with Equation [11]: Within the fuel particle 
itself, f = 1; at the phase boundary there is a discontinuity 
in f, which falls off gradually to zero as the radius increases. 

The quasi-steady theory leads to the conclusion that f de- 
pends on the ratio r/rs alone, i.e., that all the f profiles are 
similar. This means that since rg decreases with time, so 
does the width of the profile. Now since in evaluating f we 
take no account of whether the fuel is burned or unburned, 
it is clear that the area under the f profile is a measure of the 
amount of fuel in the original fuel particle; the vanishing of 
this area as time proceeds therefore contravenes the law of 
conservation of matter. 

This unsatisfactory feature results from the neglect of the 
transient term in the derivation of Equation [13] from Equa- 
tion [9]. This term becomes increasingly important (because 
of the multiplier p) as the pressure increases. In order to 
examine its influence, we now abandon the quasi-steady as- 
sumption. 


Transient Combustion of a Single Fuel Particle 


Outline 


In this section a fresh approach is made to the prediction 
of the burning rate of fuel particles introduced into an oxidant 
gas. The quasi-steady assumption is abandoned, and also 
its corollary that the rate of burning is identical with the 
rate of vaporization. This latter assumption is obviously 
useless when the phase boundary has ceased to exist. 

The new theory has two main parts. In the first, the 
variation of fuel fraction with position and time is established. 
This part makes use of a result from the theory of unsteady 
heat conduction. In the second part, the distinction is made 
between the fuel fraction and the mass fraction of unburned 
fuel; knowledge of the stoichiometry of the reaction then 
permits the first part of the theory to be used for the pre- 
diction of how the burning rate varies with time. The total 
time of burning is found incidentally. 

High speed computing machinery is needed to integrate 
the equations when full account is taken of the variation of 
the gas properties with temperature. This has not been 
done so far. However, the case of uniform properties is easy 
to evaluate analytically, and gives information which is ade- 
quate for present purposes. This is the case which is dis- 
cussed below. 

The main result is a new formula for burning time (Eq. 
[34]), which should be used in preference to Equation [22] 
at high pressures. 


Uniform-Property Solution for Fuel Fraction Distribution 


Differential equation and boundary conditions 


Since we are now not particularly concerned with the phase 
boundary, there is no advantage in using the variable b. 
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Instead the fuel fraction f will be considered. The corre- 
sponding conservation Equation [6], for uniform properties 
and spherical symmetry, becomes 


of 1a/.af 
— Dp. -—(r2—=) =0 

ot r? or ¢ [23] 
Here the convection term G- Vf has disappeared, because, with 
the initial fuel density assumed equal to the (uniform) gas 
density, G is zero throughout. This follows from consider- 


ation of mass conservation. — 
The boundary condition at large radius is 
e 
r= = fo [24] 


where it has been assumed, for generality, that the surround- 
ing gas already contains a certain amount of fuel. 

Another boundary condition is given by the fact that the 
total amount of fuel in the system is constant. This will 
be expressed by 


— fo) pdr = — fo) [25] 


where W,; is the mass of fuel injected at time zero sym- 
metrically about the point r = 0.4 

In addition, there are initial conditions (f as function of r 
at t = 0) which state how the fuel is distributed at the mo- 
ment of injection. The initial condition will be discussed 
later. 


A particular solution 


A well-known result of heat conduction theory concerns 
the distribution of temperature in an infinite solid after a 
sudden transfer of a finite amount of heat at a point in the 
solid. This can be used in the present problem. We obtain 
for the fuel fraction resulting from sudden injection of the 
quantity W, of fuelatr = 0 


9 
p(4arDt)*/2 


-fe= 
which is a solution to [23, 26 and 25], as may be verified. 


Fitting the initial conditions 


The starting condition of most practical interest is that in 
which the injected fuel is in the form of a sphere, of radius 
rs, in which f = 1, surrounded by as yet undiluted oxidant 
gas, in which f = fg. The corresponding solution for the sub- 
sequent fuel distribution can be obtained from Equation [26] 
by summing the effects of an infinite number of infinitesimal 
fuel sources, of appropriate magnitude and uniformly dis- 
tributed over the sphere of radius rs,. 

However it is easy to show that, for times of the order of 
rso?/D and greater, the fuel distribution is practically the same 
whether the initial source is concentrated at a point or is 
finite. We shall therefore use Equation [26] without further 
integration, and shall substitute 


where py, is the initial fuel sphere density, allowed here (and 
here only) to differ from that of the surrounding gas. Of 
course the use of Equation [26] is unrealistic at short times, 
because it then implies f values greater than unity, which are 
physically impossible. 

In this way we simultaneously sidestep the difficulty 
that the differential Equation [23] does not hold at the phase 

4The (1 — fe) term arises because the fuel, immediately on 


injection, comprises a volume in which f = 1; it is the part of this 
in excess of the average fuel fraction f¢ which diffuses out. 


a 
> 
he 
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boundary (if present), where there is a discontinuity in f. 
The simplifications which are being made are thus diametri- 
cally opposed to those of the quasi-steady theory; the present 
theory may be expected to become exact at just those condi- 
tions where the quasi-steady theory is inaccurate. 

Fig. 4 shows some successive distributions of f according 
to Equation [26]. The contrast with Fig. 3b is obvious. 


Variation of Flame Radius With Time: Burning Time 


In order to proceed further, it is necessary to make the 
assumption that the gas is everywhere in thermodynamic 
equilibrium. Then knowledge of the local value of f com- 
pletely specifies the local gas state, including its temperature, 
the degree of dissociation, etc. 


Stoichiometric considerations 


As an example, a fuel is considered which is capable of react- 
ing with the oxidant in only one way, and in which dissoci- 
ation is negligible. This means that at any point it is possible 
to have either unburned fuel or unburned oxidant, but not 
both, and that the mass fraction of unburned fuel m,; is given 
by 


my = (f — fst)/(1 — fst) [28] 
where f,, is the stoichiometric fuel fraction, and f> f... If 
f fu, mys = 0. 

Flame radius 


The flame is a surface separating the fuel-rich region from 
the oxidant-rich region, characterized by f = f.. Its posi- 
tion may be determined by substituting this value in Equa- 


Fig. 4 Successive fuel fraction distributions: Point-source 
theory 
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Fig. 5a Flame radius vs. time (Eq. [32]). 6 varies from 0 to 1.0 
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tion [26]. 
the injection point. 


The flame is evidently a sphere concentric with 
If its radius is r,,, we have 
in| W; (1 — fa) ] 
4Dt (far — 


29] 


This equation permits the flame radius to be evaluated :is 
a function of time. For ease of computation we introdu:e 
¢, the dimensionless flame radius, and 6, the dimensionl ss 
time, defined respectively by 


(Wy — fa)/plfet — fe)] /s 
“ 
i= 4r Dt 
Then Equation [29] becomes 


[: 


= [—(3/2r)6 In 6]'? 


This function is plotted in Fig. 5a. It shows that tie 
flame radius is zero at the moment of injection, increases tv a 
maximum at 8 = 0.368, and then diminishes to zero. This 
is illustrated by the sketches in Fig. 5b. Since the maximum 
value of ¢ is 0.419, the maximum radius attained by the flaie 
is 


stmax = 0.419(Wy(1 — fe)/p( fee — fe))'”* [33] 


This quantity therefore increases in proportion to the cube 
root of the injected mass and inversely as the cube root of the 
difference between the stoichiometric fuel fraction and the 
fuel fraction in the gas reservoir. If the latter gas consists 
of stoichiometric combustion products, the flame is infinitely 
far from the point of injection. 

Because the fuel is not really injected at a point, Equation 


[32] is in error for small 6. _ an 7 


One of the most important implications of Fig. 5a is that 
the flame radius diminishes to zero when 6 = 1. This time 
marks the instant at which all the injected fuel has finally 
been burned. The time for complete combustion ¢, is there- 


fore 
— fo)/p( fa — [34] 


4rD 


Time for complete combustion 


So the time of burning increases as the two-thirds power of 
the injected mass, and as the minus two-thirds power of the 
difference f,, — fe. In particular it appears that, if the gas 
reservoir consists of stoichiometric combustion product 
(fe = fe), the burning time is infinite; this occurs because, of 
course, there is no free oxidant in this case. 


Variation of Burning Rate With Time 


The region within the flame contains unburned fuel. Equa- 
tions [28 and 26] may be combined to enable the quantity 


NB. OROPLET DISAPPEARS 
BEFORE FLAME 


OROPLET 


Fig. 5b Appearance of burning droplet: Transient theory 
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of this fuel to be evaluated. Then the efficiency of com- 
bustion 7 can be related to time and other variables. 

If 7 is defined as unity minus the ratio of the fuel remaining 
unburned to the injected fuel quantity, we have 


Tst 
Pa f dr 
1—n=/0 [35] 
W,; 
Introducing Equations [28 and 26], Equation [35] becomes 


= 
st 
f fa — fa + (Wi 77/424 de 

Ss = 36 
— fu) 136] 
Afier some rearrangement, this equation reduces to 


n) (1 — fat) 
(1 — fe) 


3 '/2 
ot | 0) | AF 
6 \'/2 
6°/*((—In + (—In [37] 
where @ is once again the dimensionless time coordinate 
defined by Equation [31]. (1 — ) (1 — fi)/(1 — fe) is 
plotted vs. 0 in Fig. 6. 
xamination of Fig. 6 shows that the proportion of un- 
burned fuel diminishes to zero at @ = 1, which agrees with 
the conclusion about the burning time. At 6 = 0, (1 — 7) 
(1 — f..) = 1; this signifies that the inefficiency of combustion 
can be greater than 100 per cent. This conclusion is of course 
unrealistic, and arises from the use of the point-source formula 
Equation [26] without further integration; for this formula 
leads to f values on the axis which exceed unity at small 6, 
and so to negative oxidant fractions. It seems wisest to use 
only the right-hand part of Fig. 6, for practical calculations, 
and to imagine the left-hand branch to consist of a hori- 
zontal line at height (1 — f..)/(1 — fe). 
The burning rate is obtained by differentiating Equation 
[37]. We obtain 


(1 = fat) = Sut) an 27 1/2 3/2 
(1 —fg)d0 (=) [38] 


This function is also plotted in Fig. 6. It is seen to have a 


maximum at 8 = 0.0498, which in many cases will be in the — 


“unrealistic” period. 
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Fig. 6 Unburned fuel proportion and burning rate. 6 varies a ix: 
from 0 to 1.0 


NOVEMBER 1959 


Discussion 


Comparison of Present and Quasi-Steady Theories 


Comparison of flame behavior 


Figs. 3a and 5b illustrate qualitatively the differing pre- 
dictions of the two theories. According to the quasi-steady 
theory, the flame shrinks steadily in proportion to the droplet 
radius®; according to the transient theory, it increases at 
first and then decreases. Experimental data are available 
showing that, even at atmospheric pressure, droplet flames 
exhibit the latter behavior rather than the former (11). How- 
ever, if the droplet is introduced at a temperature lower 
than its wet-bulb temperature, the initial value of B is less 
than the final value; then even the quasi-steady theory will 
predict an initial increase of flame size while the droplet 
is heating up. These data cannot therefore be taken as 
confirming the importance of transient effects in the gas phase, 
although they certainly indicate their importance in the 


liquid phase. 
= 
Comparison of burning times a = 


In order to compare the burning time predicted by the 
quasi-steady theory (t))¢.s. with that of the | resent (point- 
source) solution (t,)p.s., Equations [22, 27 and 34] will be 
combined. We obtain 


9 99 (22) — fa)” 
“"\p In (1 + B) 


[39] 


Evidently, if this fraction is much greater than unity, the 
quasi-steady formula [22] should be used for burning time; 
if it is much smaller than unity, the transient formula [34] 
should be used. When the fraction is of the order of unity, 
both formulas will underestimate the burning time; in the 
absence of an analysis of the ‘mixed’ regime, it is probably 
wisest to take the sum of (fp)¢.s. and (tp) p.s. as the burning 
time. The corresponding formula for ¢, is therefore 


ty = + 2.82 (fe — fo) “/(1 — fe) | x 


p In (1 + B) 
— — fa)}*/* 
4rD “ted 


[40] 


Example 


To fix ideas, we consider a liquid oxygen droplet burning 
with kerosene at 40-atm pressure in a rocket motor with 
fe = 0.25 (fuel rich). For these propellants we have f.. = 
0.226. Other data required are: p,/p = 382, B = 134.2. 

Replacing f by 1 — f in Equation [39] (since an oxidant 
particle is in question), the ratio of burning times can now 
be determined; it is 0.88. So here the two predicted burning 
times are of the same order of magnitude for the oxidant. 
For the fuel, of course, the transient formula must be used, 
since kerosene has a critical point below 40 atm. 


Influence of pressure 
The density ratio of particle to gas appears in Equation [39] 
and causes the burning time ratio to decrease as the pressure 


5 Using the quasi-steady theory, and Equations [11 and 
27], it is easy to show that the flame radius r.; is given by 


(2 
4a Pp Tso? Pp 


In (1 + B) 
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rises, with fixed particle density; for p is proportional to pres- 
sure. The point-source formula is therefore the one to use 
at high pressures. 

There is another effect of pressure, however. As has been 
pointed out above, B tends to infinity as the critical pressure 
is approached, because the latent heat of vaporization 
vanishes. This effect causes the fraction in Equation [39], 
which reappears in Equation [40], to diminish still more. 
Once again it is the transient terms which dominate the 
solution. 

The transport property Dp is practically independent of 
pressure for gases. To show the effect of pressure on the 
burning time predicted by the point-source solution, we 
rearrange Equation [34] to give 

— fo)/(fa — 


= 
4m Dp [41] 


Since p is proportional to pressure, Equation [41] shows that 
the burning time of a particle of fixed mass increases in pro- 
portion to the cube root of the pressure. This throws light 
on the increasing relevance of the transient terms in the equa- 
tion from a different point of view, for the quasi-steady pre- 
diction of burning time is independent of pressure and is 
bound to be increasingly in error as the pressure rises. 
Influence of fg is 

In most engines, fuel particles are injected into gases con- 
taining an appreciable proportion of combustion products; 
in rockets for example, f., — fg may be of the order of 0.1. 
Now it is easy to prove that such a dilution of the oxidant gas 
has very little effect on the value of the transfer number B, 
and so still less on In(1 + B). It follows that dilution of 
oxidant gas by combustion products decreases the ratio in 
Equation [39] and again makes the burning time nearer to 
that predicted by the present theory. 


Influence of particle size 


In one respect the quasi-steady and transient formulas 
agree: The burning time is proportional to the square of 
the diameter of the injected particle, ie., to (W,)’*. So 
fine subdivision of the injected fuel is important for fast 
burning, whichever formula is appropriate. 


Influence of Approximations Made During the Analysis 


It should be noted that the analysis leading to the new 
burning time formula [34] involves a number of departures 
from reality. These include: 

1 In the differential equations of conservation, the main 
relevant assumption is that all chemical substances possess 
identical diffusion coefficients. 

2 Thermodynamic equilibrium has been supposed to 
prevail throughout. 

3 This equilibrium assumption has been supposed to 
imply that unburned fuel and oxidant cannot have finite 
concentrations at the same point; dissociation and the possi- 
bility of more than one set of reaction products have been 


neglected. 
4 D and p have been assumed uniform throughout the 
gas. 


5 The injection of a fuel particle of finite size has been 
simulated by an instantaneous point source of fuel. 

6 The mutual influence of neighboring droplets has been 
neglected. 

We will discuss the influence of each of these assumptions in 
turn: 

1 This assumption is necessary and usual. It involves 
only a small departure from reality unless very heavy or very 
light molecules are present. 

2 This assumption probably accords well with reality in 
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rocket motors owing to the rapid chemical reaction rates per- 
taining to high pressures and temperatures. 

3 This assumption represents too crude a picture. In 
rocket motors, for example, the addition of further fuel to a 
stoichiometric fuel-oxidant mixture would cause further 
reaction. However the present assumption has the advan- 
tage of clarity and approximate correctness. 

4 and 5 These assumptions are certainly incorrect, and 
cause the theoretical predictions for the early stage of burning 
to be in error; they have much less effect at larger times. 
The assumptions can be removed without difficulty if suf!:- 
cient computing effort is available to permit numerical solii- 
tion of the nonsimplified equations. 

6 This neglect, together with assumption 3, is responsib'e 
for the implication that no burning takes place if fg is already 
equal to f.:. In reality of course, fuel particles can bwin 
under these conditions because oxidant particles are prese::t 
within a finite distance. This matter has been analyzed |v 


the author and will be reported elsewhere. oat h 


Conclusions 


The quasi-steady formula for the burning time of a fuel 
(or oxidant) particle underestimates the burning time when 
the pressure is large, the transfer number becomes high and 


_ the gas condition approaches the stoichiometric. 


The recommended formula for the burning time of a fuel 
partic le, irrespective of conditions, is Equation [40]. This 
is valid far above and far below the initial point; its validity 
when both terms are of the same order requires further ex- 
amination. 

The burning time increases the nearer to stoichiometric is 
the fuel-oxidant ratio of the atmosphere surrounding the 
particle. 

It will be necessary to re-examine, in the light of the pres- 
ent theory, the conclusions drawn from studies of the role of 
particle burning in the performance of liquid fuel rocket 
motors (27 to 29). Although the present theory agrees with 
the earlier ones in showing that fine fuel subdivision is favor- 
able, it shows that the important propellant properties may 
be considerably different from those which have previously 
attracted attention. 
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Nomenclature 


dimensionless generalized mass transfer variable 

(= be — bs), transfer number 

specific heat at constant pressure 

diffusion coefficient 

fuel fraction (mass of fuel, whether burned or not, in 
unit mass of mixture) 

enthalpy 

thermal conductivity 


mass fraction of chemical element @ in mixture : 
mass of element a in unit mass of compound j 


pressure 
general conserved property 

heat transfer from gas per unit mass changing phase 
radius 

time 
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T = temperature 

m” = mass flow rate per unit area through phase boundary 
G = mass flow rate vector 

W = injected mass 

= efficiency of combustion 

= — fe)/p(fet — fa) 

p = density ; 

¢ = — Fo)/o( Fat — 

Subscripts 


= chemical element 

= complete burning 

= fuel (unburned) 

= fuel reservoir 
= gas reservoir 

chemical compound 

= mixture, mean = 
= oxidant 
= injected particle 

= point-source theory 

= quasi-steady theory 

= gas in equilibrium with surface 
= initial surface 
stoichiometric 

= transferred substance 
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Vaporization of Propellants 
Rocket Engines 


RICHARD J. PRIEM' and 
MARCUS F. HEIDMANN! 
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Cleveland, Ohio 


The theory and equations required to calculate the vaporization rates in a rocket engine, based on 
a model considering propellant vaporization as the rate controlling combustion process, are prc- 


sented. Included in the model are changes in combustion gas velocity, droplet velocity, temper- ‘ 
= ature and mass. The vaporization rates of heptane, ammonia, hydrazine, liquid fluorine and 
ah A liquid oxygen sprays were computed for various log probability distributions and mass median f 
_ drop sizes of the spray. The rates were also calculated for an engine operating with various injec- t 
a tion velocities, final gas velocities, initial propellant temperatures and chamber pressures. The . 
4 =. calculations of the per cent of propelJant vaporized are correlated with an effective chamber length. t 
The correlated results show that the percent of propellant vaporized can be increased either by de- t 
creasing the mass median drop size and the initial drop velocity or by increasing the chamber length, 8 
chamber pressure, final gas velocity and initial fuel temperature. A methodof analysis that relatesex- 1 
perimental rocket engine performance to the quantity of liquid propellant vaporized is also presented. U 
With this method experimental engine performance results are compared with the calculations. The 7 
comparison suggests that incomplete vaporization is responsible for combustor inefficiency. Ex- a 
perimental results agree with the calculations for sprays having a geometric standard deviation 
of 2.3 and a mass median drop radius of 7 70 to 240 microns, depending on the type of injector. fi 
he 
HE LARGE number of different phenomena that can good agreement with the experimental results of single-drop 
have a fundamental role in the combustion within a liquid investigations over a wide range of conditions that approach 
propellant rocket engine makes it difficult to determine and those encountered in rocket engines (9, 15 and 16). Therefore, 
study the controlling process. Some of these phenomena, as this technique (14) was used to calculate the vaporization T 
given in (1 and 2),? are atomization, heating, vaporization, rates of sprays by using a model which assumes that the N 
liquid and gaseous diffusion, combustion of droplets, and vaporization of the propellants is the rate controlling step in m 
liquid or gas phase reaction. The consideration of all of these rocket engine combustion. The analysis also assumes one- aX 
phenomena results in a large number of similarity parameters dimensional steady-state flow with no interaction between as 
(1 to 4) that should be used in designing and scaling rocket drops and no shattering of drops. Many combustor param- of 
engines. eters were varied to show how they affected the combustion th 
The inherent complexity of the problem, when all the chamber length required to vaporize the propellant. The 
phenomena are considered, make it desirable to simplify the calculations were made for five propellants: Heptane, 
problem. Simplification is obtained by neglecting those ammonia, hydrazine, liquid oxygen and liquid fluorine. From 


processes that require only a small time or distance in the 
chamber, and by considering only the phenomena that require 
the greatest time or distance. Recent studies (5 to 9) show 
that the atomization, heating, gaseous diffusion and gas 
phase reaction require only a small fraction of the rocket 
chamber length under normal operation. Liquid diffusion 
and reaction are only encountered with premixing rocket injec- 
tors using certain hypergolic propellants, as nitric acid, with 
some fuels. Therefore, vaporization and combustion of drop- 
lets remain to be evaluated. The results reported in (10) 
indicate that vaporization rate equations evaluate both the 
vaporization and the combustion of droplets. The results of 
an analytical study of the importance of propellant vaporiza- 
tion rates on rocket engine combustor efficiency are reported 
herein. 

Of the various techniques (11 to 14) for determining vapor- 
ization rates, the one given in (14) seemed most applicable to 
the conditions encountered in rocket engines. This technique 
permits consideration of the heating period of the drop, uni- 
directional diffusion and the effect of mass transfer on heat 
transfer. Furthermore, the results of the technique are in 
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Nov. 17-21, 1958. 
1 Aeronautical Research Scientist. 


2 Numbers in parentheses indicate References at end of paper. 


836 


these results a general correlation was developed to ade- 
quately cover the range of parameters studied. 

By assuming that droplet evaporation is the rate limiting 
step in rocket combustion, the performance of a propellant 
combination can be either correlated or predicted from cal- q 
culated vaporization rates. This was done for the heptane- 
liquid oxygen system. The experimental performance of one 
injector was correlated over a wide range of operating condi- 
tions, and the performance of six different injectors was pre- 
dicted and compared with experimental values. 


Theory 


The calculation technique was based on applying known 
steady-state mass, momentum and heat transfer equations to 
droplets vaporizing in a rocket engine. The first order differ- 
ential equations were set up using the model shown in Fig. 1 
for a single drop. 

A liquid droplet is shown at position x and at a small incre- 
ment later in time A@ and distance Az. During this incre- 
ment the drop velocity v changes by a small amount Av. In 
this same increment, the gas velocity u changes by Au. The 
drop velocity would start at the injection velocity and either 
slow down or speed up, as determined by the drag on the 
drop. The gas velocity is zero at the injector and increases 
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downstream as the propellants vaporize and burn. While - 


drop is moving through the increment, heat is transferred to_ 
the liquid surface at a constant rate g,, and mass is being 


transferred from the surface at a constant rate w; oxidant is _ 


added to burn the vaporizing fuel (or vice versa) instantane- 
ously at stoichiometric conditions. During the increment, the 
mass and the radius of the drop change from m, to mz and 7, 
to m2, respectively, while the temperature 7; changes by 


AT,. 
The equations used to determine the mass transfer rate are: 
err. 
w = A,Kpa (1] 
jum = —— = 2 *Re'/* 2 
a=-—In [3] 


The Ranz and Marshall (17) correlation (Eq. [2]) was used 


for the Nusselt number Nuy, since a heat transfer and mass : 
ve 


transfer correlation had been determined therein under con- 
ditions involving simultaneous heat and mass transfer. Also, 
the nature of their investigation was similar to that studied in 
this paper, although the maximum temperatures and _pres- 
sures were not as high as those used in this investigation. 
The mass transfer equation included the a term to account for 


unidirectional diffusion in the theory presented in (14 and 18). ea 


The @ factor is always greater than 1 and ranged between 1 
and 3 in these calculations. 

The rate of heat transfer to the liquid surface is obtained | 
from the following equations 


qv = AAT, — T)Z 


Nun = h2r/k [5] 
Z = 2/(e? — 1) [6] 
= Wey, s/hA, [7] 


The Ranz and Marshall (17) correlation for the heat transfer 
Nusselt number Nu, was used for the same reasons given for 
mass transfer. The heat transfer included the Z term to 
account for the sensible heat taken up by the diffusing vapor 
as required in the theory presented in (14 and 19). The value 
of Z is always less than 1 and ranged between 1 and 0.05 in 
these calculations. 
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Schematic model of fuel or oxidant drops vaporizing in 
rocket engine 


Fig. 1 


where 


gt = — 


The drop velocity was calculated from equations for aero- 
dynamic drag 
mdv 


= > {10} 


The drag coefficient obtained by Ingebo (20) was used as it 
was determined for evaporating sprays over a Reynolds num- 
ber range of 6 to 400 Oo 


Cp = 27 Re-?.84 [11 | 
No correction was added for the effect of mass or heat transfer 
on the drag coefficient. The equation was also used for 
Reynolds number as large as 2000. 

The gas velocity equation is derived from the assumption 
that stoichiometric reaction follows instantaneously after 
evaporation; thus the gas phase consists only of reaction 
products. Therefore, gas velocity, developed by mass addi- 
tion at constant temperature and pressure, is directly propor- 
tional to the fraction of propellant evaporated and is given by 
Equation [12], which was derived in (21) for a single drop 


The change in drop temperature is determined by and in (22) for sprays ps 
IT 
a q 
[8] =1- Mi ms [12] 
Table 1 Range of conditions used for calculations 
Mass 
4 median Initial 
—— Initial drop temperature, R—————— drop Geometric drop Final gas Chamber 
Hep- Ammo-  Hydra- Fluo- radius, standard velocity, velocity, pressure, 
tane nia zine Oxygen rine in. deviation in./sec in./sec lb/in.? 
400 300 0.003 2.3 1200 9600 300 
5007 40024 5007 140? 0.0037 2. 30 12007 96007 3007 
0.001 i 1200 9600 300 
0.003 600 300 
| | 1200 2400 300 
| 9600 150 
9600 600 
19,200 30000 
| | f 300 
0.009 2.3 1200 9600 30000 
700 500 700 220 220 0.003 2.3 i: 1200 3 9600 300 
4 
* Conditions used for calculations presented in Figs. 3 to 6. oo 
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For calculation purposes the propellant was assumed to be 
in five groups of drop sizes, which were arbitrarily selected to 
correspond to a log-probability distribution as described in 


dr +r 2 In oy 


The five sizes were chosen from the log probability distribu- 
tion as those equal to 10, 30, 50, 70 and 90 per cent of the mass 
in drops smaller than the drop radius. The number of drops 
in each group was so chosen that each group had 20 per cent 
of the total mass. Plots of the distribution obtained by this 
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Fig. 3 Typical droplet temperature histories for various propel- 
lants. (Conditions used in calculations given by (*) in Table 1] 


method and from Equation [13] are shown in Fig. 2 for a 
geometric standard deviation o, of 2.3. The geometric stand- 
ard deviation is a spray characteristic where 68.23 per cent of 
the mass is in the drop size range between M,, w/o, and 
o,M M- 

The vaporization equations were solved numerically by an 
iterative procedure to obtain the per cent of propellants 
vaporized in various chamber lengths. Table 1 shows the 
spray conditions as well as the design and operating para:n- 
eters used in the calculations. Physical properties used in 
the calculations for the various propellants are given in the 
appendix of (24). 

The relation between combustor efficiency n- (charact:r- 
istic exhaust velocity efficiency) and per cent of propella: ts 
vaporized was obtained by considering the exhaust to bc a 
mixture of oxidant drops, fuel drops and gaseous combusti »n 
products. The gaseous combustion products were assumed to 
be in thermodynamic equilibrium, and the volume and kine ic 
energy of the drops were assumed to be negligible. Efficier »y 
is then related to vaporized propellant by the following eqi:a- 
tion, derived in (25) 


(1) (ctn* ovr Wo wy 


Results and Discussion 


Analytical Results 


The results of the vaporization calculations are dropiet 
histories, examples of which are shown in Figs. 3 to 6 for the 
various propellants. The conditions for which these calcula- 
tions were performed are indicated by asterisks in Table 1. 
The calculated temperatures of the mass median drops are 
shown in Fig. 3. With all propellants, the temperature of the 
droplet rises to within 1 degree of a constant temperature 
(wet bulb) in less than one tenth of the distance required to 
vaporize the drop. 

Typical histories of the per cent of propellant vaporized for 
various propellants are shown in Fig. 4 for the mass median 


80 


oO 


PERCENT VAPORIZED 
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Fig. 4 Typical droplet mass histories as shown by per cent 
vaporized for various propellants. [Conditions used in calcula- 
tions given by (*) in Table 1] 
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Fig.5 Droplet velocity histories for various propellants. [Condi- : 
tions used in calculations given by (*) in Table 1] 


drops. The curves are similar for all propellants. The oxi- 
dants—oxygen and fluorine—vaporize in shorter lengths than 
the fuels—heptane, ammonia and hydrazine; the length in- 
creases in the order listed. 

Drop velocity histories of the mass median drops are given 
in Fig. 5. The drop velocity first decreases because of the 
drag produced by the low gas velocity. After reaching a 
minimum, the drop velocity is increased because of increasing 
gas velocity. 

Per cent of total mass vaporized is shown in Fig. 6. Since 
the per cent of mass vaporized (for all drops) is proportional 
to the gas velocity (Eq. [12]), the gas velocity is also shown 
in this figure. The per cent vaporized in a given length was 
greatest with oxygen and fluorine and decreased with heptane, 
ammonia and hydrazine, in that order. Comparison of the 
per cent of the mass median drops vaporized (Fig. 4) with the 
total mass vaporized (Fig. 6) shows that with heptane the 
mass median drops are completely vaporized in 9 in., but that 
only 88 per cent of the total mass is vaporized in that length. 
This difference is due to the large drops that are only partially 
vaporized. It is the large drops in the spray, therefore, that 
reduce combustor efficiency. Curves showing the effect of 
various operating parameters on the per cent vaporized for 
single heptane drops may be found in (21). The effect of 
changing the propellant spray characteristics is given in (22). 

The relation between per cent fuel and oxidant vaporized 
and combustor efficiency, as determined from Equation [14], 
is shown in Fig. 7 for heptane and oxygen. For equal frac- 
tions of fuel and oxidant vaporized, the combustor efficiency 
is equal to the per cent vaporized. Each value of efficiency, 
however, may be obtained with various relations between the 
per cent oxidant and fuel vaporized. Additional curves relat- 
ing performance to propellant vaporized for various mixture 
ratios can be found in (25) for hydrogen-fluorine, hydrogen- 
oxygen, ammonia-fluorine and JP-4-oxygen propellant com- 


binations. 
3 
Correlation of Analytical Results 


Vaporization histories similar to those in Figs. 3 to 6 were 
obtained for heptane, ammonia, hydrazine, oxygen and 
fluorine as listed in Table 1. A range in chamber pressure of 
150 to 600 psi, median drop radius of 0.001 to 0.009 in., drop- 
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Fig. 6 Typical average mass histories, as shown by per cent 
vaporized. [Conditions used in calculation given by (*) in Table 1] 


let velocity of 600 to 2400 in. per sec, and gas velocity of 2400 
to 19,200 in. per sec was covered. An attempt was made to 
fit the results for each propellant over these ranges of condi- 
tions to a single curve, usable for engine design. The correla- 
tion of results was accomplished by modifying chamber 
length with a correction factor to obtain an effective chamber 
length 


Effective length = 


1.9 1075/(1 — [15] 


A plot of the per cent mass unvaporized as a function of 
effective length for heptane sprays with a distribution of 2.3 
is shown in Fig. 8. The log-log plot is used in Fig. 8 and in 
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Fig. 7 Percentages of vaporized fuel and oxidant required to 
maintain a given C* efficiency. Propellants, heptane and oxygen; 
oxidant-fuel weight ratio, 2.2 
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Correlated results for all combustor conditions with hep- 
tane sprays 
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Fig. 9 Correlated results for all combustor conditions with vari- 
ous propellant sprays which have a particular distribution. 
(Standard deviation of 2.3) 


100-— LENGTH,| MIXTURE | FUEL FLOW |TOTAL FLOW 
IN. 
— BO lg 2 1.37 to 3.41 | 0.198 to 0.383 0.90 
L lo 4 1.19 to 3.84| 0.193 to 0.420 0.90 
so. |° 6 1.41 to 3.58 | 0.196 to 0.375 0.90 
A 8 1.40 to 3.60 | 0.196 to 0.374 0.90 
2) eo 4 1.00 to 4.05 2.60 0.565 to 1.19 
© og 4 2.5 0.148 to 0.369 0.90 
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= 4 
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Fig. 10 Combustor efficiency for various chamber lengths and 
flow conditions 
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subsequent figures to obtain better resolution in the 90 to 
100 per cent vaporized region. The combustor length required 
to vaporize a given percentage propellant increases with 
increasing drop size and injection velocity, and with decreas- 
ing gas velocity and chamber pressure. The spread in effective 
length required to vaporize a given per cent of propellant for 
various conditions was approximately +20 per cent about the 
log mean curve. 

The spray distribution of 2.3 used to obtain Fig. 8 and su)- 
sequent figures corresponds to that observed in (6) for an 


impinging-jet injector. 


The effect of using a different distri- 


bution is illustrated in Fig. 8. Additional calculations can |e 
found in (24). 

The correlation of per cent of mass unvaporized with effe:- 
tive length for various propellants is shown in Fig. 9. Only 
the log mean curve is shown. The spread in effective length 
for various conditions was always less than +20 per cent for 


all propellants. 


The curves for fluorine and oxygen are almost 


identical with the effective lengths required to vaporize a 
given percentage of propellant increasing in order with hep- 


tane, ammonia and hydrazine. 


Correlation of Combustor Data 


The correction factor to obtain an effective lengt 


(Ey 


[15]) may also be used in the correlation of experimental co1n- 


bustor data. 


Experimental heptane-oxygen combustor ef!i- 


ciencies for a given injector, in which propellant flows, mixture 
ratios, chamber pressure and chamber length were varied, are 
shown in Fig. 10. The injector used two impinging jets of 
heptane to produce a sheet and was similar to that used in 
(6). 
small injector orifices and, therefore, could be assumed to be 
completely vaporized in the shortest chamber investigated. 
The plot of combustor efficiency against chamber length (Fig. 


10) shows that a wide 
with the various conditions. 


The liquid oxygen was finely atomized by using very 


variation in efficiency was obtained 
Converting combustor efficiency 


to per cent fuel vaporized by use of Equation [14] resulted in 


the correlation of the data as shown in Fig. 11. 


Modifying 


the chamber length by the correction factor to obtain an 
effective length (Eq. [15]) produces the curve of per cent of 
fuel vaporized against effective length shown in Fig. 12. The 
absolute value of the results agrees with the correlated ana- 
lytical results shown by the shaded area in Fig. 12. The 
spread in the data was less than that of the correlated results. 


The mean drop size M,, » for the 


various data points was 


calculated from the following equation obtained from (6) 


4 = 1.12VD,V; +0108 D,|V,— Vi] [16] 


= 
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Fig. 11 Propellant vaporized for various chamber lengths and 


flow conditions 
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Table 2 Drop sizes produced by various injector types 


Mass mean Mass mean 
Injector type radius, pu radius, ref. 6 
jets 240 OF 
sheets 95 92 
triplet 70 


In this equation the diameter of the jet D; and the jet 
velocity V, are directly calculated from the experimental 
conditions. However, the gas velocity V, cannot be deter- 
mined in this manner. This velocity increases from zero at 
the injector face to the final velocity at the point where all 
the heptane is evaporated. An average value of 1200 in. per 
sec was assumed for this calculation. 


down the chamber, or it could be the velocity perturbations 
produced by small pressure fluctuations. 

The comparison of correlated experimental and analytical — 
results on the basis of propellant vaporization exhibits excel-_ 
lent agreement for an injector where cold-flow studies have 
established drop sizes. All variations in experimental per- 
formance with different operating conditions are predicted 
analytically by changes in the vaporization rate as a result of 
changes in operating conditions. The variations are a func- 
tion of the difference in drop size, injection velocity, chamber 
pressure and final gas velocity produced with different operat- 
ing conditions. 


Predicted Combustor Efficiency 


The analytical results are also sufficient to predict com- 
bustor efficiency for various injection techniques. Experi- 
mental heptane-oxygen combustor efficiencies of six injectors 
with chamber lengths of 2, 4 and 8 in. at an oxidant to fuel 
weight ratio of 2.5 have been presented in (26). 

The combustor efficiencies for these injectors in a given 
chamber length may be predicted from a knowledge of the 
drop sizes, the injection velocity and the initial temperature 
of both the oxygen and heptane, as well as environmental 
conditions of chamber pressure and final gas velocity. These 
parameters specify an effective chamber length as expressed 
by Equation [15]. The injectors used in (26) differ primarily 
in the drop size produced by the various atomization methods. 
The effective length can be determined for each injector by 
assuming that the same drop size was produced with oxygen 
and heptane for each method of atomization and by assuming 
a drop size for each method, as shown in Table 2. The only 


This velocity corre- 
sponds to the velocity attained about one quarter of the ond 


ANALYTICAL REGION 
FROM (FIG 8) 


PERCENT UNEVAPORATED 


LENGTH | MIXTURE |FUEL FLOW/TOTAL FLOW 


2 VARIABLE} VARIABLE | CONSTANT 


20 


CONSTANT VARIABLE 
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\ 2 3 4 5 678910. 
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Fig. 12 Comparison of experimental and analytical results for 
various chamber lengths and flow conditions 


check on drop size is with the sheets or impinging-jet atomizer, 
where the size agrees with that obtained by Equation [16]. 
The velocity of the gas into which the propellant is injected 
again was assumed to be 1200 in. per sec. 

A summary of the analysis to predict the combustion effi- 
ciency of the various tests made in (26) is given in Table 3. 
The combustor variables are given along with the effective 
lengths calculated for each propellant. From these effective 
lengths, the per cent of each propellant vaporized was taken 
from Fig. 9. 

In all cases where the atomization technique was the same 
for both propellants, somewhat more oxygen was vaporized 
than heptane, although the difference was not great. From 
these percentages vaporized, the combustor efficiency was 
determined by using Fig. 7. The predicted results are com- 
pared with the experimental results in the last two columns 
of Table 3. The calculated and experimental values agree 
within +4 efficiency units, which is within the accuracy of the 
analytical correlations. 

The comparison of experimental and calculated efficiencies 
in Table 3 also shows that there is a wide variation between 
per cent of oxidant or fuel vaporized and combustor effi- 
ciency. To obtain high efficiency, both propellants should be 


Table 3 Comparison of experimental and calculated performances — 
Calculated results— 
Chamber Oxidant Fuel 
Injector type length, Effective Per cent Effective Per cent Experimental? 
Oxidant Fuel in. length vaporized length vaporized Ne 
jets jets (wide) 8 1.08 1.43 38 41 
jets jets (close) 8 1.15 1.22 33 47 
sheets jets “3 5.1 1.22 33 51 
jet sheets 8 1.27 5.65 68 
sheets sheets “ee 2.10 2.15 68 
4 3.22 3.33 76 
a. 5.5 5.05 82 
triplet triplet 7 2 3.22 3.33 81 
5.1 5.25 
8 7.4 8.00 88 
Ref. 26. 
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vaporized. However, it is possible to obtain an efficiency that 
is higher than the per cent fuel vaporized if the oxidant is 
almost completely vaporized, or higher than the per cent 
oxidant vaporized if the fuel is almost completely vaporized. 
The results of the comparison also indicate that better 
atomization of the propellants is responsible for the higher 
efficiency obtained with the triplet injector. Additional 
studies on drop sizes produced by various injector techniques 
are still required to either verify or disprove this point. Such 
data would also be valuable as an aid in engine designing. 


Summary of Results 


Calculations were made to determine the vaporization rates 
of sprays in a rocket engine for five different propellants. 
The calculated results show that the combustor length re- 
quired to vaporize a given high percentage of propellants is 
the shortest with oxygen and increases in order for fluorine, 
heptane, ammonia and hydrazine. The combustor length 
required to vaporize a given percentage of propellants de- 
creases with decreasing drop size and injection velocity, and 
with increasing gas velocity and chamber pressure. 

Comparisons of experimental and calculated results showed 
excellent agreement for injector types where the drop size 
could be calculated. Variations in engine performance with 
various operating conditions were explained by determining 
the relation between performance and the per cent of propel- 
lants vaporized. The results indicate that combustor effi- 
ciency can be predicted from a model based on propellant 
vaporization as the rate controlling combustion proses. 


_ for Pure Liquid Drops,” 


Nomenclature 

A, cross-sectional area of drop, in.? 
A, = surface area of drop, in.? 

a = constant for mass distribution, dimensionless _ - 

Cp = coefficient of drag for sphere, dimensionless : 
Cp = specific heat at constant pressure, Btu/(lb) (F) 

= characteristic exhaust velocity, f 

D = diffusion coefficient, in.2/sec 

D; = jetdiameter, in. ates" 

e = Naperian base > oa 

F = fuel vaporized, per cent = 

h = heat transfer coefficient, Btu/(in.2)(sec)(F) = 

K = coefficient of mass transfer, 

k = thermal conductivity, Btu/(in.) (sec) (F) —T 

L = chamber length, in. 

M = molecular weight of propellant, lb/mole 4 hae 

= Mass median drop radius, in. 
m = mass of drop, lb ‘ 
Nuz = Nusselt number for heat transfer, dimensionless & ; 
Num = Nusselt number for mass transfer, dimensionless v 3 
n = number of drops of a particular size ' 4 
= oxidant vaporized, per cent 
= total pressure in chamber, psi ; 
p = vapor pressure of liquid, psi 

Pr = Prandtl number, dimensionless 

= sensible heating rate of drop, Btu/sec 
Qe = heat transfer rate to surface of drop, Btu/sec —s_— 

R = per cent of mass in drop smaller than r 

R = universal gas constant, lb/(mole) (F) _ 

Re = Reynolds number, dimensionless 

r = radius of drop, in. 

Sc = Schmidt number, dimensionless 

T = temperature, R 

= critical temperature of propellant, R 

T; = reduced initial propellant temperature, Tr, o/Te 

§ = mean temperature in boundary, R 

U = velocity difference between gas and drop, in./sec 

u = velocity of gas in chamber, in./sec ig 

Utin = final velocity, in./sec 

Vy = gas velocity, in./sec 

V; = jet velocity, in./sec 

v = droplet velocity, in./sec 

w = vaporization rate, lb/sec 


Wo = oxidant weight flow, lb/sec 

wy = fuel weight flow, lb/sec 

x = axial position in chamber, in. 

Z = correction factor for heat transfer, dimensionless 

a = correction factor for mass transfer, dimensionless 

n = characteristic velocity efficiency, per cent of theoretical 
6 = time, sec 

r = latent heat of vaporization, Btu/lb 

p = density, lb/in.* 

9 = geometric standard deviation, dimensionless 


Subscripts 

j = vapor 

a = ithsize drops 


O/F = oxidant-fuel weight ratio, ti 
0/F = ratio of vaporized oxidant to vaporized fuel weig! t 


flow, Otio/Friy 
th = theoretical 
= beginning of time 4 
- 
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Maximum Coasting Apogee Velocity ross: 


Lockheed Missile Systems Division 


Palo Alto, Calif. 


: te the present paper, thrust direction (angle of attack) is optimized during the second of two ad- 


At second stage burnout, a coasting ellipse is to carry a third and 
We seek the thrust direction program which will 
maximize velocity at the apogee of this coasting arc. An expression for this program is obtained, 
and assumes a form noted in several previous papers. Supposing second stage thrust to be of con- 
stant magnitude, the equations of motion there are integrable in closed form. During this phase 
only, the Earth is assumed flat; aerodynamic effects are neglected beyond first stage burnout. 
Under these assumptions, determination of the optimum trajectory is reduced to a one-parameter 
iteration, readily effected by machine or hand calculation. Forms are derived for the second varia- 
tion and also for the total variation. These expressions permit the investigation of conditions which 


jacent satellite ascent stages. 
final stage to some prescribed orbit height. 


Perrone problems of wide engineering interest very 
A often involve not the design of totally radical systems 
to perform a new task, but rather the adaptation of those al- 
ready in existence. We cite, as an example, the present use of 
large military rockets as satellite launching vehicles. The 
trajectory analyst often encounters fundamental difficulties 
in applying available booster stages to a job for which they 
were not designed. He may wish to program thrust or lift 
according to some preconceived optimizing relation, but prac- 
tical limitations can narrowly restrict his choice. A typical 
dilemma arises, which furnishes the subject of this paper. 

A certain long-range rocket is to carry some fixed final 
stage mass to a specified perigee injection altitude. Because of 
aerodynamic instabilities associated with any but the most 
insignificant angles of attack, the booster stage must ascend 
on a zero-lift trajectory, at least until detrimental effects be- 
come negligible. After a certain condition is reached, e.g., at- 


tainment of some given height, air density, time or dynamic 

pressure, the angle of attack (thrust angle) may be varied at 7 di = < — (T/m) cos a = 0 [1] 

will. Thrust magnitude is to be held constant, however, and +9 —(T/m)sina = 0 [2] 

may not be altered. Following booster engine burnout, the 

final stage is to ascend to the given injection height by means where x and y are the rectilinear position coordinates. Second 

of a gravitational coasting ellipse. The final stage engine stage motion is assumed to occur over a flat, airless Earth, 

provides a velocity increment at coast apogee, sufficient to and under the influence of a constant gravitational accelera- 

establish an orbit. Clearly, the smaller the velocity increment tion g.* T is a constant; m(t) = mo — riot, where mo and tito 

required, the greater the useful final stage payload which may are given constants, and ¢ denotes time. a@ is the angle of 

be realized.2- The problem, therefore, becomes that of choos- attack (thrust angle), referred to the z-axis. 

ing a thrust angle program during the latter portion of the eee 

boost (this latter portion we shall call the “second stage,” as B. 

ight Trajectory of a Satellite, OPULSION, Vo no 

indeed may actually be the case), so that maximum velocity May 1008, po. 502-048. ‘The author of the note discunee Sat 

be attained at the coasting apogee. : Earth assumptions. He concludes that these approximations are 
The notion of partial optimization of a composite ascent valid, providing that distance covered during this phase is small 

path was mentioned by Perkins (3)* and also by Fried (4). compared t to = a radius be that . — ty is much less 

. ‘ ’ than min, which we assume to be the case here. 

Their results, along with those of a number of other autho ay As a matter of fact, the following bounds on burnout condition 

(2-5) arose in the consideration of optimum thrust direction errors can be shown to hold 

programs for a flat, airless Earth. The model thus established f. . 

indicates a thrust angle program whose expression is quite sim- —1/2[w(t, — to)] < — S — 1/2[w(t, — to)] 

Presented at the ARS 13th Annual Meeting, New York, N. Y., — S (ye — Yo) S — to)] 

Nov. 17-21, 1958. 1/3[w%(ty — < — ys) < 1/3[wXt — bo) 5] 


1 Senior Scientist. Member ARS. 

2 We assume, as a convenience, that fuel and hardware in the 
third stage are interchangeable, supposedly possessing average 
densities which are equal. 

3 Numbers in parentheses indicate References at end of paper. 
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insure mathematically optimum paths. 
nounced Vanguard data as parameters. 
- compared with two non-optimal programs. 


Finally, a sample calculation is presented, using an- 
Results for the optimum sajeutany are obtained and 


ple. Furthermore, for constant thrust, the equations of motion 
during optimized powered flight are integrable in closed form 
and are easily adaptable to rapid calculation procedures. 
The main drawback, of course, is due to the neglect of drag 
forces and Earth’s curvature; the inclusion of these (6, 7) 
renders the problem intractable. Here, we propose to use an 
exact analysis for the first boost stage, and also for the coast- 
ing phase. The second stage is to make use of the preceding 
simplifying assumptions. This composite method combines 
ease of computation with a close adherence to physical reality. 


Equations of Motion 


Referred to a fixed rectilinear coordinate system located at 
the point of second stage ignition (the “staging point’’), the 
equations of motion for an ascent rocket, treated as a point 
mass, are (see Fig. 1) 


where the quantities with a bar superscript are corrected values 
obtained by perturbing a a seco {1] and [2] with terms of order 
1/R, in the gravitational force expansion. Here, Rw* = 4g, 


and we assume jj > 0 in (to, ts). 
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COASTING 
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Nomenclature and conventions 


Zero-lift trajectory characteristics for the first boost stage 
are assumed known, and these may include drag and 
Earth curvature effects, if so desired. For the present problem 
we suppose that staging occurs after a given time of first 
stage flight. The velocity components, height and mass at 
staging then are known functions of some first stage parame- 
ter, such as yo, the staging height. This one-parameter set of 
burnout values provides initial conditions for the second stage 
motion, which is to be optimized. At second stage burnout, 
an inverse-square gravitational field is assumed, and the burn- 
out conditions here completely determine an ellipse whose 
apogee is interpreted to be the orbital injection height. 


The Variational Problem 


Maximizing tangential velocity at a coast apogee of given 
radial distance from the Earth’s center is equivalent to 
maximizing angular momentum at that distance. And since 
the coasting ellipse preserves angular momentum, this is, in 
turn, equivalent to maximizing angular momentum at second 
stage burnout. Because of the assumption that (1/R.)(x»—20) 
is negligible, this latter quantity may be written as 


My = (Re + yo) ae [3] 


where R, is the length of a radius vector from Earth’s center 
to the x,y origin. y and z are burnout values, referred to the 
coordinate system of the previous section. 

Since the coasting ellipse also preserves total energy, we 
have 

2 (R. + Yo) 2(Re Ya)? (Re Ya) 

(R. + ya) is the specified orbital injection height. 

Mathematically speaking, the problem is to find that func- 
tion a(t) which maximizes expression [4] subject to condi- 
tions [1, 2 and 4], where (R, + yz) and the first stage burnout 
characteristics are known. This problem is treated within 
the framework of the variational calculus, and an analysis is 
presented in the Appendix. We describe here the principal 
results, and a computational scheme to be employed in the 
solution. 
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Discussion of the Analytical Results 4 | 


In support of the conclusions of (2 to 5), we find the optimal 
thrust angle program to be characterized by the expression 


tana = b —at [5] 


where a and 6 are constants to be determined. Equations 
{1 and 2] which must also be satisfied, add four more con- 
stants to be found, making a total of six. 

On the other hand, we have the parameter-varied set of two 
staging values (yo) and yo(yo), plus the value of yo itself. 
Two relations at burnout are afforded by the energy equation 
[4] and the burnout transversality condition (Eq. [A-20]) 


tan @ 
tan 0 


_ 
th Qn» + ap 


[6] 


where = + and ax = 
Finally, the sixth condition is given by the staging trans- 
versality relation for a fixed staging time to) (Eq. [A-24]) 


dyo dyo Jy 
which guarantees that the complete solution is stationary 
over all neighboring choices for first as well as second stage 
trajectories. 

The solution is therefore determined in principle, and pre- 
sumably is calculable for any choice of orbit height. How- 
ever, the transcendental nature of the equations defining a 
solution precludes any direct approach. It proves more 
expedient to plot orbit height as a function of starting condi- 
tions, and to extract the required result by interpolation. For 
every choice of yo there will result a corresponding orbit height 
(R. + ya), and an apogee velocity Va, which will be maximum 
for that height. We propose to calculate (R. + y.) and Va as 
functions of yo, and to select from these curves the desired 
values. The one-guess calculation proceeds as follows. Choose 


5 Referring to Equation [6], it can be shown that a, is the op- 
tima] thrust-angle turning rate for maximum burnout momentum, 
neglecting the energy consideration [4]; conversely, a is th: 
optimal turning rate for maximum burnout energy, if momentum: 
considerations are neglected (2). 


ARS Journal! 


& 
( 
] 
( 
a 
> 
> 
7 
] 
t 
s 


J 


a value for yo. Next, guess the value of b = tan a. The cor- 
responding value of a is obtained from Equation [7], since the 
derivatives involved are assumed known, at least numerically, 
if not analytically. For the value of yo selected, z and y 
are also known. Now the integrated equations of motion are 
used (2,4) 


tan [(@ — 
tan [(a@ — ax)/2] 


= + COS In [8a] 


tan [(a@o — a,)/2] 


the question of maximality of V. reduces to determining the 
sign of the right-hand side of Equation [11]. If this latter 
quantity is negative, V. is mathematically a local maximum 


and the trajectory is to be considered optimal. 


From physical considerations, the integral term 


lb T 
— sec a ba? dt 
™ 


jp = [Yo — g(t — &)] + 


A 


(Re + yw) = [((Re + yo) + Yolts — to) — g(t» — t)?] — mp {sin a,x In 
mo 


c is the effective exhaust velocity, considered constant in the 
vacuum we have assumed. ax = tan! [a(mo/mo) + 6], and 
is interpreted as being that thrust angle which would finally 
result, were the entire second stage mass to be burned as 
fuel.6 a, y and y, obtained here must satisfy Equation [6]. 
I) this is not the case, guess another value for b and repeat 
the calculations until Equation [6] is satisfied. The 
correct values of a, y and y, thus found are used to 
evaluate M, in Equation [3] and (R. + ya) in Equation [4]. 
Finally, V. is obtained from the momentum equation 2 


M, = (R. + Ya)Va = M, [9] 


(2. + ya) and V, are plotted against yo first, and then against 
each other, so that coast apogee velocity may be obtained 
directly for any given injection height. The scheme is quite 
simple, and very rapid if an automatic computer is used. 
I:ven by hand, a few hours labor will usually suffice for the 
calculations. 


Second Variation and the Question of Maximality 


In consequence of the results derived in the Appendix, we 
find (Eq. [A-29]) 


tan | Re + dyo 2! 2! J, 
{in + 


By Equation [9], since (R. + ya) is fixed, 6M, = 6V. = 0, 
and 62M, = (Re + ya)6?Vc. Equation [10] reduces to 


tan a dio 
&v, = — b 
| p*) tan 2! 


tan @ . . 


tan Q, 
+ ay? 


Total velocity increment due to optimizing stems from three 
factors. The first represents a contribution from staging 
point errors, resulting from non-optimal use of the initial 
booster stage. The second term arises when a non-optimal 
second stage thrust direction program is employed. Finally, 
the last term is due to a difference in coasting arc lengths be- 
tween the optimal and any comparison path. : 
For the range of values considered here 


p? = (Re + yw)/(Re + ya)? ~ 1 


and the term (1 — p*)(tan a»/tan @, )is negligible in compari- 
son with unity. The coefficient of 6?V, is positive then, and 


6 We omit consideration of the infinite accelerations which 
would theoretically be produced with a vanishing mass. 
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d*y 
~ 
1, ™m 


tan [(ao + [sb] 
tan + 
tan — ax)/2] tan [(Q@ + m/2)/2) 


tan [((a, — ax)/2] + 4/2)/2] 


[see — see a] 
a 
We therefore conclude 


in Equation [11] must be negative. 
that the condition 


d*io +1 dijo] bye? 
dyo? 


dyo?_| 2! \ 


tan 


[E? + 


(i? — bye? 
(Re + 


is sufficient for a maximum V,. Inclusion of the first brack- 
eted term in condition [12] guarantees a compound maximum, 
implying that the entire trajectory, containing both boost 
phases plus the coasting are, is the optimal solution. 


Furthermore, let us now suppose 
| 
dy¢? 


dys? 


: seca 6a? ant 


m 


0 1s [10] 
— > 
f (Re + yw)? 
seca 6a? as 
[11] 


(Re + ya Re + ys) 


to be negative. Rough numerical calculations seem to indi- 
cate that the assumption is justified for several existing 
boosters, and there seems good cause for its acceptance. With 
this conjecture valid, the condition 


tan 


a [13] 
b 


[E? + (au? — >0 


becomes sufficient for a maximum. The inequality holds if 


— yo?) > 0 


and 
tan a,/tan 6, > 0 


1 
Re 
| 
) 
a 
r 
t 
1 iz 
Ss 
| 
= 
845 


both of which seem to be the case in most physically sensible The first bracketed term on the right-hand side of Equation 
instances. [19] may be written as 


Total Variation; Direct Comparison With Non- [aAyo + At» + bAg) = 


P (= + dy? 2! dyo* + dyo* 3! 
Consider Equations [1 and 2] for the special case when they 
define motion along the optimal trajectory. Let us adjoin to by virtue of the staging relations zo = o(yo), yo = go(yo) and 
them the multipliers \; and \, of the Appendix, and integrate the staging transversality condition [7]. 
the sum Ai¢: + Azd2. There results It is interesting to compare expression [11] for the secon 
= t wT variation of V, with Equation [19] for the total variation in 
—ay-—z-—ytana| = —sec a dt [14] cluding terms O[Ay,*/(R. + y)?]. The apparent discrepancy 
. ad to between the two equations represented by 
Now consider Equations [1 and 2] when they define any ad- : ‘ 
missible comparison trajectory (variables on this path are (p? — 1)AV.? 
denoted by the superscript tilde, ~). We again adjoin ), iy 2i tan A 7 
and Xo, add and integrate 
. “ disappears when higher order variations are considered. A 
-—ay—#-—g7tana| = f — (cos @ + tana sin @)dt any rate, this latter term seems to be dispensable, considerin 
to tm that both (p? — 1) and |AV.|/2p%, are usually quite small i: 
[15] magnitude. For |Aa| small compared to 1 radian and term. 
By subtracting Equation [14] from Equation [15], we find (=. 
aAy + At + Aytana| =aAy+Az+bAy| 
th to 
[16] negligible, we may write da Aa, dyo > Ayo, dy, > Ay, 
— sec a[l — cos Aa]dt Both formulas then agree entirely, i.e... = AVa. For 
practical purposes the second variation (Eq. [11]) provides « 


li good approximation to the total variation when each 6( ) is 
replaced by the corresponding A( ). 
(y—y) = Ay = + d+... When |Aa!< 1 radian, the integral term in Equation [19] 
be written 


and similarly for the other variables.’ A(_) is the total varia- 
tion, and refers to the actual numerical difference between tb T rs 
— sec a[l — cos Aaldt = — sec? a; e[(Aa?)] 


optimal and non-optimal paths. “ 
Since 
A(uv) = (u + Au)Av + vAu a +1 = where a = ai = a. €[(Aa?*)] is the mean-square value of 


_ Aa, for any given non-optimal program, weighted rather 


strongly toward t, by the term = (7'/m) cos a, which in- 
1 Au creases rapidly with time. Because of this factor, relatively 

7 ; ae eee little penalty is suffered if a small, finite time | ired i 
a a +hn) ___ little penalty is suffered if a small, finite time be required in 
practice to adjust the pitch angle immediately after staging. 
Equations [3 and 4] provide (after considerable manipulation) Generally speaking, a non-optimum program which steadily 
the relations approaches the neighborhood of the optimum is a favorable 
substitute; zero-lift programs seem to fall into this category, 

Rigs M, — inAy i. (17] which apparently explains their effectiveness (2). 


[(Re + yo) + Awl] 


Ay E | AM;? Ay? goR.? ) 

+ yo) ins Re + yo(Re + 2y(Re + yo)? yl Re + yo) 2 Re + 


Here we have assumed that terms 


etc. are negligible. 
Equations [17 and 18] inserted into Equation [16] produce 


1 tana 
— 4 (p? — 1)AV.2 + p -(1- p?< [aAyo + Azo + — — sec a[l — cos 
19 
7 Equation [16] for the total variation may also be obtained ; een = j 
by each successive ‘variation of + and forming 


the infinite sum. 
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Illustrative Example 


A convenient example is furnished by the Vanguard satel- 
lite launching vehicle, formerly in use for IGY research. 
The rocket is composed of three stages; the first two are 
sequential boost phases of the type under discussion. After a 
coasting period which carries the final stage to its apogee, a 
third rocket supplies the necessary velocity increment to 
establish an orbit. According to publicly released nominal 
data,’ the following constants were used 


A 


= 129.5sec c = 8694 fps 


mo = 149 sl. m = 37.3 sl. 


Staging data were available only for a single trajectory; these 
are 


yw =19X 10 ft = = 3889 fps 

Since parameter-varied staging information was unavailable, 
the following two assumptions were applied and seemed ac- 
curate enough for present purposes: 

1 zo + gy = const = 2(3889), for all staging points in 
the vicinity of the known one. 

2 Total energy is assumed invariant for all staging points 
in the vicinity of the known one. 

Using these assumptions, reasonable expressions were ob- 
tainable for the derivatives appearing in Equation [7]. The 
results are plotted in Fig. 2, where flight with the optimum 
second stage program is compared with two non-optimal cases, 
zero lift and constant thrust direction. The same first stage 
trajectory family and identical second stage approximations 
were used for all three curves. With the two foregoing ap- 


proximations valid, the term 

» - 
d*y 
dy»? 


dyo* 


mentioned in the previous section, is always negative. 


Appendix 


The variational method is applied, with customary nota- 
tions, in the following manner (8). From Equations [1 


and 2] 
[Al 


7 bg: = by — (T/m) cosa dba = 0 [A-2] 


Multiplier functions ,(1) and A2(t) are introduced, and the 
product F; = Ai gi, + = O is formed. The quantities 
and \, are as yet undetermined in form, and subject only to 
the restriction that each be piecewise continuous between fo 
and ¢,, the staging and burnout times. Now F; is integrated 
between fo and 


ty to T 
Fidt = pide + + Da sin 
t J to 


The terms \,6¢ and Ady are each integrated twice by parts. 


= 6% + (T/m) sina da 


[A-3] 
Az cos alba dt 


We find 
= — | + J dxdt | 
to ito to 


Equation [A-3] becomes 


8T am indebted to S. B. Kramer of Lockheed MSD for these 
figures. 
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Fig. 2 Coasting apogee velocity vs. orbital injection height for 
different thrust programs y 


+ ody — Ai — = 


t 
an an 


7 
+ Ady + [A; sin — cos dt 
ty 


Equation [A-4] must hold for any vile of 6x, dy, dz, dy at 
each limit, in particular, when all eight of these values are 


zero. Therefore 
j [Arsin — Az cos al dau [A-7] 
t 


[A-4] 


Now we choose 


must hold. 

However, of the three variables 6x, dy, da appearing in 
Equations [A-1 and A-2], one may be considered inde- 
pendent—say da. If da is arbitrary, the Fundamental 
Lemma (8) applied to Equation [A-7] guarantees that 


Ai sina — A, cosa =0 


tana = 2/\ : [A-8] 


in (to, ty). 
We have the particularly simple solutions to Equations [A-5 


and A-6] 


Ae = b — at [A-9] 
=d+et [A-10} 

from which Equation [A-8] shows 
tan a = (b — at)/(d + ct) [A-11] 


| 
ea 
= 
Or 
Ww 


to be the optimal program. This result is to be found in 
(2-5). 
Returning to Equation [A-4] in its general form, we have 


to 


to 


[A-12] 


Equation [A-12] must hold for any and all values of either 
side, so that each side must be identically zero. Consider 
the left-hand side 


=0 


+ — Arde — Ardy [A-13] 


We have, however, at point ¢, the two relations Equations 


q 
Therefore Equation [A-21] is rewritten as 


[A-23 


— dé — Ardy] — [hey — — Ars] 


However, under the assumption of the problem at hand that 


ty be fixed, Equation [A-23] becomes a 
= 0 [A-24 | 


The second variation of M, is now evaluated. We form 


*° dyo 


di» 


dyo 


= + = 67% + + seca ba? = 0) 


where \,; and ), are defined as previously. Performing int« 
grations analogous to those of Equation [A-3] ff., we obtai: 


[3 and 4]. Since M, is to be a maximum, 6M, = 0, and the 
first variations of these expressions may be written as 4 T 
y — — — seca da7dt A-25 
= 6M, =0=(R,+ +- Yo [A-14] ty 2! to m 
= + goRe bys which reduces to 
(Re + ys)? @ bye? 
From Equations [A-14 and A-15] we have 2! 
(Ra 2! Jum mm 
; 1 goR? Likewise, the second variations of Equations [3 and 4] are 
= + E [A-17] found (using Eqs. [A-16 and A-17]) to be 
these latter two relations into Equation [A-13], it &M, — (Re + = — Sys" [A-27 
M, 1 1 ( y ] 
— — = + 4° + — - [A-28 


a, | 6M, 
1 
| tan Re + y 
becomes 
Avs (: _ oh? ) (s 
——}) - r by. — 
Avda, = 0 [A-18] 
6x, and 6y, are arbitrary, so that we find 
[A-19] 
a tan a, gol.” ) 
= A-20 


Equations [A-5 and A-6] are homogeneous in A; and ), and 
their derivatives, we may assign any arbitrary value to either 
multiplier at any point in (t, &). We choose Ay = d = 1. 

Now let us return to Equation [A-12]. The right-hand side 
IS 


Note here that the operator 6( ) refers to changes from one 
comparison are to another at a constant time t. However, 
the staging time t may itself vary from one first stage arc 
to another, depending upon the particular staging condition 
chosen. In the general case then, the complete displacement 
resulting from following the staging point from one arc to a 
neighboring one is 


d( 


hody — — = 0 [A-21] 
to 


~ 


= ) + [(d/aty( 
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t 
[E? + (2? writ 


Equations [A-26, A-27 and A-28] are solved for 67, fe 


Oye? 
(Re + ya)(Re + yo) 


= (Re + yo)/(Re + yu). 


where tan = and p 
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NHE MENTION of an Earth-period or 24- shin satellite 

immediately brings to mind the often discussed and much 
publicized vision of some type of vehicle which will remain 
suspended over a point on the Earth’s equator as if connected 
to some form of “sky hook.” This is possible if the following 
conditions can be met: 

1 The orbital plane is exactly coincident with the Earth’s 
equatorial plane. 

2 The orbit is exactly circular. 

3 The altitude is 22,240 miles, corresponding to a period 
exactly equal to the Earth’s rotating about its axis (i.e., one 
sidereal day). 

4 The satellite is not subjected to perturbing forces which 
might alter its path. 

These conditions would be impossible to establish with per- 
fect precision, which brings up the question of what happens 
in the imperfect or general case applying to all Earth-period 
satellites. 

Let us assume that the satellite period is perfectly equal to 
one sidereal day. This assumption defines the semimajor 
axis of the elliptical orbit, or the radius in the circular case, as 
equal to 6.61 Earth radii (see Appendix A). The period, 
being a function of only the semimajor axis, depends on the 
initial conditions of velocity magnitude and attitude at 
launch. The remaining orbital parameters depend on velocity 
magnitude and orientation, altitude, position and time of 
launch. The first question of interest is the relative ground 
track, which is described by a line between the Earth’s center 
and the satellite, intersecting the Earth’s surface. Six orbital 
elements define the position of any satellite in space at any 
time. Since in the present case we are interested only in the 
relative motion of the satellite with respect to some arbitrary 
point on the Earth’s surface at any arbitrary time, two of these 
elements, defining space orientation and time of passage, may 
be ignored. The four remaining elements which will define 
the relative motion and time are: 


P = period (23"56” 4.091°) 
= eccentricity of the elliptical orbit 
7 = angle of inclination of the orbital plane with the Earth’s 


equatorial plane 
angular longitude of perigee measured from the line of 
nodes 


(For the “sky hook,” e = 


ll 


0, 7 = 0, and w is meaningless.) 


Inclined Circular Orbit (e = 0, i + 0) 


Now consider a case where the orbit is circular; the period 
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_ An Earth-period satellite is a satellite which has an orbital period of revolution exactly equal to 
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the Earth’s period of rotation on its axis (i.e., one sidereal day). As originally conceived, the term 
implies an idealized equatorial and circular orbit where no relative motion exists between the 
Earth and satellite. This paper discusses the type of relative motion which will exist in a practical 
case where the orbital plane is inclined, the motion is elliptical, and perturbations due to the moon, 
sun and Earth’s oblateness are considered. 
various different orbital parameters (using exaggerated values) on the satellite ground track. Some 
possible practical uses are mentioned, and it is suggested that the use of an inclined orbit may have 
some merit for observation purposes. 


Several illustrations are given showing the effect of 


is equal to one sidereal day; ~ orbital plane is inclined at 
some angle 7 to the Earth’s equatorial plane; and no per- 
turbing forces exist. That is: 


= 23"56” 4.091° 
a = 26,200 statute miles — 
e = 0 


By a simple transformation of spherical coordinates (given 
in Appendix A) the relative Earth motion, or ground track, of 
the satellite may be described by the deviation in longitude 
and latitude from a reference point on the equator. Since the 
period is the same as the Earth’s rotational period, the satel- 
lite must return to the same point from which it started over 
the Earth at the end of one period. For convenience, the 
latitude and longitude of this point may be taken as the 
reference. A plot of the deviation from the equinoctial ref- 
erence point is shown in Fig. 1 for several values of i. The 
ground tracks representing the relative motion are shown to 
be “figure eight”’ paths over the Earth’s surface. The lobes 
of the “eight” shrink to a point for 7 = 0, and expand to 
occupy an entire hemisphere for a polar orbit. 

Fig. 2 shows a plot of the maximum longitudinal deviation 
as a function of inclination ear 7 for a circular orbit. 


The ground tracks of elliptical orbits are not so easily 
visualized as in the preceding case, since two additional 
parameters shape the apparent path (P and 7 having been as- 
sumed). These are: Eccentricity e and angular longitude 
of perigee w, measured from the ascending node. In Fig. 3, 
Q defines the line of nodes, w is the angular longitude of 
perigee, v is the angular displacement of the satellite from 
perigee at time ¢ (i.e., time is assumed to start from perigee 
crossing), and uw is the angular displacement measured from 
the line of nodes. It is convenient to take the Earth meridian 
at perigee crossing as a reference, corresponding to t = 0. 

The general effect of eccentricity variation is to bend or 
distort the lobes of the basic figure eight pattern, as shown in 
Fig. 4. For these cases, the line of apsides has been chosen 
to lie in the equatorial plane, i.e., w = 0. Since north-south 
symmetry exists, only one lobe (north or south) is shown, for 
clarity. 

The effect of rotating the line of apsides (i.e., the angular 
longitude of perigee w) is to cause a displacement in latitude 
of the crossover or nodal point as illustrated in Fig. 5. The 
dashed curve shows that the nodal point is shifted from the 
reference perigee meridian and the equatorial plane in a 
similar figure eight locus as w is varied. 
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Perturbing Forces 


In the preceding case, no perturbing forces have been as- 
sumed. Because a 24-hr satellite is at a distance from the 
Earth of 6.61 Earth radii, the gravitational pull of the moon 
and sun may be expected to be more significant than with 
lower altitude satellites. However, the more significant 
perturbing force is still that due to the oblateness of the 
Earth. 

The effect of atmospheric drag has been assumed to be 
negligible. The sun and moon will exert an attraction on 
both the Earth and the satellite. Thus, both will accelerate 
toward the sun and moon. The difference between the ac- 
celeration of the Earth and the acceleration of the satellite is 
the disturbing acceleration which perturbs the orbit. In 
Appendix B, the disturbing accelerations due to the sun, the 
moon and the Earth’s oblateness have been calculated using 
a first-order approximation. For a low altitude satellite, the 
disturbing forces are about the same as the tidal forces exerted 
on the Earth, but the disturbing acceleration on a 24-hr 
satellite due to the sun is about 200 times as great as this 
tidal acceleration. Similarly, the moon will cause an ac- 
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Fig. 2 Maximum departure from meridian for circular orbit 


— (e = 0) ae 


celeration 230 times greater. These perturbations will cause 
the plane of the orbit to regress and the path to depart from 
circularity or ellipticity. In the worst case, where the Earth, 
satellite, moon and sun are in line, or in conjunction, this de- 
parture will be about 3600 ft. The combined regression rate 
due to sun and moon is 1.44 X 10~® radian per sec, resulting 
in one revolution of the orbital plane about the ecliptic in 139 
years. The maximum regression rate due to the Earth’s 
oblateness will be 2.74 X 10~® radians per sec about the celes- 
tial pole, which, when added to the rate due to the sun and 
moon, will result in approximately one revolution of the or- 
bital plane in 49 years, or 7.5 deg per year. Thisisa worst-case 
value, since the period will be greater for higher inclination 
angles owing to the variation of the rate due to the Earth’s 
oblateness with cos 

Summarizing, the effects of perturbations on a 24-hr sate!- 
lite are quite small relative to those encountered by a low alti- 
tude satellite. The most significant result of these perturba- 
tions will be to cause a westward drift of the ground track of 
the order of 7.5 deg per year. In theory, this drift might be 
removed by lowering the period to slightly less than a sidere:! 
day, creating a counteracting eastward drift. 


ARS JourNAL 


> 
4 
> 
] 
( 
- ( 
~ 
> 
inclination vari 
~ 
> 
- - @ 
] 
f 
. 
a4 
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u=wty .* 
: w = Angular longitude of perigee 
~—_ v = Angular displacement of satellite 
from perigee at any time, t 
' 


u = Angular displacement of satellite 
from line of nodes 


Fig. 3 Spherical coordinate system 


Possible Applications * 


By definition, a 24-hr satellite is a satellite selected to - 


minimize relative motion with respect to the Earth As we 
have seen, this motion may be nearly zero for the equatorial 
orbit, or may range over an entire hemisphere for the 90-deg 
inclined polar orbit. The zenith positions of all are confined 
to one hemisphere of the Earth. The equatorial-plane 
satellite will only be within line-of-sight for points on the 
Earth within a globular sector +81.3 deg latitude and +81.3 
deg longitude (relative to the satellite meridian); whereas the 
polar satellite will be within line-of-sight of almost every 
point on the Earth for at least one time during the period. 
This implies that for observation purposes, an inclined orbit 
might be more desirable due to its figure eight wandering 


ground track. By properly phasing two or more satellites, | 
the northern (or southern) latitudes and, in particular, the | 


polar regions may be kept under constant observance. This, 
of course, is not possible with the equatorial orbit, because of 
the 81.3-deg limitation. 

For purposes of communication, it would seem that the 
equatorial orbit is more desirable in order to minimize the 
Karth-antenna pointing and tracking problem. 

The ability of an optical sensor, such as TV or photography, 
to resolve ground detail at high altitudes depends greatly 
upon the contrasting luminosities and background. Because 
of reflection and light scattering from layer discontinuities 
above the tropopause, the brightness and contrast of ground 
areas will be veiled by the higher level returned light. How- 
ever, due to the low atmospheric density, there is no 
reason to believe this effect would be any worse for a 22,000- 
mile orbit than for a 500-mile orbit. Accordingly, sufficient 
aperture and focal length should provide comparable ability 
to resolve ground detail. 

At 22,200 miles altitude, 100 ft on the Earth will subtend 
an angular are of about 0.1615 sec. The use of Dawes’ em- 
pirical criterion for aperture and resolution results in a re- 
quired optical aperture of 28 in. to resolve two points sepa- 
‘ated by this distance. In terms of wave length, this aperture 
‘s roughly 1.42 X 10°. Accordingly, for the same 100-ft 
resolution, an infrared sensor would need an aperture of 112 
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in., and a 3-cm radar would require 42.6 km. For higher val- 
ues of resolvable angular arc, the required aperture is in- 
versely proportional. 

These values are given to suggest the practical limitations 
involved in any application for detailed ground area observa- 
tion. Such large aperture sizes may not be required to detect 
point sources of light (or heat). In that case, detection will 
depend upon the intensity of the point source, upon back- 
ground contrast and upon the square of the aperture. 
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a 
where 

k = 0.0743651 

P = 1436.067 min 


Circular Orbit (e = 0, i + 0) 


Let n equal the Earth’s angular rotation rate about its axis. 
Then n = 27/P (where P = 25* 56" 4.091*). For a circular 
orbit, the angular rate uw (see Fig. 3), is constant and by defini- 


tion 
cm 
[1] 


Now, assuming a spherical coordinate system as shown in Fig. 
3, where 


u=nt 


@ = angle between r and the projection of r on the Earth’s 
equatorial plane 

radius vector from the center of the Earth to the 
satellite 

angle measured in equatorial plane between the line of 
nodes and the projection of r on the equatorial plane 


ll 


ll 


Then, by trigonometry 
cos As = cos u/cos 
id 
As = cos cos u/cos@ 
sing = sinusini 
cos @ = 


As = cos~ cos u/+/1 — sin? u sin? 7 


— sin? u sin? i — sin? sin? 7 


Designating the angle in the Earth’s equatorial plane 
through which a meridian moves in a time t as Ag = nt, then 
the longitudinal deviation Dy, of the satellite from any given 
Earth reference meridian is 


Dy = Ag — As = nt — cos~ cos (nt)/+/1 — sin? (u) sin?i [2] 


where, from Equation [1], u = (nt). 
By trigonometry, the deviation in latitude is given by 
D, = = sin~' sin u sin7 [3] 
Fore = 0 


@ = sin™ sin (nt) sinz 


Differentiating Equation [2] with respect to (nt) and setting 
it equal to zero, a maximum occurs when 


snAg=+ 1/1/(1 + cos?) 


Dmax = cos~! 2 +/cos i/(1 + cos 7)? [4] 


and 


Elliptical Orbit 


The analysis is essentially the same as in the circular case, 
see (1), except that instead of u = nt 


u=wt+v=w + tan-“+/1 — e?sin E)/(cos E —e) [5] 


where £ is the eccentric anomaly defined by Kepler’s equa- 


2?The reader is referred to reference (1) for explanation of 
astronomical terms and notation used herein. 
3 Numbers in parentheses indicate References at end of paper. 
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tion 


E—esinE = nt 


Therefore, Equations [2 and 3] are the same except that wu is 
given by Equation [5]. 


The following approximations have been given by Lyman 
Spitzer Jr. (2), and are based on the methods of Hill (3), Brown 
(4) and Brouwer (5). They assume atmospheric drag is 
negligible, the satellite orbit is nearly circular, and perturba- 
tions due to other planets are insignificant. In the following 
examples, the maximum values of the perturbing effects arv 
calculated to a first-order approximation. 

Assuming that the Earth, satellite and the sun (or moon) 
are in conjunction and that an z-axis is taken along a line 
joining the three bodies, the acceleration of the satellite rela- 
tive to the Earth is given by 


_ GM, GM, 
Rs? 
where 
R= distance between Earth’s center and satellite 7 


R; = distance between sun and Earth 
= mass of the Earth 

Ms; = mass of the sun 

Gy = gravitational constant 


The disturbing acceleration is the sum of the last two terms. 
Designating y as the ratio of the disturbing acceleration to the 
Earth’s acceleration, and expanding (Rs — R)? 


M, (R\3 3R\ 


(neglecting the higher order terms). 
numbers 
Ms; = 1.98 X gm 
M, = 5.97 X 1077 gm 
R = 26,200 miles : 
R; = 93 X 10'miles 
and for the moon heh fl 
Ry = 2.39 X 10° miles 
My = 7.32 X 10% gm 


y (24-hr satellite) = 1.48 x 10-5(sun) 
= 3.72 1075(moon) 

(500-mile satellite) = 0.73 1077 (sun) 
= 16 X 1077 (moon) 


an approximation to G. W. Hill’s equation for departure from 
circularity p is given by 


= 2575 ft (moon) 


At conjunction of Earth satellite, moon andsun 
Pmax = 3600 ft 


The regression rate of the line of nodes referred to the 
ecliptic plane and the ecliptic pole is given by 


| 
Total = 14.25 x 10~” radians/sec (at conjunction) 


The regression of the line of nodes due to the Earth’s oblate- 
ness referred to the equatorial plane and the celestial pole is 


approximated by 
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2 3 
Q = 3 = 4.05 X 107 (sun) (» = 


10.2 (moon) 


Ee 


3 
= 6.61 Earth radii 
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where 
3(C — A)/2M = 1.64 x 
by measurement, and a = mean radius of the Earth. Ac- 
cordingly 
Q = 2.73 X 10-*radians/see 


due to the Earth’s oblateness. The regression of the line of 
nodes due to the sun and moon has been calculated for the 
ecliptic plane. This regression rate referred to the equatorial 
plane is then 


Q = 1.425 X 10-* cos 23.5 deg = 1.31 X 10~® radians/see 


Wire 


dynamic performance and control. 


its characteristics. 


f URING the last 10 years several authors have considered 

the application of nuclear energy to propulsion systems 
for space missions. Although a number of exotic propulsion 
schemes have been proposed, most of the reported work has 
been concerned with the more straightforward approach of 
solid structured heat exchanger systems (1,2). Since such 
systems are intrinsically temperature limited, low molecular 
weight working fluids, such as hydrogen, have generally been 
considered (3,4) for high performance. 

It has been shown that the nuclear rocket engine offers a 
substantial performance advantage over chemical systems; 
however, many engineering and operational problems must 
be solved if this advantage is to be realized. McGrath (5) 
has pointed out the need for a systematic development plan in 
order to achieve this goal. An extensive subsystem develop- 
ment program must precede engine hot firing, because of the 
hazardous and perhaps infrequent nature of nuclear engine 
hot firings. However, experience with more conventional 
systems has shown that many component problems arise from 
the system environment and require solution both by real 
and simulated system operation. In order to accomplish 
simulated system operation, it is necessary to properly model 
system characteristics. The simulation hardware may con- 


Presented at the ARS 13th Annual Meeting, New York, N 
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methods of simulation are discussed. 
applications, is used to develop system dynamic performance criteria. 


Dynamic Analysis of a Nuclear 
Rocket Engine System 


The total regression rate for small inclinations is then 4.04 
10~® radians per sec or one revolution in 49 years. 
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An important consideration in the design of an operational nuclear engine is that of the engine 


Chemical systems have generally been designed for fixed cali- 


bration rather than automatic control of performance. 
signed to be sufficiently sluggish so that manual control is feasible. 
neither approach is suitable for a nuclear rocket engine system. 

exchanger-type nuclear engine is described together with nonlinear dynamic equations describing 
Distributed parameter effects in the heat transfer processes are presented, and 
The analysis, together with considerations of nuclear vehicle 


Nuclear reactors have usually been de- 
This paper will indicate that 
A typical solid structured heat 


sist of mechanical, hydraulic or electronic models; however, in 
all cases a realistic mathematical system model must exist in 
order to permit transformation of real system characteristics 
into a simulator performance specification. 

The purpose of this paper is to present one such mathemati- 
cal model which has been used to study system dynamic and 
control characteristics with an electronic analog computer. 
System operational and control problem areas which have 
been developed during the study are discussed. 


System Operational and Control Requirements _ 


Applications 

In order to determine some of the problem areas which 
must be considered in the design of nuclear rocket systems, it 
is advantageous to speculate on possible applications of nu- 
clear vehicles. Table 1 illustrates vehicle performance re- 
quirements for typical missions. The two lower specific im- 
pulse values (J, = 300 see and J, = 400 sec) are in the range 
discussed by Tor mey (6) for the more energetic chemical pro- 
pellants. The higher values are in the range calculated by 
Rosenblum (4) for possible nuclear working fluids. Evi- 
dently, Earth orbits can be achieved with a single-stage 
nuclear vehicle, but for longer journeys one might eject the 
booster after reaching orbital or escape velocities. However, 
for less energetic space maneuvers, such as lunar capture, one 
would probably not wish to discard a nuclear engine after 
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time is given by (8) 
Table 1 Ideal mass ratio requirements for space maneuvers, ; ts 
Ideal velocity R F/W 
increment 
Maneuver (15), fps | Average specific impulse, sec As shown in Table 1, mainstage will last for less than 10 or 


20 min. Longer burning times occur if the engine is a sus- 
tainer type, which does not overcome the full force of gravity. 

Curorr. Engine thrust decay or cutoff requirements dc- 
pend on the type of mission. However, fission fragment afte- 
heat defines the procedure of reactor shutdown, because by 
inserting all of the control elements, immediate neutron lev '] 
depletion is impossible due to background radiation. An 
Earth takeoff, with booster ejection after reaching some pr-- 
scribed velocity, requires a rapid thrust decay for separatio .. 
Hence propellant flow rate must be reduced to zero, and t! e 
fission fragment afterheat will vaporize the engine in a fev 


300 400 800 1200 + # 1600 
enter 350 mile Earth 
orbit 27,000 16.4 8.13 2.85 2.015 1.692 
enter lunar transfer 10,000 2.82 2.18 1.476 1.295 1.215 
lunar capture and 
landing 19,100 2.57 2.02 1.425 1.266 1.194 
enter Mars transfer 11,600 3.32 2.46 1.57 1.35 1.253 
Mars capture 6125 1.89 1.61 1.27 1.172 1.127 
minimum energy to 
Mars orbit and re- 
turn to Earth or- 
bit 63,500 710 104 11.6 5.16 3.42 


every firing. Not only would such operations be expensive, 
but also inefficient from an overall vehicle performance 
standpoint, since the heavy nuclear engine would not be used 
to its full capability. 

From these considerations, it is clear that for some types of 
missions, it is desirable to be able to start the engine in space 
and cool the reactor after firing. Proper functioning of pro- 
pellant pumps also requires close control of propellant _ 


cal properties on extended journeys. 


Engine Operation 


In order to typify the of. nuelea ear ‘engine, 
consider the possible operating sequence described here. 

IniT1AL Reactor Srartup. After the reactor is critical, 
its power must be increased by many decades before a tem- 
perature rise is felt. In the range of low neutron flux levels, 
stabilizing feedback reactivities are absent, and one would 
proceed slowly with startup over an interval of perhaps 10 
min. For example, to increase power by 10 decades from 
source strength on a reactor period (e-folding time) of 20 sec, 
about 8 min is required. The reactivity required is a rela- 
tively safe value of 30 cents (7). During this time neutron 
level indications are required; however, the control rod posi- 
tion could be manually adjusted for booster and manned ap- 
plications. 

Reactor Warmup. As the reactor approaches the rated 
power range, a period of a minute or more is required for the 
fuel elements to reach the operating temperature. The exact 
time would be fixed by reactor thermal time constants and 
thermal stress criteria. Neutron detectors (ion chambers) 
would supply signals for the reactor warmup rate control, and 
reactor temperature indications would initiate the controlled 
engine transition phase of operation. 

EneineE Transition. As the core nears mainstage tem- 
perature, coolant must be supplied. The turbopump would 
be started at a low power condition, perhaps from a chemical 
gas generator. After propellant values are opened, the reactor 
power would be programmed to mainstage in close synchroni- 
zation with the propellant flow. To conserve the propel- 
lant, this would be done as quickly as would be practical. 
During this transient and most critical period of reactor opera- 
tions, accurate and reliable neutron level, propellant flow and 
temperature indications are required for control. 

MarnstaGE Operation. After achieving rated power, 
steady-state operation will not be immediately achieved due to 
long term thermal transients in the structure, fuel burnup, 
poison buildup, etc. Depending on the relative values of 
component tolerances and internal feedback effects, automa- 
tic control may be required during mainstage operation, but 
this system may be very simple compared to that of the start- 
up period. Mainstage duration is less than J, in seconds be- 
cause for a constant thrust, Earth takeoff vehicle, burning 
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thermal lag times. The second and more difficult form of 
cutoff, if one desires to re-use a nuclear engine, requires pro- 
pellant to continuously cool the reactor until power is sufii- 
ciently low so that another form of cooling, such as therm il 
radiation, can be employed. In order to conserve propellant, 
it would be desirable to throttle the propellant flow so as ‘o 
keep the exit temperature as high as practical. In the latt-r 
case, the expellant would not entirely be a weight penalt., 
since thrust would be generated w hich might be be used for 
vernier trajectory correction. ~ 


Control Requirements 


From these general considerations, some of the control 
hardware requirements may be outlined. 

INSTRUMENTATION. In order to assure satisfactory engine 
operation during start, mainstage and cutoff, transducers are 
required to independently measure the following variables: 
Neutron level (power), temperature and reactor flow rate. 

Obviously, many other measurements are also required in a 
development program. 

ControL Computers. An integrated engine control 
system requires a control computer which accepts and properly 
operates on the subsystem performance information in such 
a way as to satisfactorily synchronize component operation. 

Finat Controut ELeMEnNts. In our example, the control 
computer generates signals to properly position the following 
subsystem control elements: 

1 Reactivity control elements. These might be moving 
poison rods, or flows of poison or moderating fluids within the 
reactor. 

2 Turbine control valve. The gross buildup of propellant 
flow is here determined by programming the pump input 
power. 

3 Main propellant values. Initial priming of the pro- 
pellant passages (with which may be associated a positive 
reactivity insertion) as well as the proper switching from 
warmup to transition is controlled by the main propellant 
valves downstream of the pump. 

The requirement that the hardware must operate reliably 
in an environment of high nuclear radiation and heat genera- 
tion poses several novel design problems with even the sim- 


plest of controls components. 

Engine Concept 

An illustrative nuclear engine system schematic, such «s 

has been discussed by Rosenblum et al. (4) and McGrath (5) 
is typified in Fig. 1. It is similar in concept to a chemical 
engine except that the combustion chamber is replaced by a 
nuclear reactor. The working fluid, stored in a low pressure 


tank to save weight, is pumped.to high pressures and initial!y 
used to regeneratively cool various engine components which 
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System Representation 
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are heated due to the absorption of thermal and nuclear radia- 
tion fluxes. The fluid then enters the reactor heating chan- 
nels and finally exhausts at supersonic velocities through a 
conventional de Laval nozzle. A portion of the heated work- 
ing fluid is either bled off from the main stream or passed 
through a separate heating region in order to provide hot 
gases to power the turbopump. A summary of lumped pa- 
rameter equations used to describe these processes is shown 
in Table 2. 


Reactor Neutronics 


The configuration of a solid structured, gas cooled reactor 
suitable for rocket propulsion will generally be established 


Table System equations 
I NEUTRONICS 


TANH I+ S/we 


8 P= (Te—Tro) + (Ma-Mro) 


2-7 Cr 


9 Qy= V E,)N 


IT HEAT TRANSFER ( ONE LUMPED PARAMETER SECTION) 
10 M, Cy (1, -7,) i. (T,+T2) 


1) Qs= h (Tw—0.5T, —O.5Te) 


IZ REACTOR — THRUST CHAMBER PRESSURES 
M, 


SM, 
16. 


IZ FEED SYSTEM 
PUMP ; 


17 Mp= M, +My 
18 Pp = Po+dp K, Nz + Ke Me Npt ks Me 


I9 Le= dp Ky MpNpt 
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from heat transfer and structural design considerations. 
However, if the thrust level is low, then reactor criticality 
also must be considered. The steady-state design concept 
used depends strongly on the performance requirements; 
hence, many trial and error designs are required before the 
final concept is evolved. Fortunately, the kinetic features of 
the fission process are sufficiently similar in many reactors 
so that a fairly simple model produces useful results. 

The speed with which a reactor can change power level is 
determined primarily by the magnitude of reactivity changes 
produced by temperature, density, control poisons, and by 
the yield and decay constants of certain fission fragments 
which during decay emit “delayed” neutrons. The delayed 
neutrons permit power changes sufficiently slow to achieve 
control. 

The reactor kinetic equations include the neutron-level 
equation (7) 

dN p- 


—= + UC: + 8 
dt i=1 
and the precursor concentration equation 
oc; iN orc; 
= BiN AC; + 
ot l oy? 
rate of change production radioactive diffusion 


of concentration rate decay loss (10) 
These equations imply the following assumptions: 

1 The rate of change of neutron level is the sum of the in- 
crease of prompt neutrons produced in fission and the buildup 


of delayed neutron emission from fission product decay. 
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Nuclear rocket engine system schematic 
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2 The spatial distribution of power is independent of time; 
only the amplitude changes. 

3 Neutrons may be assigned averaged properties, and 
delayed neutrons have the same properties as those produced 
directly in fission. If the delayed neutron parameters are 
adjusted to the energy spectrum of a particular reactor, the 
latter assumption need not be made. 

4 A fraction of the delayed neutron precursors diffuse 
from the fuel elements at elevated temperatures and are 
immediately swept out by the propellant stream. 

Power transients are of principal interest in many dynamic 
studies, since temperature changes are relatively slow. Thus, 
it is instructive to consider the case where the diffusion ‘‘con- 
stant”’ D is either constant or varies slowly with time. 

Then, taking the Laplace transform of the precursor con- 
centration equation with the initial conditions equal to zero 


@C(y,S) A+S8 BiN(S) 
Ci i(y, S 


Ci(y, = 


The solution is 


(i + 8) + sy 


At the boundary, the eiinatiiilleds is zero (9). 


_ Bi N(S) 
( 
exp ( D y}] + exp { — y 


; Averaging the precursor concentration over the region 


E aWVs ( D2 


+ 8) 

Although this transcendental form is not suitable for elec- 
tronic analog simulation work, its exact frequency response 
plots can be obtained by digital computer techniques. The 
results can then be approximated by a series of lag-lead-lag 
approximations which converge faster than a direct applica- 


tion of finite difference methods. 
Feedback Effects of Reactivity 
The variable principally used for the discussion of neutron 
kinetics is the reactivity. In the present definition it has the 
units of dollars, where 1 dollar represents the fraction of 
neutrons arising from fission product decay. Essentially, it 
measures the excess of neutrons produced over that required 
to maintain a chain reaction. As long as the reactivity is less 
than 1 dollar, changes of neutron level can be brought about 
sufficiently slow (since it is necessary to wait for precursor 


decay) to achieve control. If the net reactivity input to the 
reactor exceeds 1 dollar, the power will rise very rapidly, and 


controller response requirements become unreasonable. 
Hence, it is necessary to limit rapid changes to less than this 
value. 


Origin of Reactivity Changes 


In addition to the precursor loss, reactivity changes arise 
from four principal phenomena in a nuclear engine. 

ConTrot ELEMENTS. These may be in the conventional 
form of movable neutron absorbing materials (boron or cad- 
mium) or fissionable or moderating material. Other possi- 
bilities include injection of poison into the working fluid, 
variation of the reflector and circulation of poison material in 
the core. 
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Gross TEMPERATURE CHANGES IN MODERATING MATERIAL. 
Changes in moderator temperatures affect criticality by 
modifying the neutron energy spectrum, thereby altering 
fission cross sections. Due to the high temperatures required 
in a nuclear rocket reactor, this effect is of major importance 
during startup. 

PROPELLANT Density CHANGES IN CoRE AND STRUCTURE, 
Three of the possible nuclear working fluids which have becn 
considered (4) are rich in hydrogen, a good neutron modera- 
tor. In the liquid state, hydrogen, ammonia and water have 
hydrogen atomic densities differing only by a factor of 2 or 3. 
Since light water moderated reactors may be critical at sizes 
as small as 12 in. (10), it is clear that the introduction of 
small amounts of dense propellant into a nuclear rocket reac- 
tor, which may be undermoderated to conserve weight (2) 
can have an important reactivity contribution. 

THERMAL EXPANSION AND PRESSURE EFFECTS ON REACTOR 
Geometry. Neutron leakage from the reactor boundary is 
altered by changes in geometry, thus altering criticality. 

These reactivity effects are summarized as linear coefli- 
cients in Equation [8] of Table 2. In general, the effects are 
nonlinear and must be evaluated from multigroup reactor cal- 
culations or critical experiments. Even then, the values can- 
not be sufficiently well relied upon to predict a reactivity vs. 
time curve which would yield a proper power buildup at 
start. Hence, the programming of an automatic power con- 
trol system is required to automatically define the proper 
reactivity sequence. 


Heat Transfer 


The difference between chemical and nuclear engines is most 
strikingly illustrated by the distributed parameter effects of 
heat transfer. In a chemical engine, energy is imparted to the 
propellant gases in times measured in milliseconds (8). The 
dynamic effect is frequently neglected or only crudely ap- 
proximated in controls analysis studies, since it is of high fre- 
quency compared to many control element responses. In a 
nuclear engine, on the other hand, energy transfer to the pro- 
pellant may take place over a region of several feet (4). 
Thus, in order to obtain a temperature rise in the propellant 
gases, it is necessary to wait until a large mass of solid material 
has warmed. Although large and rapid reactivity changes 
may cause violent power surges to occur under malfunction 
conditions, damage to the reactor fuel elements generally will 
be associated with the resulting temperature rise. Assuming a 
reactor can be designed to possess a negative temperature 
coefficient, then an important parameter in stability and con- 
trol is the thermal lag times of the fuel elements, that is, the 
time required for temperatures to reach 63 per cent of equilib- 
rium values after a step change in power. In deriving this 
parameter, the following assumptions are employed: 

1 Fuel element thermal conductivity is constant in time. 

2 One-dimensional heat transfer, that is, temperature 
gradients are small except in the direction of heat flow to the 
propellant. 

3 Uniform heat generation in the radial direction. 

4 Single phase, turbulent and subsonic flow of propellant 
through the reactor core. 

The following basic equations are thereby implied: 
element heat transfer and enthalpy storage (11) 
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_ Coolant heat transfer and enthalpy storage (12) 


oT, 
Coolant mass storage (12) ~~ 
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Coolant pressure drop and inertia (12) 


M? 
Ox ( gd,A ) 


Gas law (12) 


In order to eliminate the space dimension in the fuel 
element heat transfer and enthalpy storage equation and 
render it more suitable for simulation studies, the tech- 
nique of iterated Laplace transform will be employed. To 
simplify the treatment, initial conditions are neglected; this 
is equivalent to taking perturbations about a steady state. 
Transforming the time domain 


@T(y, duCu 8) = _ QAS) 
dy? Ku Ku 


Since heat generation is independent of y, the solution is 


8) = 22) + ¢, (V/aS y) + C2 exp sy) 
where 


a = dyCy/Kyu 


Assuming boundary conditions of heat fluxes Qa and Qs 
leaving the surface y = 0 and y = b respectively, then 


Qu/Ku = C, VaS — Vas 


—C; V a8 exp (Vas b) + ae — 
exp V a8 b) 


Solving simultaneously and integrating over y to obtain the 
average element temperatures 


Qe Ky = 


The wall temperatures are then given by 
+ Hie 
3Ky 


= 
where 


H 


H, = (6/78) [1 — WrS esch (~/78)] 


Both transfer function gains are unity at low frequencies 
and they begin dropping at 10/7 radians per sec, where 


T= dyCmub?/Ky 


For example: Taking 6 = 0.1 in., which is in the range sug- 
gested by Levoy and Newgard (2), then for graphite 


= 0.2 sec 


As a later example shows, the overall thermal time sae. 


stants may be appreciably greater than this value; these 
transcendental transfer functions may be set equal to unity in 
some cases. 

Since the coolant heat transfer and pressure drop and gas 
law equations are nonlinear, it is not possible to separate the 
space dimension; hence, one must use finite difference tech- 
niques in treating the axial temperature distribution. If 
the steady-state temperature distribution is known, it is 
frequently possible to obtain a method (such as logarithmic 
averaging) which converges more rapidly; however, the 


arithmetic averaging of properties in a region has been found 
the partial space 


suitable for dynamic studies. Here, 
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; 0.5 05 


°*. 
= (3/78) coth (~/78) 
and 


derivatives are replaced by finite differences, and time 
derivatives are taken on average properties. For this reason, 
the simplified, lumped parameter equations of Table 2 repre- 
sent only one increment into which the reactor is mathe- 
matically separated. In general, several regions are required 
to satisfactorily represent the distributed parameter heat 
transfer effects. 

In order to obtain an order of magnitude estimate of the 
overall thermal time constant, we use Equations [10 to 13] of 
Table 2 to obtain a linearized relationship between perturba- 
tions in reactor power level and temperature. The effect of 
fuel element temperature gradient in Equation [13] is also 
neglected. Perturbing Equations [10 to 13], we have 


— 0.57; — 0.5T,]\ T, 


6T, _ 1 ( Qo 
Ts S dyCy I~. Q: 


Solving simultaneously 

— = — — | /(1+7,8 

Tp Tp ) Q: 2T Mr 


= Os 


This response is that of a simple lag; that is, temperatures 
reach 63 per cent of final steady-state values at rp seconds 
after a step change in power. 

Taking, for example, a graphite element at 4000 F and a 
power density of 100 megawatt/ft* which is in the range men- 
tioned by Levoy and Newgard (2), then the thermal lag one 


might expect is 
(103) (0.5) (4000) (1055) ; 
Tp = = 2.2 sec 


108 


By comparison, the Calder Hall power reactor (2) en at 


11.5 kw/ft? and 7, — 7; = 353 F which yields - 


Tp = 1330 sec 

The need for automatic control in operating a nuclear engine 
thereby becomes obvious: In 20 min a human operator can 
easily (and reliably) set operating temperature, but in a few 
seconds, only a go-no-go type decision is possible. If the 
reactor design possesses a negative temperature coefficient, 
then a system this fast will normally be very stable at main- 
stage due to rapid negative reactivity changes caused by tem- 
perature variations; however, since the thermal time con- 
stant is inversely proportional to power, the temperature feed- 
back needed for stability will be slow and therefore ineffective 
at low power levels. 


~ 22 min 


_ Feed System 

The basic feed system components illustrated in Fig. 1 are 
the pump, turbine, hot gas throttle valve and main propellant 
control valves. Normally, the latter valves are required only 
to properly synchronize initial flow buildup. At rated power, 
the hot gas valve could be used as a thrust regulating device, 
if desired. The principal dynamic effects to be considered 
have been discussed by Smith (14) and by Gore and Carroll 
(16) and will be only briefly outlined here. 

The pump is assumed to obey the affinity laws, and empiri- 
cal or theoretical data are used to obtain the effect of flow 
variations. The turbine nozzles are assumed choked, and 
only one turbine stage is represented by the equations of 
Table 2. If the turbine has more than one stage, one may 
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_ either use several sets of equations or obtain an empirical 
map to show the effects of speed, temperature and flow on 
performance. The basic low frequency feed system time con- 
stant is associated with the rotary moment of inertia of tur- 
bopump moving parts. Other high frequency effects arise 
from the fluid line inertias and the gas compressibility effects. 
A feed system analog mechanization diagram is shown in 
Fig. 7. 

Since these phenomena are fluid mechanical] in nature, it is 
clear that the overall feed system response may be made 
sufficiently fast as to be dynamically decoupled from the reac- 
tor thermal response. An example of the application of con- 
~ trol synthesis techniques to improving feed system response is 

shown in the root locus construction (17) of Fig. 4. For 
generality, the response is nondimensionalized in terms of the 
turbopump break frequency. The controller possibilities in- 
clude: 

1 Proportional, in which case the steady-state error will 
be large for gain values low enough to avoid instability due to 
high frequency poles. 

2 Integral, in which case the steady-state error is zero, but 
little improvement in dynamic performance is possible with- 
out oscillatory response due to high frequency poles. 

3 Integral with lead compensation. Here, zero steady- 
state error is obtained, and the effect of anticipation is a con- 
siderable improvement in the dynamic response. 


Interactions 


The system block diagram of Fig. 2 illustrates some of the 
basic subsystem interaction problems. The basic feature of 
system response is that the reactor power level is not uniquely 
determined by the propellant flow. In a chemical rocket, suit- 
able orificing of propellant supply lines to obtain proper pro- 


pellant flow rates automatically fixes thrust. In a nuclear J 7 
rocket engine, power perturbations effect both reactor and 
feed system operation due to temperature changes thereby 
produced. Increased temperatures will result in a rise of tur- 
bine output horsepower, if not controlled. This rise will even- 
tually be felt back at the reactor due to increase pump flow 
rate, with a resulting decrease in reactor temperature. 
Hence, the feedback loop of temperature through the feed 
system should be stable if the temperature coefficient. pre- 


dominates. However, if the positive reactivity effect due to dis 
propellant density changes is large, an increase of reactor U? 
power may feed through the turbopump to produce an in- Fo 


creased amount of propellant in the reactor, which in turn 
raises power farther. Thus, if the density coefficient is large, | 
not only does the system become much more sensitive to dis- pel 
turbances, but also the possibility of overall system instability J '! 
exists. The need for designing a nuclear rocket reactor to be J Sid 
both insensitive to density changes and to minimize the tie: 
amount of propellant in the core is thus apparent. 


Control Loops 

In view of the preceding discussion, it is clear that two dis- 
tinct control loops must be considered. The intrinsic speed 
with which the reactor power and the propellant flow can 
respond requires that these two parameters be automatically 
controlled at least during the dangerous period of engine 
transition. Temperature changes take place relatively slowly 
with respect to the speed of response of these two loops. 
Hence it is possible to consider pre-setting the desired maxi- 
mum values of power and flow before each run. The main- 
stage tolerance situation will then dictate whether a rather 
slow and simple automatic adjustment of temperatures is 
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needed during flight, or if a short preflight calibration test will 
suffice. 


Conclusions 


Considerations of nuclear engine applications and opera- | 


tional requirements indicate the need for advanced control 
methods if satisfactory and safe engine performance is to be 
obtained. To achieve the high degree of component reli- 
ability required to permit a realistic engine development pro- 
grain to succeed, extensive controls and other component test- 
ing under simulated system operation is required. This paper 
has presented a mathematical model which can be used to 
generate a simple, but useful, electronic simulation of system 
dynamic characteristics. 


Appendix 
Example of Simulation and Servo Analysis 
Techniques in Studying System Response 


Simulation of Kinetics Equations 

A useful mathematical technique which permits transforma- 
tion of equations into scaled form suitable for electronic simu- 
lation and simultaneously allows application of linearized 
servo analysis tools to the system, is that of nondimensional, 
nonlinear perturbation. Consider the basic neutronic kinetic 
equations without precursor loss 

i=1 


q N= 
Ci = (BN/l) — 


Let N = No (1 + n/100) and C; = Cio (1 + C,/100), and re- 
place p by p/100 where the new reactivity is in cents (or volts 
on an analog computer). The nominal steady-state values are 
related by 


6 


(p — 1) (BN/I) + > AC 


Cio = 
Then, the perturbed equations are _ 
1 pn 
(6+ B+ ) --* 


The values of the delayed neutron parameters have been 
discussed by Keepin (13) and are summarized in Table 3 for 
235, 


Form of Open Loop Neutronics Transfer Function _ 


Elimination of the nonlinear term in the previous equation 
permits a linearized transfer function to be formed by simul- 
taneous solution of the seven equations. For simplicity, con- 
sider only one delayed group with suitable averaged proper- 
ties. Then 


Table 3 Delayed neutron parameters (10) 
Delayed Half life Relative yield Decay constant 
group SEC B:/B ht = 0.693/11/2 
1 54.3 0.036 0.0128 
2.18 0.409 0.318 
5 0.135 1.51 
6 0.13 0.018 5.33 
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The transfer function is then 


1+ S/d 


; ‘The response is basically integral, modified by a lead-lag 


term. With six delayed groups, there would be six lead-lag 
terms. 

In order to evaluate a suitable averaged decay constant 
X, one method is to make the integral gains equal (7); then 


1 


Ae, 
Example of Response 


As an example, consider a prompt lag //8 = 0.01 sec, and a 
negative temperature coefficient of —5 cents per per cent 
temperature change. The simplified, one-group neutron 
transfer function is then 


These values have been chosen for convenience (a 0.01-uf 
capacitor and a X5 amplifier input were available), however, 
they are in the range considered representative by Schultz. 
With the previous estimate of a 2-sec thermal lag time (break 
frequency 0.5 radian/sec) the open loop transfer function 
becomes 


(1 + S/0.0805) 
(1 + S/0.5) (1 + 8/100) 


The closed loop response becomes clear from root locus con- 
sideration. At low loop gains, the integrating pole moves to 
the first zero at 0.08 radian per sec, thus stabilizing the sys- 
tem and yielding a flat response below the thermal lag fre- 
quency. The heat transfer and prompt lag poles at 0.5 and 
100 radians per sec move together and eventually break 
away from the real axis at high gains with oscillatory response. 
Normally, this situation would only be obtained in water 
moderated reactors in which the negative coefficient of reactiv- 
ity due to density changes is very large. For a nuclear rocket 
reactor, this oscillatory condition might be obtained if de- 
layed neutron loss were sufficiently high. 

For our example the heat transfer pole simply moves to 
higher frequencies yielding the following approximate closed 
loop response to external reactivity changes , 7 


pr _ 0.4025 
p 


0 0.2 (1 + S/0.0805) 


p (1 + S/0.069) (1 + S/3.04) (1 + 8/97) | 


0.2 (1 + 8/0.5) (1 + $/0.0805) 
p (1 +8/0.069) (1 + S/3.04) (1 + S/97) 


These transfer functions are shown in the form of frequency 
response plots in Fig. 3. The typical lead-lag response of 
power to reactivity is evident from the figure. At frequencies 
below the thermal lag frequency, the response is flat, but when 
disturbances are inserted at frequencies above the thermal lag 
frequency, the gain goes up as the stabilizing effect of the ther- 
mal lag disappears. The gain difference between low and 
high frequency response is closely associated with the tem- 
perature coefficient, thus providing a convenient experi- 
mental method for measuring the latter quantity. 

The effect of nonlinearities on response is evident in the an- 
alog computer results of Fig. 5. These results were obtained 
using the mechanization diagram of Fig. 6 and the previous 
data. For simplicity, the section representing precursor loss 
was omitted. At low reactivity inputs, the response is linear 
and as expected; however, as large positive reactivities are 
inserted, very large power surges occur before the tempera- 
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(+ 8/100) 
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ture feedback can become effective. Although it appears that 
| = ae —__———we the existence of a negative temperature coefficient will even- 
i tually stabilize these power excursions, it is clear that they 
are undesirable from the point of view of other components. 
Hence, reactivity must be closely controlled to guard agaist 
any sudden, large changes. 
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Nomenclature 


‘ Ans = reactor heat transfer area, in.? 
7 2 Ar = nozzle throat area, in.? 
Ar = reactor flow area, in.? 
Ayg = hot gas valve effective area, in.? 
_ Ary = turbine nozzle throat area, in.? 
Ar = main valve effective area, in.? 
b = fuel element thickness, in. 
: a = characteristic length for diffusion, in. 
a C = turbine nozzle spouting velocity, ips 
c* = nozzle characteristic velocity, ips 
jwd = turbine characteristic velocity, ips 
i Ci = ith group precursor concentration 
© Cpr = specific heat of turbine working fluid, Btu/lb R 
CG, = specific heat of propellant, Btu/Ib R 
: Cu = specific heat of fuel element, Btu/lb R 
5 d, = reactor coolant density, lb/in.* 
d, = propellant density, Ib/in.* 
svstew dy = fuel element density, lb/in.* 
¥ | D = precursor diffusion coefficient, cm?/sec 
= E; = energy yield per fission, Btu/fission 
wh F = vehicle thrust, lb 7 
Fig. 4 Feed system bode plot with lead compensated integral f = friction factor 
control GR = gear ratio 
g = acceleration due to gravity, in./sec” 
a h = heat transfer coefficient, Btu/R in.? sec 
= effective pump inertia, lb-in.-sec? 
43% AveRAce FUEL ELEMENT TEMPERATURE | Is = specific impulse, sec 
Ir = effective turbine inertia, lb-in.-sec? 
a ae : Ky = coefficient in pump H-Q curve fit, sec?/in. 
= coefficient in pump H-Q curve fit, sec?/in.? 
Ks = coefficient in pump H-Q curve fit, sec?/in.5 
| Ky = coefficient in pump torque curve, sec.2/in.4 
| Ks = coefficient in pump torque curve, in.-sec.2> 2 
NEGATIVE Positive Ke coefficient in pump torque curve, sec*/in.? 
+2% 1 K;, Ks = coefficients in turbine efficiency curve fit 
| NEUTRON LEVEL Ku = fuel element thermal conductivity, Btu/in. R sec 
= main propellant line inertness, sec?/in.? 
| Lr =: turbine torque, in.-lb i 
—— L, = required pump output torque, in.-lb 
a ie L = length of lumped parameter reactor section, in. 4 
TWO CENTS STEP IN Pp, l = effective neutron lifetime, sec 
m = hydraulic radius, in. 
AVERAGE FUEL| ELEMENT TEMPERATURE Mp = mass of propellant in reactor, lb 
Myo = nominal propellant mass in reactor, lb 
M, = reactor weight flow rate, lb/sec 
| M, = pump mass flow rate, lb/see 
| | | Mr = turbine mass flow rate, lb/sec 
| POSITIVE Nz = density of fissionable atoms, atom/in.*’ 
N = neutron density, neutron/in.* 
= pump speed, radian/sec 
Nr = turbine speed, radian/sec 
P, = reactor inlet pressure, psia i 
P, = reactor outlet pressure, psia s: 
| | | | Pp = pump discharge pressure, psia © 
—] Ge Ps = pump suction pressure, psia 
Pr = turbine inlet pressure, psia ; 
ONE DOLLAR STEP IN p, ae = turbine exhaust pressure, psia raW 
Reactor power and average fuel element temperature Py = bootstrap sensing pressure, psia 
response to step inputs of reference reactivity Qs = fuel element surface heat flux, Btu/in.? sec 
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_ Fig. 6 Reactor neutronics analog mechanization diagram 
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= reactor heat generation rate, Btu/in.’ sec 
= turbine radius, in. 

= effective resistance of main propellant line, sec?/in.* 
= vehicle mass ratio 

= gasconstant, Btu/bR 
= Laplace variable, sec™! 
= gas temperature i 
= fuel element temperature, R 
= nominal fuel element temperature, R 
= propellant temperature at reactor inlet, R 
= propellant temperature at reactor outlet, R 
= fuel element wall temperature, R 
= turbine inlet temperature, R 
= average neutron velocity, in./sec 
= reactor flow volume, in.® 

= vehicle weight, lb 

= total delayed neutron yield 

= neutron yield of the ‘th precursor group 
= temperature coefficient of reactivity, dollar/R 

= propellant mass coefficient of reactivity, dollar/Ib 

= specific heat ratio 

= fission cross section, barns 

= fuel element temperature distribution lag, sec i 
= thrust chamber storage lag, sec 

= reference reactivity, dollar = 

= reactivity, dollar 
= ith precursor decay constant, sec! 


| | 
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Applicability of Euler’s Dynamical 
Equations to Rocket Motion’ 


JAMES W. ELLIS? and CHARLES W. McARTHUR? 


Florida State University, Tallahassee, Fla. 


Let S denote a system of particles, such as a rocket, 
where angular momentum is transferred from the system 
by mass particles leaving the system. Rosser and asso- 
ciates have shown that the moment of the exterior forces 
acting on S is equal to the time rate of change of the mo- 
ment of momentum of S plus the rate at which the par- 
ticles that are leaving S are transferring moment of mo- 
mentum out of S, with all angular momenta being taken 
about the center of mass of S. From this principle, a set 
of equations for the rotational motion of a rocket with 
varying moments of inertia and center of mass is derived. 
This set of equations is a generalization of the standard 
Euler dynamical equations. The new terms which appear 
in them are seen to be negligible only under certain condi- 
tions. The set of equations reduces to the standard Euler 
dynamical equations when the mass and center of mass 
of the rocket are constant. 


NUMBER of mathematical studies of the rotational mo- 

tion of rockets are based on Euler’s dynamical equations. 
These equations have been derived for rigid bodies with con- 
stant moments of inertia about principal axes. Inasmuch 
as the moments of inertia of a rocket vary as its fuel burns, 
the question arises as to what extent the use of the Euler dy- 
namical equations is justified in studies of rocket motion. 
Jarmolow (1)4 claims to give ‘a heuristic proof of the appli- 
cability of these equations.” Using a principle of mechanics 
stated and proved by Rosser (2) and associates we derive, 
under assumptions given below, equations which are a general- 
ization of the standard Euler dynamical equations. Condi- 
tions under which the new terms which appear in them may 
be neglected are discussed in the conclusion. 


Assumptions and Notation 


We will discuss the equations of motion of a rocket which is 
losing mass through k openings in its base. The following 
simplifying assumptions will be made: 

1 The ejected particles are not given any angular ve- 
locity relative to the rocket. 

2 The system of openings through which the particles are 
ejected lies in a plane perpendicular to the longitudinal axis, 
and has 2n-gonal rotational symmetry about that axis for 
some n > 2. Symmetrically placed openings will be assumed 
to lose mass at equal rates. 
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3 The center of mass will be considered to always lie on 
the longitudinal axis of the rocket. 

4 All the mass leaving through an opening will be con- 
sidered as leaving through its center. 

The unit vectors %1, %2, Xs, will denote a right-handed 
orthogonal coordinate system, fixed in the rocket, with %, along 
the longitudinal axis and directed forward, and %2, ¥; in the 
plane perpendicular to this axis through the center of mass. 
In addition, the following symbols will be used: 


B = —Bzx, = the vector from the center of mass to C_ 
C = the center of symmetry of the system of openings 
in the base. (C is the intersection of the base 
with the longitudinal axis) 
h = angular momentum of the rocket about its center 
of mass 
[A]boay = time derivative of h with respect to the system 
Bo, 
[h]ipace = time derivative of h with respect to an inertial sys- 
tem 
I; = magnitude of moment of inertia about the #; axis, 
t = 1, 2, 3 (note that J, = J;) 
ki = radius of gyration of rocket about ; axis, 7 = 
2: 3 (ke ks) 
M = mass of rocket and contents at time ¢ 
m = mass of particles which have left the rocket by 
time ¢ 
mi = mass of those particles lost through the 7th opening 
at the instant ¢ 
k k 
= (> mee) [( p> 
i=1 i=1 
7; = vector from C to the center of the 7th opening 
3 
e = )°7,%; = sum of the exterior torques (about the 
i=1 
center of mass) acting on the rocket 
3 
w = wiki = angular velocity of the rocket, or of the 
$=] 
moving axes, with respect to an inertial system 
Z = the rate at which angular momentum is being trans- 
ferred out of the system by loss of particles 


Generalized Euler Equations 3 


Rosser (2) and associates prove that i . 
4 ots. 
T= 
Since %, 2, 3 are principal axes of inertia, it follows that 
From our assumptions it follows that 
= mq m( B2 + (WeXe W3X3) [3] 
It is well-known that 


[h |space = Jioay + WX h [4] 
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Also,wehave 


1; = = + M(kj*) 


j=1,2,3 since M =—m. 
Using Equations [2, 3, 4 and 5], Equation [1] yields the 
following system of coordinate equations 


Tw, + U3 — + mq? — + = Ti 
+ (I, — Is)wsw, + + $9? — + M(k2*)' we =T2 
+ (U2 — + m(B? + — + M(k2?)'ws = 73 


Rosser (2) has derived these equations in the special case 
Ww, = w; = 0 and with the openings equidistant from the 
center of mass. 


Conclusions 


Equations [6] differ from the standard Euler equations 
for nonvarying inertia 


Tw, + (Us — = 


by the addition of the terms m(q? — ki?)w, + M(k,?)’w,, 
m (B? + 49? — ke?)we + and m(B? + gq? — + 
M(ke*)"ws, respectively. These are the terms disregarded by 
Jarmolow (1). We note that making these terms negligible 
becomes a design problem. If the radii of gyration change 
slowly then the terms M(k;?)'w;, 7 = 1, 2,3 may be neglected. 
If, in addition, the rocket is so constructed that q? is near k,?, 
and B? + 49? is near k,?, these terms may be virtually neg- 
lected and the standard Euler Equations [7] used in place 
of the more complex set [6]. 

In the usual case, in which the direction of w is close to that 
of the longitudinal axis, only the first of these requirements 
(that q? be near k,?) is really relevant, due to the smallness 
of we and ws. 


1 Jarmolow, K., ‘‘Dynamics of a Spinning Rocket with Varying Inertia 
and Applied Moment,” J. Appl. Phys. vol. 28, pt. 1, 1857, pp. 308-313. 

2 Rosser, J. B., Newton, R. R. and Gross, G. L., ‘‘Mathematical Theory 
of Rocket Flight,’”” McGraw-Hill Book Co., Inc., New York, 1947, p. 15 
and Appendix 2. 
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Some Extensions of the Hohmann 
Transfer Maneuver 


T. N. EDELBAUM! 
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N 1925, Hohmann demonstrated (1)? that tangential im- 
pulses at the apsides of a cotangential transfer orbit 
minimize the fuel consumption for transfer between coplanar 
circular orbits. It has since been recognized that such 
maneuvers also minimize fuel consumption for transfer from 
a circular orbit to coplanar orbits which satisfy any desired 
functional relationship between energy and angular momen- 
tum and which have any desired orientation. Examples of 
these functional relationships are fixed energy (or major axis 
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or period), fixed angular momentum (or latus rectum), fixed 
apogee or perigee radii, and fixed eccentricity and major axis. 

In determining what combinations of tangential apsidal im- 
pulses are optimum for any given problem, it has been cus- 
tomary to consider only one- or two-impulse maneuvers. It 
is the purpose of this note to point out some examples of two- 
and three-impulse maneuvers which are more economical of 
fuel (or characteristic velocity) than the corresponding one- 
and two-impulse maneuvers. 

The first example is a comparison of the classic 
Hohmann transfer and a three-impulse transfer between 
circular orbits. For orbits which are more than 12 times larger 
(or smaller) than the original orbit, fuel can be saved hy 
transferring beyond the final orbit, applying a retrothrust at 
the apogee of this ellipse and entering the final orbit at the 


ov 
Yo | 


20 50 100 200 500 1000 


Fig. 1 


Characteristic velocity requirements for transfer between 
distant circular orbits 
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Fig. 2 Characteristic velocity requirements for a solar probe 
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ments for hyberbolic escape 


perigee of a second transfer ellipse. The further away the 
apogee of the transfer ellipses, the greater is the fuel saving. 
The absolute minimum fuel consumption necessitates escaping 
to infinity along a parabola, using an infinitesimal impulse at 
infinity to cause the vehicle to return along a new parabola 
which is tangent to the final orbit, and entering the final orbit 
at the perigee of this parabola. A sketch of this maneuver 
and its characteristic velocity savings are shown in Fig. 1. 
As the maximum possible savings in characteristic velocity 
are only about 8 per cent, the resultant fuel saving probably 
does not justify the added time and complexity of such a 
maneuver. However the next example represents a case 


where the fuel savings are so large that they more than justify 
the additional time and complexity. 

This example represents a solar probe in which it is de- 
sired to come as close as possible to the sun to sample its 
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atmosphere. It is possible to approach the sun to within one 
tenth to one hundredth of an astronomical unit (A.U.) before 
the probe melts. However this is a formidable propulsion 
problem. The conventional method of impacting the sun is 
to apply a large enough retrothrust to completely stop the 
vehicle in its motion around the sun. This requires a charac- 
teristic velocity of about 100,000 fps, more than twice as 
much as is required to escape the solar system. This propul- 
sion requirement can be greatly reduced by first applying a 
forward thrust to transfer beyond the Earth’s orbit and 
stopping the vehicle at the aphelion of the resulting ellipse. 
Fig. 2 shows this path, and the characteristic velocities re- 
quired to transfer from the Earth’s orbit to within various 
distances of the sun as a function of the aphelion distance of 
the transfer ellipses. The two-impulse maneuver is optimum 
for approaching to less than 0.2 A.U. For closer approaches, 
the characteristic velocity is a monotonically decreasing func- 
tion of the maximum excursion from the sun. Even going out 
only to the orbit of Mars produces appreciable savings in 
characteristic velocity. 

The third, and most complicated, example is concerned with 
maximizing the energy of hyperbolic orbits so as to maximize 
velocity at infinity. Lawden (2) has analyzed a two-impulse 
maneuver for this case. However a three-impulse maneuver 
can provide appreciable savings over both the one- and two- 
impulse cases. The maneuver consists of first going outside 
the original orbit with a forward thrust, passing inside the 
original orbit with a retrothrust at apogee, and finally escap- 
ing with a forward thrust at the perigee of the second trans- 
fer ellipse. In general, this maneuver will be limited by a 
minimum attainable orbital radius. Figs. 3 and 4 represent 
two typical cases, minimum radii of two tenths and eight 
tenths times the original orbital radius. In each case the 
velocity at infinity is plotted vs. the characteristic velocity, 
for various maximum excursions. The figures show that it is 
always possible to save fuel with this maneuver by going out 
far enough initially. The magnitude of the fuel saving de- 
pends upon how closely the origin can be approached. The 
theoretical optimum represents escaping to infinity along a 
parabola, using an infinitesimal impulse to fall back to the 
origin, and using another infinitesimal impulse at the origin 
to obtain any desired velocity at infinity. By closely ap- 
proaching this optimum it is possible to obtain very large 
velocities at infinity with characteristic velocities only slightly 
greater than that needed for escape. 
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Comparison of Minimum Energy Paths 
and Apogee Designation Paths 


WAYNE TEMPELMAN! 
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A comparison of minimum energy paths and apogee 
designation paths reveals that for distances greater than 
three times the injection distance there is little difference. 
Distances less than twice the injection distance reveal 
energy paths, 


significant differences. For minimum 


corresponding to each velocity there is a reachable en- 
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velope in space which constitutes an ellipse. 


Also, con- 


stant injection angles are found to be represented in space 
by straight lines, equally spaced, originating from the 
initial point. 


HE ANALYSIS is based on the usual simplified two-body 

assumptions. The formulas that are not derived are easily 
derived from a number of sources.? All distances are dimen- 
sionless, with the initial distance being unity. 


Minimum Energy Paths 

The locus of points representing the far focus is obtained by 
noting that the apogee angle is between the initial line and 
the major axis, and that the linear eccentricity extends to the 
far focus (see Fig. 2) 


tan 0, - 
y 


2 Singer, S. F. and Wentworth, R. C., ‘““A Method for Calculat- 
ing Impact Points of Ballistic Rockets: Convenient Representa- 
tions,’”’ JET PROPULSION, vol. 28, 1958, pp. 684-687. 


: 


This is a circle with origin at R = 1, 6 = 0 deg, and a radius 
of Vo2/(2 Vo"). 

Considering that for an ellipse the distance from a focus ‘o 
a point on the ellipse plus the distance from the point to the 
other focus equals the major axis, the maximizing of R for a 


give 6 occurs when gh = fh. Considering Fig. 2 
fh =9F -fl 


V,2 
R = (R? — 2R cos 6 + 1) 2 — Ve 


4V02/(4 
~ 1 — (2 — Vo2)/(2 + Vo2) 


- - - -APOGEE DESIGNATION PATHS 
MINIMUM ENERGY PATHS 


} 
Fig. 3 Reachable envelopes 
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The envelope of reachable points in R, @ space for a given 
velocity is therefore an ellipse with c = 3, putting the origin 
of the ellipse at R = 4, 0 = 0 deg (see Fig. 3). The axes are 


2 + Vo? hy, 


a= 


Considering the equation of line fT, the velocity components 
can be found using Equation [1], which lead to the injection 
angle 


_ Reos@ — 1 + (R? — 2R cos + 1)” 
[RU + cos 0)(R? — 2R cos + 1)'2]'/2 
B= — cos 6))'/2 
2R cos 6 + 1)'/4 
R cos @ — 1 + (R? — 2R cos 6 + 1)'7 
Yo = 
R sin 6 R 


Evuations [4] result in a plot of straight lines in R, @ space for 
a: onstant yo, originating from R = 1, 6 = 0 deg (see Fig. 4). 
Tie proof involves proving that the relation between the in- 
jection angle and the angle w (see Fig. 2). is given by 


w = 27) — 90 deg - [5] 
The relation w = f(R, 6) is found from Fig. 2 C 
w = sin [(R cos 8 — 1)/(R? — 2R cos 6 + 1)'] (6) 10° 


- “APOGEE DESIGNATION PATHS 
MINIMUM EXESGY PATHS 


Using trigonometry 
sin (270 — 90 deg) = —(1 — tan? yo)/(1 + tan? Yo) [7] 
By inserting Equations [4] into Equation [7], [6] is seen to be | 
true, proving Equation [5]. 
The formulas given in Equations [2 and 4] can be applied | 
for R = 1 (the solution for this special case can be obtained 


Fig. 4 Injection angles 


directly by maximizing the apogee angle) with the results — 
4 — 4V,? — Vo! 160 
= | —1+ (21 
1 + cos 6 
= Vo 2 ve) 80 
(2 — Vo2)'/# 40 
Yo = tan“"(1 — 
Fig. 5 Minimum energy paths for R = 


The general conic section formula can be written for the — 


coordinates of the apogee point 
' ee relationships for the semilatus rectum and the semimajor 


a(l — 6) axis 
2 
The eccentricity and semimajor axis are derived from Equa- .4 e™ ® a = (2 — B? — 1/a)7 


tions [8] 
R(R — cos R-1 \3 
R-1 R(R — cos 8) \2R-1-cosOL \R—cos8 
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Table 1 


Zero deg 


Special conditions for ¢ 
90 deg (MEP) 


= 0, 90 and 180 deg 
90 deg (ADP) 180 deg 


+ 1): 
R \ 


— 1) 
R 


AR — | 


= 
a. 


— 
R 


(R? + 1)'2 -1 


(2k? — 2k + aR 
R rs é + i) 


+ 
Ri/ ( 2k ) 
0 
(R? + ARI 
R? + 1)'/2 — 
e 
2 + R — + 2k 1 


R(R — cos 


R(R — cos 8) [1 -( R-1 \ 
R — cos 8 


R(R — cos 


2R — 1 — cos 8 


6) 


The injection angle is 
2R cos 6 -1 R — cos 6 R(R bie 


2R — cos 6 — 1 
cos 


—1 
R(R — cos 6) 


2R —cos# R — cos 


The reachable envelope for a given velocity is found using 


Equations [8 and 9] = 


It is seen from Figs. 3 and 4 that the solutions are identical 
for 6 = 0 and 180 deg. The results for these two cases along 
with @ = 90 deg are summarized in Table 1 


R(2 — V,*) — 


cos 6 = 


Nomenclature 
a = semimajor axis . 
b = semiminor axis ~ 7 
c = linear eccentricity 
= eccentricity 
f = far focus 
F = focus at center of attracting body 
g,h = points (see Fig. 2) 
I = injection point ‘- 
= semilatus rectum 
R = distance 
v = velocity 
V = dimensionless velocity 
a = dimensionless radial velocity __ 
8 = dimensionless normal velocity 
6 = polar coordinate (measured from F/]) | 
@ = angle (see Fig. 2) - 7 
Subscripts 
A = apogee 
ADP = apogee designation paths 
c = circular 
MEP = minimum energy paths 
0 = initial condition 
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New Ion Source for Propulsion of 


Space Vehicles 


BANTA, 
R. F. KING 


J. W. JOHNSON and 


Cc. D. MOAK,' H. E., 


Oak Ridge National Laboratory, Oak Ridge, Tenn. 


OME recent developmental work on high intensity ion 
sources for use in accelerators indicates that certain 
modifications of the duo-plasmatron source, originally de- 
scribed by von Ardenne in 1948 (1),? may prove to be well- 
adapted for low thrust, high specific impulse ion propulsion of 
satellites and space vehicles. 

Fig. 1 shows the configuration of a small duo-plasmatron 
now being used at Oak Ridge National Laboratory (2). 
Originally built for 10-milliamp proton service, this source has 
delivered up to the 50-milliamp capability of the power supply 
used in the tests. Throughout this current range, the source 
can deliver up to 100 amp of protons per cm? of exit aperture 
area, in agreement with figures taken by von Ardenne for a 


proton current of 1 amp. Power consumption is only 400 w, 
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and gas efficiency is very nearly unity. For one or more 
amperes of protons at accelerator voltages of 10 kv or more, 
beam power would constitute 96 per cent of the total power 
consumption of the propulsion system. Total source weight 
is approximately 5 lb. 

It has been shown (3) that optimization of the propulsion 
system design calls for a choice of exhaust velocities which 
makes most efficient use of the power supply and propellant 
loads. In some cases the design calls for programmed exhaust 
velocity (4). In most instances the use of heavy ions with 
accelerator voltage of some tens of kilovolts is dictated. The 
use of light ions would dictate very low accelerator voltage 
and such short gap spacings as to make the accelerator rather 
small and fragile. Experience with He and N; as feed gases 
indicate that no great difficulty should be experienced with 
heavy-ion operation of the duo-plasmatron. The problems to 


be encountered with space charge neutralization and beam 
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MAGNET COIL 2000 amp turns — 
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—+—0140in. 
-95° TO 100° 
30°" 
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Fig. 1 Duo-plasmatron ion source 7 


expansion should be no different with this than with any other 
source of ions for this type of service. 

Clearly, the weight of the accelerator can be small com- 
pared to that of the power supply and propellant loads so that 
power supply weight now becomes the most important item 
to be minimized. Even with power sources presently availa- 
ble it appears that several applications of ion propulsion are 
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ADOLPH B. AMSTER! and JOSEPH B. LEVY? 


U. S. Naval Ordnance Laboratory, Silver Spring, Md. 


The presence of bubbles is known to sensitize liquid pro- 
pellants. The magnitude of the effect is determined by 
the composition of the gas in the bubble; this in turn de- 
pends upon the solubility of various gases in the liquid. 
This paper presents data for the solubility, at ambient 
temperature, of air, No and O2 in n-propyl nitrate. Ab- 
normal solubility of oxygen and the consequences thereof 


are noted. 
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Fig. 2 Solubility of air in in n-propyl hee. 


n-Propyl Nitrate 


HE SOLUBILITY of different gases in n-propyl nitrate 
(NPN) could conceivably be an important factor in the 
determination of the safest ways to use this propellant. If, 
for example, in its use a pressurization by some gas were to be 
followed by depressurization, the evolution of gas in the form 
of bubbles or foam would constitute a potential hazard (1, 2).* 
The degree of such gas evolution depends in part on the gas 
solubility; such data seemed to be lacking at present. 


Procedure 


An apparatus has been completed with which a study of the 
solubility in NPN has been made of various gases at pressures 
below 1 atm and at ambient temperature. The equipment 
has been adapted from that described elsewhere (3). 

The apparatus is shown in Fig. 1. In operation a weighed 
quantity of NPN is put into the flask F. Dissolved gases ::re 
removed by alternately freezing and thawing the NPN, and 
pumping on the frozen material until a constant residual prvs- 
sure is noted. With stopcocks 2, 3 and 4 open, gas is ad- 
mitted to 1 until the desired pressure P is reached; 2 anc 3 
are closed, and F is shaken vigorously enough so that the 
liquid within is turbulent. Work elsewhere (4) has shown that 
a negligible quantity of gas dissolves prior to shaking. The 
change in level of the mercury in capillary C indicates the 
change in pressure within that section confined by the flask 
and flexible spiral up to stopcock 2 and mercury level L (in 
the widest part of the bulb). The volume of this section V 
is known,‘ and it is a relatively simple matter to calculate the 
solubility of the gas. The final expression for doing this is 


= V’AP/WR 


3 Numbers in parentheses indicate References at end of paper. 
4 See (3) for the method of obtaining V. 
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where 

S, = moles gas absorbed per gm. NPN at equilibrium pres- 
sure P, 

V’ = total volume V minus volume of NPN 

P,. = total pressure at equilibrium minus NPN vapor pres- 
sure 

AP = pressure change during equilibration be « 

W = weight of NPN used a - 

R = gas constant : 

7 = absolute temperature 


In our calculations we have used values for the density 
and vapor pressure of NPN given by Thomson (5). 


Results 


The solubilities of air, oxygen and nitrogen in NPN have 
been studied at room temperature and over a pressure range 
of spproximately 200 mm to 1 atm. 
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Fig. 4 Solubility of nitrogen in n-propyl nitrate 
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Fig. 5 Comparison of solubilities of air, O. and N» in n-propyl 
nitrate 


Technical grades of oxygen and of water-pumped nitrogen 
were used. The gases, assaying better than 99.5 per cent 
purity, were passed through a drying tube before being ad- 
mitted to the system. The air was similarly dried. 

Results are plotted in Figs. 2, 3 and 4 as S, vs. Pe. 

A least-squares treatment of the data has given the fol- 
lowing equations 

air S, = 5.787 X 10-°*P, + 4.265 X 10-4P,? 

nitrogen S, = 7.371 X 10-°P. — 7.277 X 1073P,? 
Sp, = 7.515 X Pe + 6.012 X 107}? P,? 


~“I 


oxygen 


The resultant curves are superimposed on the data in Figs. 
2, 3 and 4 and are collected in Fig. 5 where they may be com- 
pared with each other. In addition, there is a curve for a 
hypothetical mixture consisting of 80 per cent nitrogen and 20 
per cent oxygen. This is predicted from the solubility of each 
gas at its own partial pressure. 


Significance of Results 


It is seen that the oxygen solubility is significantly greater 
than that of nitrogen. Furthermore, the solubility of air is 
greater than that predicted on the basis of the oxygen and 
nitrogen data. The higher solubility of oxygen indicates an 
unusual, but not unexpected, interaction between it and 
NPN. A possible interpretation of the abnormal air solubility 
might be that the dissolved air is enriched in oxygen, at the 
expense of the nitrogen. This is further supported by the 
apparently linear behavior of the nitrogen which indicates, for 
this gas, no interaction with the solvent. 

If it is true, as seems to be the case, that nonideal absorp- 
tion occurs, then it should be borne in mind that the data 
presented here for air are dependent upon the final gas to 
liquid ratio as well as the equilibrium pressure. This restric- 
tion does not apply to nitrogen and oxygen. 

A consequence of this phenomenon is that bubbles formed 
by the release of pressure during the pumping of NPN are apt 
to be significantly enriched in oxygen. This would, of course, 
present a distinct hazard. This same unexpected solubility 
phenomenon is observed with water. At atmospheric pres- 
sure air dissolved in water is about 34 per cent oxygen (6). 
It must be strongly recommended in the absence of contra- 
dictory information, that an inert gas be used when pumping 


the liquid NPN. 
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Magnetogasdynamic Flow Regimes 


The prediction of plasma flow subject to a magnetic field 
is of considerable interest to the engineer designing mag- 
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PHREE dimensionless parameters are commonly used in 


magnetogasdynamics. These are the magnetic Reynolds 
number (1),3 the Hartmann number (1) and the Lundquist 
number (2), and are defined respectively as follows 


Rem = VLy,,~” Oy {1] 


Lu = [3] 


In addition to these dimensionless parameters at least two 
others should be considered because they conveniently define 
the flow conditions. Both of them are based on the Alfvén 
wave propagation velocity V4. 

In our hypothetical model we assume that the plasma is 
described by the thermal] equation of state p = pRT. It is 
also predicated that the electrical conductivity of the plasma 
o — ©, and that the permeability » = 1. Finally, we con- 
sider steady, one-dimensional, continuum flow. Thus, if the 
plasma flow is in the z-direction and a homogeneous magnetic 
field Hy is applied in the z-direction, then an electric current 
flowing in the y-direction is produced. This in turn produces 
a variable magnetic field H’ in the direction of flow. Alfvén 
(3) has shown that the magnetic disturbances propagate with 
the Alfvén velocity V, defined as 


Hy Bo 

It has been shown by deHoffmann and Teller (4) and by 


Dungey (5) that the speed a,, with which small disturbances 
in magnetogasdynamic systems propagate is given by 


ay =a + Va [5] 


We therefore suggest in this note that a magnetic Mach 
number ,, be defined as 


My = V/(a + Va) [6] 
An associated magnetic Mach angle a, may then be defined 
as 

sin~! ay = 1/My [7] 


It follows that for the same plasma flow with and without a 
magnetic field, the Mach numbers and Mach angles are 


My <M 
Simple laboratory experiments can be conceived, to verify 
Equations [8 and 9]. This would also make possible the con- 
venient measurement of the speed of small disturbances in 
magnetogasdynamics. 
For studies in astrophysics, Cowling (1) has used the S-ratio 
which we shall call the Cowling number C defined as 


C=+~-— [10] 


In Equation [10], the numerator represents the energy asso- 
ciated with the propagation of a magnetic disturbance, and 
the denominator represents the kinetic energy of the fluid. 
In effect, C is analogous to 1/M? which is used extensively by 
aeronautical engineers. It seems more consistent, however, 
to write the ratio 1/VC. We shall do so, and in honor of the 
eminent engineer Theodore von K4rm4n, we shall call this 


3’ Numbers in parentheses indicate References at end of paper. 
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[11] 


It is evident from Equation [11] that the Karman number is 
an indication of the magnetic field dissipation. The applica- 
tions of this in the analysis of magneto-aerodynamic devi:es 
are obvious. Thus, if K < 1, the disturbances in the field :.re 
propagated as magnetohydrodynamic waves without prodiic- 
ing any appreciable changes in the field. On the other hand, 
when K > 1, the magnetic disturbances produce distortions at 
such a rapid rate that they cannot be dissipated as magne ‘o- 
gasdynamic waves. Therefore, we conclude that magne‘o- 
gasdynamic problems are of engineering interest when A is 
small. 
Logarithmic differentiation of Equation [11] gives i 


It follows that particularly at high velocities the Karn in 
number is a measure of the density and thus of the co n- 
pressibility. 

The logarithmic differential of the Mach number is (see | or 


example (6 

example (6)) 
[13 
M 2T 


dM dV 


If we combine Equations [12 and 13] and perform the 
algebraic manipulations, we arrive at an equation relating thie 
Mach number and the K4rm4n number 


2dM dK dp 


2M+(y-—1)M? 2p 


Equation [14 | will be helpful in predicting the density varia- 
tion in magnetogasdynamic flow. Thus we have compression 


re 


M 
[2+ (y — 


2dM 
OM DM? K 


or when 


M 


In magnetogasdynamics, the continuity equation remains 
unchanged, but the momentum and the energy equations 
must be modified. (See for example (7).) Recalling the 
definition of the Alfvén wave propagation velocity, the energy 
equation for adiabatic flow becomes 


a? + V4? = constant [19] 


2 
ve + 


4 At the completion time of this note, it came to the authors’ 
attention that in a lecture at Amsterdam (to be published in the 
Proceedings of the IX Annual Congress of the International 
Astronautical Federation; Lange and Springer), entitled ‘Mag- 
neto Fluid Dynamics in Relation to Space Flight Problems,’ Dr. 
von Kérmdén himself suggested the use of this parameter and 
designated it as the magnetic Mach number. 
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We can represent Equation [19] geometrically as in Fig. 1 in 
which several magnetogasdynamic regimes are recognized. 
Thus in regime A we have subsonic flow with no appreciable 
magnetic effects, whereas in regime B the flow is supersonic, 
though still without important magnetic effects. In regimes 
C and D the magnetogasdynamic phenomena are important. 
The flow in regime C is supersonic; in D it is subsonic. 

In classical gasdynamics, the equations resulting from the 
various combinations of the mass and momentum equations 
are called the Rayleigh equations. Analogous magnetic 
Rayleigh equations may be derived by combining the classical 
continuity equation with the magnetogasdynamic momentum 
equation. If we do this, and use in the algebraic manipula- 
tions the equation of state for a perfect gas as well as the defini- 
tious of the Mach and Kaérmdn numbers, we get the following 
magnetic Rayleigh equations 


+ + 1/2K,?) 
Pp 1 + + 1/2K,2) 


[20] 


T, + yM,%(1 + 1/2K,*)] 


[21] 


(Vmax, 0,0) V 


ic flow regimes 


Magnetogasdynam 


[22] 

For the case when K — ©, there is no magnetogasdynamic 

effect, and the equations reduce to their usual counterparts in 


classical gasdynamics. 


Nomenclature 


a = acoustic velocity 


aa = magnetic acoustic velocity 

B = magnetic field vector @4 » 

C = S, Cowling number _ 

H = magnetic intensity 

Ha = Hartmann number 

K = Ké&rmén number 

L = characteristic length 

Iu = Lundquist number 

M = Mach number 

My = magnetic Mach number 

Pp = pressure — 
Rem = magnetic Reynolds number =~ 
8 = entropy an 
T = temperature 

V = flow velocity 
V4 = Alfvén velocity ; 
a = Mach angle ag 
ay = magnetic Mach angle bd 

¥ = specific heat ratio 

n = fluid dynamic viscosity 

= permeability 
Ym = Magnetic viscosity 

p = density A 
o = electrical conductivity 
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Technical Comments 


iF Comment on “Generalized Trajectorie 
for Free-Falling Bodies of High Drag” 
and “Vertical Re-Entry Into Earth’s 
_ Atmosphere for Both Light and Heavy 
Bodies” 


te 


JERZY H. TEISSEYRE! 


Technical University in Wroclaw, Wroclaw, Poland 


N THE April? and May® 1958 issues of JET PROPULSION, 
there appeared investigations concerning a free-falling 
_ body coming from very great height through the Earth’s at- 
_ Earth’s acceleration due to gravity is constant and equal to 
9.81 m/sec? or 32.2 fps. 
_ | shall discuss here a rapid method of calculation where the 
Earth’s variable acceleration is taken into account and dif- 
ferent drag coefficients can be inserted in the calculus. More- 
over, the exponential density equation p, = poe” can be 
considered with variable coefficient k as is really the case. 
The atmospheric data assumed are given in Table 1 and 
also represented in Fig. 2. 

_ Considering Fig. 1 and assuming the following nomencla- 


m = mass 

g = acceleration of gravity on the Earth’s surface 

gx = g(R/R + x)? acceleration of gravity at the height x 
above the Earth’s surface 

v = velocity, m/sec 

p = airdensity, kGsec?/m‘ ge 

A = cross-sectional area of body, m? 


= 


Earth and the adopted designations 
1 Received April 27, 1959. a 
‘Chairman, Polish Astronautical Society, Wroclaw. 


* Turnacliff, R. D. and Hartnett, J. P., ‘““Generalized Trajec- 
tories for Free-Falling Bodies of High Drag,’’ Jer PRopuLsIoN, 
vol. 28, April 1958, pp. 263-266. 

Linnel, R. D., rtical Re-Entry Into the Earth’s Atmosphere 
for Both Light and Heavy Bodies,’’ Jer Propu.sion, vol. 28. 
May 1958, pp. 329-330. 
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Table 1 Atmospheric data of free-falling bodies 
Height, m X 10 Ratio density 
Po 


~ 
x 
x10 
83 X 1076 
x 107-8 
x 
x 10~° 
90 
13 X 107% 
70 X 107 


OR 
w 


* The unequal n 
meters. 


umbers result from converting feet into 


Ca = drag coefficient = drag/(1/2)pv2A 


we get the following equations of a vertically falling body 


R 2 I 
where D = (1/2) prCqAv? is the force of drag. Thus 
or 


( 
#=-g 

where constant 


150 200 250 300 Km 


50 100 
Fig. 2 Adopted characteristics of the atmosphere 


RS JouRNAL 


| 
0,887 
0.1625 
91.44 
195.3 
243.84 
296.5 
370.64 
t 
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| 
: 2 m 
<a 
13000 : 
;. 3 10° 3 ~ 
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8000 
10 7000 


ta 


[3] 


FAL 


Table 2 Example of calculation 

u = 0%, 

h, km z,km Az, km 2gR?/(r — z)? k, km 2E- /ehk = 2G-e* 2Ge*?-u du/dz m/sec? 

100 0 0 19.017 1/8 0.00000933 0 19.017 0 

95 5 5 19.04 1/7.9 0.00001327 1.261 17.78 95085 

90 a, 19.07 1/7.8 0.0000250 4.35 14.72 173985 

85 15 5 19.10 ae 0.0000380 9.415 9.685 247585 

84 16 1 19.112 We Gees 0.0000433 11.135 977 257270 
Putting also z = h — 2 (Fig. 1) we have: e*” = e~-&, conditions du/dz + 2gR?/r? = 0, from which du/dz can be 


where / is the height from gy the body rt to fall. At 
this height the velocity v = z = 0. Putting G = E-e~™ 
(which is variable, k being v. siiihe as well) we get 


and since 


dv dv dz 


dt dz dt 


and putting v? = u, we get finally 


+ 
dz (r 


This equation can be readily solved by step-by-step calcula- 
tion if we know the initial conditions. Now the initial condi- 


tions are: u = 0 for x = horz = 0. So we have for initial 
acceler. (mfsec’) 
100 10-8-6-4202 4 6 610 

90 

80 

70 

60 


02 03 04 05 08 G7 08 09! 13 


Fig. 3 Speed and acceleration of falling bodies (parachute and | 
missile) 


TIME OF 


300 sec. 


Fig. 4 Time of fall for parachute and missile 
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calculated, and for an assumed Az (say of 5 km) the value of 
Uaz = (du/dz) Azcan be obtained. Putting this ws,into Equa- 
tion [5] and assuming again a certain Az, we find another 
value of du/dz and also another ua, = u. For purposes of 
quick calculation, Table 2 is useful. F rom Table 2 one can 
also calculate the accelerations of the falling body and the 
time of fall. 

Figs. 3 and 4 show the results of the calculations for two 
different bodies: One of a very large drag and small mass, 
(a parachute), with 2E = po (C,A/m) = 2.5/m—1|, and the 
second one for a small drag and large mass (a missile), with 
2E = 3 X 10-%)m~|. The cases of medial values 2E = 0.003, 
0.3 and 1 have been also considered. 

The balance of energy of the falling body is given by the 
equation 


v,2dz (per unit of mass) 


[6] 


1 h—- 
energy = —2z)+ v2? + ¢ 


where g, = e*?-e~-™. The first term gives us the potential 
energy, the second term kinetic energy, and the third one 
represents the work necessary to overcome the drag of the 


body. 


The condition of constant energy (d energy = 0) gives 
us 
—gdz + vdv + Egy, dz = 0 7] 


from which we can get the increase of sveed per unit of tra- 


jectory. 


Comments on “Recent Advances in 


Cryogenic Engineering” 


Ss. L. BRAGG! 
F Rolls Royce Ltd., Derby, England 


WOULD like to comment on the article “Recent Advances 

in Cryogenic Engineering” by Robert B. Jacobs which 
appeared in the April 1959 issue of ARS Journat, and in 
particular to disagree with Dr. Jacobs’ statement on page 
249 that ‘Dynamic pressure measurement (in a cryogenic 
liquid) offers no difficulty if the transducer temperature is 
known.” 

It appears to be normal practice in cryogenic engineering 
to keep transducers at ambient temperature so that an 
accurate calibration may be maintained. Unless, therefore, 
a water-warmed diaphragm type of pickup (such as the 
Li-Liu) is used, it is necessary to connect the transducer 


to the tapping point by a pipe some 10-in. long. Under 
Received June 9, 1959. 
1 Chief Performance Engineer, Rocket Engine Division. 
Member ARS. 
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6 Fay, J., J. Chem. Phys., vol 20, 1952, p. 942. 
™ . y conditions, this Pee will be full of ape, oe ept in the 7 Chu, B. T., ‘Vibration in a Gaseous Column behind a Strong Detona- 
a now the re in the main liquid line denly rises, orthwestern Univ., Aug. 1955, p. 95. : 
preem 8 Manson, N., Comptes Rendus del’ Académie des Sciences, Paris, vol. 222, 
liquid will flow into the tapping line, and the vapor there will 1946, p. 46. 


be compressed. However, this compression will occur much = 
= more slowly than if it were a gas pressure that was being = ug: _ ae | 


measured. This is because the volume inflow under a given 
‘s = pressure difference between main line and tapping is inversely 


La proportional to the square root of the density. The response ql 


of the transducer, which records the pressure of the vapor, 
will therefore lag the pressure change of the liquid in the line. Simple Method of Estimating the Varia- 
Finally, there may even be an uncomfortable overshoot, due tion of Impact Points for Ballistic 


to the momentum of the liquid in the tapping being suddenly > Missil 
arrested by the buffer of residual vapor at the transducer. 
Thus, measurement of transient pressures in cryogenic SAMUEL W. FORDYCE! 


liquids does produce an extra difficulty in that the response 
_ of a transducer which is not actually mounted at the tapping Space Technology Laboratories, Inc., Los Angeles, Calif. 
point will be considerably slower than it would be in either 
a completely gaseous or a completely liquid system. al 
PREVIOUS article? derived the following expression for 
2 the free flight range of a ballistic missile 


S = 2Rz sin (Yz? sin 20/2.) 


Commenis on “Recent Advances in a 
where 
Gaseous Detonation Re = radius of the Earth 
G. RIBAUD! = dimensionless velocity = vz/v. 
= cutoff speed of the missile 
: Faculté des Sciences, Paris, France : Ve = orbital speed at the = altitude (Az from the 
Earth’s center) = V 
q APERS by Déring and Schén (1)? concerning mixtures of GMz = constant, 1.405 X 10" ft® sec? 
2 CH,-O, and C,N,-O, should be added to the bibliography 6 = angle between the velocity vector at cutoff and the 
given by Gross and Oppenheim in their recent article in this local horizontal 
Journal (2). Déring and Schén avoided the theoretical € = orbital eccentricity, given by 
errors in their calculations indicated by Gross and Oppenheim, [1 — cos? (2 — (2) 
and even though the authors of (1) used outdated thermo- % 
dynamic values (3), their results are noteworthy. To my Suppose a small increment of velocity Av is added linearly 
knowledge, these theoretical errors have up to now been to the velocity vector at cutoff. The revised cutoff speed is 
avoided only by Manson (4), Déring and Schon (8), and Eisen, ‘ 
Gross and Rivlin (5). vp’ = vp + Av af [3] 
their discussion of non-Chapman-Jouguet detonations, The dimensionless velocity becomes 
Gross and Oppenheim indicated that the theory permitting j qj 
prediction of the frequency of the spin, which is often visible Yx' = Yi + (Av/»,) (4] 


in photographs of the waves and which characterizes their 

instability, was developed by Fay (6) and Chu (7). Ac- 

tually, this theory had already been discussed by Manson 

_ (8). Fay cited Manson’s work and noted some differences : 2AvY_z cos? 1 — Yz?) 

the treatment of the phenomenon. However, Chu later [1 ] 
showed that these differences were minor ones. 


References 


and since Y,>>Av/»,, the revised orbital eccentricity is 
given approximately by 


[5] 


ey, 


The revised angular range is given approximately by 
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231. 
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Ramo-Wooldridge | 
Laboratories 


in Canoga Park 


...an environment dedicated to 
technological research and development 


academic environment necessary for creative work. The 
new Laboratories will be the West Coast headquarters 
of Thompson Ramo Wooldridge Inc. as well as house 
the Ramo-Wooldridge division of TRW. 

The Ramo-Wooldridge Laboratories are engaged in 
the broad fields of electronic systems technology, com- 
puters, and data processing. Outstanding opportunities 


The new Ramo-Wooldridge Laboratories in Canoga 
Park, California, will provide an excellent environment 
for scientists and engineers engaged in technological 
research and development. Because of the high degree 
of scientific and engineering effort involved in Ramo- 
Wooldridge programs, technically trained people are 
assigned a more dominant role in the management of 
the organization than is customary. exist for scientists and engineers. 
The ninety-acre landscaped site, with modern build- For specific information on current openings write 
ings grouped around a central mall, contributes to the to Mr. D. L. Pyke. 
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Hypersonic Flow Theory, by Wallace D. 
Hayes and Ronald F. Probstein, Aca- 
demic Press, New York, 1959, 464 + xiv 
pp. $11.50. 

Reviewed by Artuur A. Kovitz 
Northwestern University 


DE: This volume is the fifth in a series on ap- 
mathematics and mechanics, editor- 
in-chief N. Frenkiel, prepared under the 
auspices of the Applied Physics Labora- 
tory, The Johns Hopkins University. This 
= has had occasion to use the first 
er third volumes, namely “Gas Dynam- 
s,”’ by K. Oswatitsch, and ““Mathemat- 
ical Theory of C Yompressible Fluid Flow,” 
by R. von Mises. To those also familiar 
with the Oswatitsch and von Mises efforts, 
it should be quite clear that any companion 


ment. It is stated at the outset that this 
latest addition does so and excellently 
compliments its forerunners. 

An outline of its contents should suggest 
Je the broad coverage that has been achieved. 
Beginning with small disturbance and 
: Newtonian theory, constant density solu- 
A tions for various shapes are discussed. 
-_ The theory of thin shock layers followed 
by other methods for blunt body and 
locally supersonic flows receive individual 


chapters. The concluding three chapters 
: deal with viscous effects, including inter- 
actions, molecule and rarefied gas flows. 
___‘The introductory chapter is not one to 
_ be skimmed lightly. After a few initial 
_ remarks on the ‘definition’ of hypersonic 
flow, a section is devoted to the assump- 
tions underlying inviscid hypersonic flow 
theory. Five fundamental assumptions 
J are listed and discussed with respect to 
_ their use in particular theories. The re- 
_ maining chapters on inviscid flow continu- 
_ ally refer to these explicitly stated proposi- 
_ tions as a basis for development of vari- 
- ous analytical schemes. It is decidedly 
_ brought home to the reader that a clear 
grasp of the fundamentals is essential to 
_ his subsequent understanding of the re- 
maining work. 
_ The treatment in the remaining sections 
of the first chapter illustrates one of the 
 leitmotifs of the following chapters on 
inviscid hypersonic flow. Wherever pos- 
i sible, existing theories applicable to per- 
‘ feet gases are formulated in terms of a 
general ideal (nonviscous) fluid. The nor- 
~ mal shock is used to illustrate the behavior 
of the density ratio, and several interesting 
- calculations are made to aid the reader in 
- acquiring an understanding of its varia- 
tions. Oswatitsch’s result, that for a per- 
fect gas the flow pattern becomes indepen- 
_ dent of free stream Mach number as this 
‘parameter becomes very large, is shown to 
be valid in a deeper sense for a general 
fluid. This is stated as a “Mach number 
independence principle.’”’ The first chap- 
ter is concluded with a brief discussion of 
real-fluid effects. 


Small disturbance and — 


Ali Bulent Cambel, Northwestern University, Associate Editor 


sonic similitude are treated in the second 
chapter. A feature of this chapter is a con- 
cise, general description of the concept of 
similitude itself. Some early Russian 
work on small-disturbance theory is dis- 
cussed. 

Newtonian theory receives detailed, 
deep attention. It begins with a highly 
diverting excursion into the pertinent 
parts of ‘Principia Mathematica.’ The 
limitations of Newton’s mode] are clearly 
pinpointed. Busemann’s centrifugal cor- 
rection to the sine-squared pressure for- 
mula is pointed out as an essential modifi- 
cation of the Newtonian flow model. Opti- 
mum shapes and shock layer structure in 
Newtonian theory are given treatments, 
some parts of which have not previously 
been reported in the literature. 

A chapter on the constant density ap- 
proximation is devoted to solutions for 
several simple bodies. Included are the 
wedge, cone, circular cylinder and sphere. 

Thin shock layers are given a chapter of 
unique explanation. The senior autheor’s 
artesian-well and shower-bath analogies 
are discussed in connection with the gen- 
eral problem of finding the shock shape, 
given the body shape. Several approxi- 
mate methods of attack are outlined in- 
cluding successive approximations, con- 
stant-streamtube area, and inevitably, 
variable-streamtube area. 

The last two chapters dealing with the 
ideal (nonviscous) fluid present numerical 
methods for treating blunt body and 
locally supersonic flows. Relaxation tech- 
niques, integral methods and the inverse 
approach are detailed for finding the shape 
and location of the detached bow shock 
given the body or vice versa. The method 
of characteristics is reviewed for a general 
fluid. Special attention is given to flows 
with equilibrium or frozen composition. 

Two chapters are devoted to viscous 
flows and viscous interactions. The chap- 
ter on viscous flows is a comprehensive 
but highly condensed review of the intense 
efforts of the last five or six years to under- 
stand boundary layer flows with relaxa- 
tion, diffusion and chemical reaction. The 
boundary layer equations are quoted in 
complete generality and transformed under 
Mangler and Howarth-Dorodnitsyn trans- 
formations as generalized recently by Lees. 
Similar laminar boundary layer solutions 
are described, and the limitations in find- 
ing them are pointed out. 

The chapter on viscous interactions con- 
tains much of the second author’s valuable 
contributions to the field in addition to 
what appears to be an excellent review. 
Much of the work done at Princeton Uni- 
versity’s Gas Dynamics Laboratory in the 
helium hypersonic wind tunnel is given 
well-deserved attention. 

The final chapter is on free molecule and 
rarefied gas flows. It is suggested that the 
criterion for rarefied vs. continuum meth- 
ods must be based on local flow conditions. 
A revised classification of flow regimes is 


presented that attempts to delineate more 
clearly the various regions, from con- 
tinuum flow at high Reynolds number to 
free molecule flows at very low Reynolds 
number. Continuum solutions for rare |ied 
gas flows are discussed after which ‘ree 
molecule transfer theory is outlined. 

Two helpful organizational features of 
the book should be mentioned. First, he 
cited references are listed alphabetic: lly 
according to author, with titles, in ne 
section. Second, a symbol index is | ro- 
vided with the page number where ‘he 
symbol first appeared. 

It has already been remarked that one 
of the aims of the authors’ presentation 
has been generality and completenvss. 
Thus, the work takes on the character «f a 
treatise, written in a sophisticated «nd 
terse style. As such, it can present com- 
prehension difficulties to readers not ‘00 
familiar with hypersonic flow. This vol- 
ume is certainly not of a textbook nature. 
The subject is complex, and one mist 
bring to it a rather well-rounded back- 
ground in order to make good use of the 
presentation. 

This reviewer has enjoyed encountering 
a book written by two authors well-pre- 
pared for the task of putting into one 
place the results of fairly recent research 
activity. The authors themselves have 
contributed significantly to the field. 
They write with a deep understanding of 
the present state of the art, and the cer- 
tain knowledge that present research 
activity, their own included, will demand 
future books on hypersonic flow. How- 
ever, the nature of the present volume 
makes it an important contribution to the 
contemporary scene as well as a book that 
will remain vital for many future years. 


Some Problems in Chemical Kinetics and 
Reactivity, Vol. 2, by N. N. Semenov, 
translated from the Russian by Michael 
Boudart, Princeton University Press, 
Princeton, N. J., 1959, 330 pp. $4.50 

Reviewed by 8S. 8. PENNER 
California Institute of Technology 


In the present volume, the following 
topics are treated: Competition between 
chain reactions and reactions between 
saturated molecules (including explicit 
discussion of such diverse topics as forma- 
tion of NO from and decomposition 
of di-iodo-ethane and of alkyl bromides, 
cracking of hydrocarbons, etc.), thermal 
explosions, chain ignition (e.g., in the 
oxidation of phosphorus and of sulfur, 
ignition limits in the oxidation of phos- 
phine, silane, ete.), chain ignition in 
H.-O, mixtures, chain interactions incliid- 
ing cold flame propagation in lean oxygen- 
carbon disulfide mixtures, and chain reac- 
tions with degenerate branching (e-g., 
oxidation of CH, H2S, liquid-phase oxila- 
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GAS PURIFICATION 
GAS SEPARATION 
AIR SEPARATION 
_-—-—- CONVERSION OF GASEOUS 0; and 
N2 TO LIQUID PRODUCT 


and fe METAL TREATMENT 


LIQUID AIR PRODUCTION 
LIQUID NITROGEN PRODUCTION 
PRESERVATION OF BIOLOGICALS | 
COLD MACHINING 
All of the applications listed—and more—stem from the one unit— LOW TEMPERATURE PULVERIZATION 
THE NORELCO CRYOGENERATOR. Available in four different capacities, a 7 
this remarkable unit consists of a single-stage refrigeration machine designed SHRINK FITTING 
for operation in the range of minus 320 degrees F. to ambient temperatures. 
As a basic cryogenic tool, the cryogenerator is completely adaptable to a 
tremendous number of applications and is in fact being widely used by a host of 
manufacturers and laboratories throughout the country. 
Since the operating principle is identical for all models, specific applications 
of each depend primarily on the capacity of the system rather than the form of 
equipment. The four models, designated A, B, C and D, have respective 
refrigeration capacities of 3,000, 12,000, 40,000 and 160,000 BTU/hr. at 
minus 320 degrees Fahrenheit. 
“Packaged” units consisting of a Model Cryogenerator with 
special headers or accessory equipment are also available. 


Research facilities at the Norelco Cryogenics Application Labora- 
tory are available for consultation on industrial problems involy- 
ing any low temperature application. Write today for illustrated 
brochures and complete technical information. 


AIR LIQUEFIER — This unit liquefies air at atmos- 
pheric pressure — thereby insuring a contaminant-free 
product. Rate of production is 7% liters per hour. 


_ NITROGEN GENERATOR — Comprised of a Model ENVIRONMENTAL CHAMBER — A specially insu- 
A Cryogenerator acting as a nitrogen condenser and lated storage chest is used in conjunction with a Model 
an atmospheric pressure separation column. Rate of A Cryogenerator. This permits the production and 
production is 6 liters per hour over an uninterrupted maintenance of temperatures down to minus 320 de- 
period up to 200 hours. grees F, for environmental testing. - 
Lreleo NORTH AMERICAN PHILIPS COMPANY, Inc. 
CRYOGENICS Division 100 Stevens Avenue, Mount Vernon,N. Y. 
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At IBM, scientists, mathematicians and engineers are currently engaged in many 
unique assignments in advanced physics. For example: 


1. Photoelasticity is being studied to determine strain propagation as a 

_ function of time in solids. This work is designed to clarify laws governing 
impulse strain waves, as well as to determine dynamic stresses that exist 
in impact and other phenomena. 


2. Fundamentals of the magnetic properties of materials are under study. The 
_ immediate interests are spin-wave investigations, particularly those re- 
a. lating to the switching times of remagnetization and spin-wave propaga- 
ees: in polycrystalline materials. 


3. Optical techniques, with organic and inorganic thin films, are being used to 
examine the properties of absorption phenomena in surfaces and substrates. 


You'll find a wide range and diversity of career areas at IBM. Advancement 
is rapid, due in part to the demands of a constantly expanding program of 
research and development, and in part to a long-continued record of promoting 
from within, based on individual merit and achievement. Depending on 
your field of interest and on your abilities, you may work alone or as a member 
of a small team. You’ll be given all the ices assistance you need. 
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CAREERS AVAILABLE 
IN THESE AREAS: 


SOME ASSIGNMENTS 
NOW OPEN INCLUDE: 


Mathematics 
Mechanics 
Nuclear Physics 
Optics 
Semiconductors 
Solid State 


Acoustics 
Aerodynamics 
Astrophysics 
Cryogenics 
Inertial Control 
Magnetics 


B.S. or Advanced Degree in Physics, Electrical 
Engineering, Engineering Physics, or Mathematics 
—plus proven ability to assume a high degree of 
technical responsibility in your sphere of interest. 


Investigation of theoretical device design and ex- 
perimental proof of feasibility for ultrahigh-speed 
transistor; P-N junction technology; surface stud- 
ies; optimization of semiconductor fabrication. 


Creating techniques for solving boundary value 
problems, such as field theory of structures 
(Green’s function methods, for example), or stress 
wave propagation in solids. 


Study of flow and pressure distributions in hydro- 
dynamic transients, to provide information for 
the development of a comprehensive theory of 
short duration (pulse), high Reynold’s number 
fluid-flow systems. 


Determining dynamic behavior of micromagnetic 
fields by applying low-level, high-frequency circuit 
knowledge for small, closely coupled circuits as a 
result of Hall effect studies. 


Study of basic scientific phenomena in the field 


of energy conversion; design of experiments for 
the investigation of electrical discharge in various 
environments and investigation of distribution 
and loss of energy. The objective is to develop 
energy conversion techniques that exceed the 
limited performance and efficiency of existing 
conversion methods. 


Physical-mathematical studies to determine how 
complex control problems may be solved with 
digital computers. Applications in missile systems 
and special-purpose computer systems. 


Development of magnetic devices for use in logic 
switching and computer ‘‘memory.”’ Experience in 
fundamentals of nonlinear magnetics necessary. 


Write, outlining background and interests, to: 
Mr. R. E. Rodgers, Dept. 572K 

IBM Corporation 

590 Madison Avenue, 
New York 22, N.Y. 


IBM. 


tion of hydrocarbons). Volume 2 con- 
tains brief appendixes dealing with the 
method of the activated complex (written 
by M. I. Temkin), quantum-mechanical 
calculations of the activation energy (writ- 
ten in collaboration with N. D. Sokolov), 
and selected additions to Volume 1 which 
were sent by the author to the translator 
after publication of Volume 1. 

To the kineticist, the present volume 
constitutes an invaluable up-to-date com- 
pilation by the author who published a now 
classical treatise on chain reactions nearly 
25 years ago; to the specialist in combus- 
tion or propulsion, the present work offers 
essential reference material in kinetics 
which, after suitable simplification, is 
sometimes useful in the quantitative de- 
scription of combustion processes. Com- 
bustion theory per se is not considered in 
Semenov’s book. 


Ordinary Differential Equations, by Wil- 
fred Kaplan, Addison-Wesley Publish- 
ing Co., Inc., Reading, Mass., 1958, 
534 pp. $8.50. 

Reviewed by E. A. TRABANT 7 
Purdue University 


The first three chapters of the text are 
at an elementary level necessary for stud- 
ents beginning in differential equations. 
However, in these first three chapters ma- 
terial is included which contributes to the 
mathematical maturity of the student who 
has already had a brief introduction to dif- 
ferential equations—that is, existence 
theorems, step-by-step integration, one-to- 
one correspondence, concept of output and 
input, and superposition, etc. 

The text is consistent with the form of 
introducing basic concepts and theory be- 
fore a technique of solution. General 
theorems are stated; their importance and 
applications noted, but most proofs are 
delayed until the last chapter of the book 
which is titled “Fundamental Theory.” 
An exception to this form is the chapter 
on phase-plane analysis where there is a 
good collection of stability theorems with 
proofs. 

A chapter on numerical methods is in- 
cluded in which application of the methods 
to digital computers is mentioned. Ana- 
log computers are referred to, but no work 
or further discussion is presented. 

The approach used in discussing linear 
independent equations in the chapter on 
linear equations of arbitrary order is ex- 
cellent. Some may feel that the summary 
presented on complex functions and the 
Laplace transforms is too brief. However, 
the material included is sufficient for the 
purposes of the text. 

A rather long and detailed treatment of 
simultaneous linear differential equations 
is made. An appendix here introduces 
vectors and the use of matrices. 

The examples and exercises included 
throughout the text are fine, particularly 
those with physical interpretation and 
application. 

The text will be useful to all those inter- 
ested in the application of differential 
equations to engineering problems, par- 
ticularly in the area of servomechanisms 


and automatic control. In addition, it pro- 
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Recent Advances in the Engineering Sci- 
ences, Proceedings of Conference on 
Science and Technology for Deans of 
Engineering, Purdue University, Mc- 
Graw-Hill Book Co., Inc., New York, 
= 1958, 265 pp. $4.75. 

at Reviewed by Cart A. ANDERSON 

~ Institute of Technology (Air University) 


Although this book is made up of the 
proceedings of a conference held in Sep- 
tember 1957, the topics covered are no 
less important today, two years later, than 
when originally presented. Certainly, 
engineering education is changing and 
growing to keep up with new develop- 
ments in science and technology, but the 
changes and growth generally do not occur 
precipitately. Rather, they occur through 
an envolutionary process attended by con- 
siderable discussion, self-study and, per- 
haps, dissent on the part of engineering 
educators. This book can be a definite 
aid to these educators as they work to im- 
prove their curricula. 

The conference was directed toward the 
problems of the engineering college admin- 
istrators whose management duties often 
make it impossible for them to keep up- 
to-date on technical developments, but 
who are, nevertheless, responsible for cur- 
ricular offerings. Each of the technical 
topics was covered in two parts: The first 
a survey of the state of the art in the field, 
and the second a discussion of how develop- 
ments in the field might affect engineering 
education. With different speakers for 
each part, a certain amount of redundancy 
was evident in the papers. Even more 
redundancy would have been acceptable, 
however, if the post-presentation discus- 


sions could have been included in the book. 
The contributions of the administrators 
attending the meeting would probably 
have made quite valuable reading for 
their fellows and their co-workers. Un- 
fortunately, as stated in the preface, it 
was not feasible to include the discussions 
in the book. 

Seven technical fields were touched upon 
by 14 invited authorities in these fields. 
J. G. Truxal surveyed automation and 
automatic control, stressing the growing 
use of computers in control systems, the 
still new adaptive controllers, and the 
types of models used in control design. 
The companion piece by G. 8. Brown was 
directed more toward engineering science 
in general and a possible reorganization of 
engineering education. Operations re- 
search and systems engineering were cov- 
ered by H. H. Goode and D. G. Malcolm 
in papers which went a long way toward 
defining these new fields and locating 
them within the engineering spectrum. 
M. Tribus organized the subject of 
thermodynamics and pointed out the im- 
portance of modern studies of irreversible 
processes; N. A. Hall argued for a better 
grounding in thermodynamics for all engi- 
neers in order that they might better 
assimilate new ideas in their specialties. 
E. R. G. Eckert related the transport 
processes of mass, momentum and heat 
transfer with appropriate limitations and 
quite specifically outlined the problems for 
research. R. B. Bird agreed with Eckert 
in recommending that teaching of engi- 
neering be organized by basic concepts 
rather than specific applications. Bird 
emphasized use of his equations of change 
as a keystone to understanding the trans- 
port processes. S. McLain saw nuclear 
engineering as a broader field than did 
H. J. Gomberg, who kept it more to re- 


actor work. Gomberg’s hope for gaining 
better understanding of energy utilization 
from better understanding of the physical 
mechanism of nuclear energy release was 
interesting. In solid-state physics and 
engineering materials, J. E. Goldman and 
G. Murphy saw a science replacing an 
art in the development of materials. 
Goldman detailed specific technical ad- 
vances, and Murphy detailed specitic 
course requirements. The seventh techni- 
cal field, computer development and 
application, was discussed by S. Alexander 
and J. P. Nash. Alexander emphasiz:d 
data processing properties of compute:s, 
Nash saw a change in engineering educa- 
tion as a result of computers changing 
methods of engineering analysis. 

R. W. Kettler developed an eighth ari a, 
administrative problems in engineering, in 
a paper “Recent Trends in the Cost of 
Engineering Education.” He had several 
suggestions for meeting the problem of ris- 
ing costs while maintaining the quality of 
engineering education. 

As is evident from this listing, maiy 
areas of current engineering interest «are 
covered in the book. The approach uscd 
is such that primary interest in its con- 
tents is to be expected from educators. [t 
is certainly worth their time to consider 
the ideas presented. Many may argue 
that the unified organization of the engi- 
neering sciences which is obvious to the 
authority may be too difficult for the 
beginner, although this reviewer would 
have appreciated more of the unified ap- 
proach as a student. For engineers out of 
the educational sphere, the book can 
serve as a good tool for broadening their 
understanding of modern technology, and, 
it is hoped, stimulate their interest in the 
engineering colleges from which their 
future co-workers must come. 


Injector and method for fuel injection 
(2,881,588). W. H. Goss, Silver Spring, 
Md., assignor to the U.S. Navy. 

In a ramijet-type guided missile, a 
fuel injector having a chambered wind- 
shield and an apertured strut. Fuel 
discharged from the windshield and strut 
against the airstream is mixed to produce 
a uniform combustible mixture. 


Jet engine silencer (2,882,991). 
Killian, Hayesville, N.C. 

A bulbous housing with a rear opening 
surrounding an exhaust tube. The 


C. R. 


Eprror’s Norte: Patents listed above 
were selected from the Official Gazette of 
the U.S. Patent Office. Printed copies 
of patents may be obtained from the 
Commissioner of Patents, Washington 25, 
D. C., at a cost of 25 cents each; design 
patents, 10 cents. 


882 


inner surface of the housing is transversely 
corrugated to provide baffles for reflecting 


sound waves. A portion of the exhaust 
tube in the housing is spaced to form a 
low pass acoustical filter for attenuating 
the sound waves. 

Low drag exhaust silencer (2,882,992), 
G. F. Hausmann, Glastonbury, Conn., 
assignor to United Aireraft C orp. 

Outer casing forming a main passage 
for a jet exhaust. A central inner body 
terminates in a hollow conical member 
rotating within the casing. Small openings 
in the member emit fluid flow at a reduced 
noise level. 

Silencer (2,882,993). K. 
Arlington, Texas, assignor 
Vought Aircraft, Inc. 

Fluid flow conduit in an air condi- 
tioning system having a jet of air near 
sonic velocity. Equally spaced flat 


T. Murty, 
to Chance 


plates normal to the conduit longitudinal 
axis have orifices which are larger in 


George F. McLaughlin, Contributor 


number and size than those of the pre- 
ceding plate, reducing the acoustical 
noise level in the system. 
Powerplant incorporating a dynamic com- 
pressor (2,883,828). A. R. Howell, 
Farnborough, England, assignor to the 
British Government. 

Rocket combustion chamber discharging 
into a hollow rotor with reaction nozzles 


positioned tangentially on the rotor 
blades. Reaction of gases from the 
nozzles drives a compressor mounted 


within a tubular main duct having a 
propulsion nozzle downstream. 

Fuze (2,883,933). R. B. King, Altadena, 
Calif., assignor to the U.S. Navy. 

Gas pressure from a rocket motor 
releases a firing pin to an armed position. 
A spring loaded acceleration responsive 
means intercepts a shutter separating the 
pin from the detonator to defer its move- 
ment for a limited time to the armed 
position. 
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GET IT FAST, 
GET IT CLEAR 
ON LINO-WRIT 4 


Du Pont Lino-Writ 4 is the fastest, toughest, whitest photorecording 
paper you can use. It accommodates writing speeds as high as 5000 
cps at maximum amplitude; has exceptional exposing and processing 
latitude; its ultra-thin, all-rag stock has unusual wet or dry strength. 
Our other outstanding papers include: 
Lino-Writ 1. Orthochromatic paper designed for test requirements 
at lower writing speeds. 
Lino-Writ 2. Moderately high-speed orthochromatic paper for use in 
intermediate speed recording requirements. 
Lino-Writ 3. For high-frequency oscillographic traces without loss 
of detail. 

These three papers are available in two types: Type B, standard- 


This advertisement was prepared exclusively by Phototypography. 


Better Things for Better Living . . . through Chemistry 


weight stock when greater opacity is re- 
quired and Type W, a thin 100% all-rag stock. 


Lino-Writ Rapid Processing Chemicals 
Kit is Du Pont’s high-speed, low-odor for- 
mulation designed for rapid stabilization 
processing of all Du Pont photorecording 
papers. 

For more information on our complete 
photorecording line, write: E. I. du Pont 
de Nemours & Co. (Inc.), Photo Products 
Department, 2430-A Nemours, Wilming- 
ton 98, Delaware. In Canada: Du Pont of 
Canada Limited, Toronto. 


REG, U.S. PAT.OFE 
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Vertical takeoff and landing apparatus 
(2,884,633). B.Stahmer, Ohaha, Neb. 


~~ 


A pair of multibladed rotors set into 
the wing of an aircraft. Rotors are driven 
by the jet stream of a jet engine striking 
buckets on the rim. When the rotors are 
not operating, they close the circular 
opening in the wing, providing additional 
lift. Forward thrust is obtained during 
flight when the engines are mechanically 
moved to bypass the rotors. 


Afterburner igniter (2,885,857). W. F. 
Hemlock, W. Hartford, Conn., assignor 
to the U. 8. Air Force. 

Pump housing with two chambers. 
Fuel is pumped from one chamber to the 
engine while fuel is supplied to the other. 
Reversing means then pumps fuel from 
the second chamber while fuel is supplied 
to the first chamber. Turbine discharge 
pressure controls the action. 


Aircraft jet system with adjustable nozzle 
(2,886,262). P. J. Fletcher, Farnborough, 
England, assignor to Power Jets (Research 

and Development) Ltd. anil 


2,886,262 


Bon! 

Nozzle for diverting a jet discharge 
smoothly from ene direction to another. 
Jet can be diverted from rearward to 


downward to produce an upthrust for 
assisting takeoff or slow speed landing. 


Sound absorbing housing (2,884,086). 
J.T. Welbourn, Pittsburgh, Pa. 

Jet engine testing structure comprising 
spaced rows of perforated hollow cylinders 
containing loose granular material to 
absorb and deflect sound. 


Fluid jet sustained and controlled air- 
~~ (2,887,284). W. P. Perry, Trenton, 
IN. od. 

Air pressure fan mounted in a convey- 
ance having a roof. Nozzles discharge 
air pressure across the opposite edges of 
areas of the roof surface to create a 
vortex surrounding and a vacuum within 
the areas to create lift. 
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Water-cooled pitching moment balance 
(2,885,890). L. L. Liccini and R. H. 
Cornett, Silver Spring, Md., assignors 
to the U. S. Navy. 

Strain-gage balance for a hypersonic 


wind tunnel in which air speeds of Ma 

to 10 are reached. Cooling water 
passes through the balance to keep 
at the temperature at which it was 
calibrated. 


Windows for admitting solar radiation 
(2,888,007). H. Z. Tabor, Jerusalem, 
Israel, assignor to the State of Israel. 


Solar trap with transparent plate oa, : 


having on the surface a number of lenses 
adapted to focus the solar radiation 
into an enclosed space where the foci of 
all lenses are located in one plane. 


Method of preparing globular propellant 
powder (2,888,713). R. L. Cook and 
E. A. Andrew, Berkeley, Mo., assignors 
to Olin Mathieson Chemical Corp. 

A lacquer of nitrocellulose subdividea 


into ciscrete particles in suspension in a> 
nonsolvent liquid. Suspension is main- — 


tained under superatmospheric pressure 
sufficient to prevent boiling of the solvent 
for the nitrocellulose. 


Partial pressure suit (2,886,027). J. P. 
Henry, Yellow Springs, Ohio, assignor to 


the U. S. Air Force. 


Suit supplying uniform all-over pressure 
to an aviator to balance the necessary 


rise in intrathoracic pressure applied 
by a special helmet, and to prevent blood 
pooling and plasma loss. Suit prevents 
dangerous shifts of liquids and gases 
within the body. 


High altitude balloon for meteorological 
use (2,886,263). D. M. Ferguson, Mar- 


blehead, Mass., assignor to the U. S. Air 
Force. 


Inclusion of a relatively small inner 
balloon within the envelope of a high al- 
titude balloon. Means between the enve- 
lopes permits gradual transfer of the 
initially-charged gas from the inflated 
inner envelope to the initially-deflated 
outer envelope during ascension. 


Stall roll control for vertical takeoff 
airplane (2,886,264). D. B. Seager, Bur- 
bank, Calif., assignor to Lockheed Air- 
craft Corp. 


2,886,264 


Control for maintaining wing loading 
substantially symmetrical at high angles 
of attack at the threshold of a fully 
developed stall condition. Wing ports 
direct fluid onto the upper surface of the 
wing to establish a high lift critical angl« 
of attack. 
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Jet Propulsion Engines 


Advanced Propulsion Systems. A Pre- 
liminary Study, by H. F. Calcote, Aero- 
Chem Res. Labs., Inc., Tech. Pub. 6, July 
1958, 91 pp. 

Factors Influencing the Application of 
Unconventional Rocket Engine Systems, 
hy George S. Gill and Leo Kusak, SAE 
Preprint 548, March 31—April 3, 1959, 15 
pp. 

A Graphical Solution for Solid-Pro- 
pellant Motor Design, by Robert N. 
Witte, Army Rocket and Guided Missile 
Agency, ARGMA Rep. 6R15, Jan. 
i959, 23 pp. 

Thrust Magnitude Control of Solid 
Propellant Rocket Motors by Mechanical 
Means, by J. S. Gates and 8. L. Pinto, 
SAE Preprint 59T, March 31-April 3, 
1959, 19 pp. 

Non-Air Breathing Auxiliary Power- 
plants, by D. L. Cochran, A. T. Biehl, 
D. R. Sawle, M. R. Gustavson and A. M. 
Taylor, SAE Preprint 53T March 31- 
(pril 3, 1959, 23 pp. 

Optimum Exhaust Velocity Program- 
ming and Propulsion Efficiency, by Robert 
Fox, J. Astron. Sci., vol. 6, no. 1, Spring 
1959, pp. 13-16. 

Combustion Instability in Liquid Pro- 
pellant Rocket Motors, by David T. 
Harrje, Princeton University, Dept. Aeron. 
Engng., Rep. 216x, May 1959, 20 pp. 

Theoretical Studies of the Performance 
of Heat Engines Using Pressure Waves, 
by J. C. Logan Jr., Cornell Aeronautical 
Laboratory Rep. DD-799-A-1, Dec. 1954, 
196 pp. (Declassified by authority of 
ONR: 429: CAP: III, June 3, 1959.) 

Experimental Investigation of a Reso- 
nant Wave Engine, by A. Hettzberg and A. 
Russo, Cornell Aeronautical Laboratory, 
Rep. DD-799-a-2, Dec. 1954, 54 pp. 
(Declassified by authority of ONR: 429: 
CAP: ITI, June 3, 1959.) 

Prospects for Thermoelectricity, by H. 
Gettings, Missiles and Rockets, June 22, 
1959, pp. 29-31. 

Rocket Turbines, by Otis E. Lancaster 
and Carlton J. Bates, 1958 Gas Turbine 
Report ASME, 1959, 141 pp., pp. 66-75. 

Nuclear Progress, by Richard P. God- 
win and Edward Dennison, 1958 Gas 
Turbine Rep. ASME 1959, 141 pp., pp. 
134-141. 

Rocket Thrust Termination Transients, 
by Howard C. Rodean, ARS Journat, 
vol. 29, June 1959, pp. 406-409. 

Suppression of Combustion Instability 
in the Plunger Jet Rocket Engine, by 
Herbert Shieber, ARS Journat, vol. 
29, June 1959, pp. 446-447. 

Development of H.O, Rocket Motors, 
by D. Hurden, Raketentechnik und Raum- 
fahrtforschung, vol. 3, no. 1, Jan.—March 
1959, pp. 9-13. (In German.) 


Eprror’s Note: Contributions from Pro- 
fessors E. R. G. Eckert, J. P. Hartnett, T. 
F. Irvine Jr. and P. J. Schneider of the 
Heat Transfer Laboratory, University of 
Minnesota, are gratefully acknowledged. 
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The James Forrestal Research Center, Princeton University 


Aerodynamics of Jet 
Propelled Missiles 


Slender Body of Revolution of Hyper- 
sonic Speeds—‘‘An Approximation to the 
Determination of the Flow Field,” by R. 
Capiaux, General Research in Flight 
Sciences, Jan. 1958-Jan. 1959, vol. IT, 
Fluid Mechanics, Inviscid Flows, Lockheed 
Aircraft Corp., Missiles and Space Divi- 
sion, Rep. LMSD-48381, Jan. 1959, 163 
pp., pp. 101-108. 

A Linearized Method of Characteristics 
for Determining Body Aerodynamics 
Forces at Angle of Attack, by J. Der, 
General Research in Flight Sciences, Jan. 
1958-Jan. 1959, vol. II, Fluid Mechanics, 
Inviscid Flows, Lockheed Aircraft Corp., 
Missiles and Space Division, Rep. LMSD- 
48381, Jan. 1959, pp. 109-138. 

Use of Drag Modulation to Reduce 
Deceleration Loads During Atmospheric 
Entry, by Richard L. Phillips and Clarence 
B. Cohen, ARS Journat, vol. 29, June 
1959, pp. 414-422. 

Effect of Thrust Termination Process 
Upon Range Dispersion of a Ballistic 
Missile, by Arnold J. Kelly, ARS Jour- 
NAL, vol. 29, June 1959, pp. 432-440. 

Theoretical and Experimental Analysis 
of a Cowling as a Means of a Drag Re- 
duction for an Axisymmetric Center Body, 
by Marian Visich Jr. and Anthony Mar- 
telluci, ARS Journau, vol. 29, June 
1959, pp. 447-449. 

Electrically Charged Bodies Moving 
in the Earth’s Magnetic Field, by W. W. 
Fain and B. J. Greer, ARS JouRNAL, 
vol. 29, June 1959, pp. 451-453. 

Fluid-Dynamic Properties of Some 
Simple Sharp and Blunt Nosed Shapes at 
Mach Numbers from 16 to 24 in Helium 
Flow, by Arthur Henderson Jr. and 
Patrick J. Johnston, NASA Memo 5-8- 
59L, June 1959, 63 pp. 

Phugoid Motion of Hypervelocity Ve- 
hicles, by G. Breaux, Symposium on High- 
Speed Aerodynamics and Structures, 3rd, 
San Diego, March 25-27, 1958 (Unclassi- 
fied Papers), vol. 3, pp. 11-16. 

Aerodynamic Coefficients on Yawed 
Slender Bodies at Hypersonic Speeds, by 
L. Trilling and J. W. Clark, Symposium 
on High-Speed Aerodynamics and Struc- 
tures, 3rd, San Diego, March 25-27, 1958 
(Unclassified Papers), vol. 3, pp. 71-87. 

Inviscid Hypersonic Flow Around Blunt 
Bodies, by Marcel Vinokur, General Re- 
search in Flight Sciences, Jan. 1958- 
Jan. 1959, vol. II, Fluid Mechanics, In- 
viscid Flows, Lockheed Aircraft Corp., 
Missiles and Space Division, Rep. LMSD 
48381, Jan. 1959, pp. 1-39. 

Supersonic Flow Around Blunt Bodies, 
an Extension of Van Dyke’s Technique, 
by M. Vinokur and R. W. Sanders, 
General Research in Flight Sciences, Jan. 
1958-Jan. 1959, vol. IT, Fluid Mechanics, 
Inviscid Flows, Lockheed Aircraft Corp., 
Missiles and Space Division, Rep. LMSD 
48381, Jan. 1959, pp. 40-56. 

Surface Pressure Distribution on Blunt 
Bodies at Supersonic Speeds, by R. M. 
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RESEARCH 
OPPORTUNITIES 
In 


SPACE 
SCIENCES 


The Space Technology 
Operations of Aeronu- 
tronic has immediate 
need for engineers and 
scientists who are inter- 
ested in working in the 
stimulating and highly 
diversified field of space 
sciences. This West 
Coast division of Ford 
Motor Company has the 
newest facilities and most 
advanced equipment for 
carrying out highly tech- 
nical work—challenging 
creative work that is ex- 
ceptionally rewarding to 
qualified men. 


Positions are at Aero- 
nutronic’s new $22 mil- 
lion Research Center, be- 
ing completed at New- 
port Beach, Southern 
California. Here, over- 
looking famous Newport 
Harbor and the Pacific 
Ocean, relaxed California 
living can be enjoyed free 
of big-city congestion, 
yet most of the important 
cultural and educational 
centers are just a short 
drive away 


AREAS OF INTEREST 
VEHICLE TECHNOLOGY 


Aerodynamic design and 
testing 

Rocket Nozzle and re-entry 
materials 

High temperature chemical 
kinetics 


Combustion thermodynamics 

High temperature structural 
plastics and ceramics 

Advanced structures 


SYSTEMS DEVELOPMENT 


Aerothermodynamics 

Re-entry programs 

High temperature heat 
transfer 

Penetration systems 

Hyper environmental test 
systems 


ELECTRONICS AND 
ASTRO SCIENCES 


Astro navigation 

Space communications and 
communication satellites 

Instrumentation, telemetering 
and data reduction 

Space environmental physics 

Advanced techniques and sys- 
tem studies 


Qualified applicants are invited to 
send resumes and inquiries to 
Mr. R. W. Speich, Aeronutronic, 
Dept. 19, Box 451, Newport 
Beach, California. 


AERONUTRONIC 


A Division of 
FORD MOTOR COMPANY 
Newport Beach 
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e Mark, General Research in Flight Sciences, 


an. 1958-Jan. 1959, vol. II, Fluid Me- 
chanics, Inviscid Flows, Lockheed Air- 
craft Corp., Missiles and Space Division, 
Rep. LMSD-48381, Jan. 1959, pp. 57-62. 

An Extension of Second-Order Shock 
Expansion Theory, by R. Capiaux, General 
Research in Flight Sciences, Jan. 1958- 
Jan. 1959, vol. II, Fluid Mechanics, Invis- 
cid Flows, Lockheed Aircraft Corp., Missiles 
and Space Division, Rep. LMSD-48381, 
Jan. 1959, pp. 63-99. 

Variational Approach to the Re-Entry of 
a Ballistic Missile—Part I, by Angelo 
Miele and Carlos R. Cavoti, Purdue Uni- 
versity, School of Aeron. Engng., Rep. 
A-59-1, Jan. 1959, 57 pp. 

Surface Pressure Distribution at Hyper- 
sonic Speeds for Blunt Delta Wings at 
Angle of Attack, by Marcus O. Creager, 
NASA Memo 5-12-59 A, May 1959, 15 pp. 

Aerodynamic Field Near the Stagnation 
Point of a Blunt Symmetrical Airfoil in 
Hypersonic Flow With Near Unity, by 
Aldo Muggia, Torino, Politechnico Labo- 
ratorio di Meccanica Applicata, TN 13 
(AFOSR TN 59-339; ASTIA AD 213660), 
Feb. 1959, 15 pp. 

The Motion of an Earth Satellite on Re- 
Entry to the Atmosphere, by T. Non- 
weiler, Astronautica Acta, vol. 5, no. 1, 
1959, pp. 40-62. 

Nose-Drag in Free roo Flow and 
Its Minimization, by H. S. Tan, J. Aero/- 
_ Sci., vol. 26, June 1959, pp. 360- 


Heat Transfer 
and Fluid Flow 


Study of Supersonic Flows With Chemi- 
cal Reactions, II: Supersonic Nozzle Flow 
with a Reacting Gas Mixture at High 
Reaction Concentration, by Peter P. 
Wegener, Calif. Inst. Tech., Jet Prop. 
Lab. Progr. Rep. 20-370, Dec. "1958, 29 pp. 

Experimental Investigation of Coaxial 
Jet Mixing of Two Subsonic Streams at 
Various Temperature, Mach Number, and 
Diameter Ratios for Three Configurations, 
by Richard R. Burley and Lively Bryant, 
NASA Memo 12-21-58E, Feb. 1959, 29 
pp. 

An Introductory Discussion of Magneto- 
hydrodynamics, by A. R. Kantrowitz 
and H. E. Petschek, Avco Mfg. Corp., 
Avco Res. Lab. Div., Res: Rep. 16, May 
1957, 15 pp. 

Supersonic Two-dimensional Magneto- 
hydrodynamic Flow, by F. Fishman, J. 
Lothrop, R. Patrick and H. Petschek, 
Avco Mfg. Corp., Avco Res. Lab. Div., 
Res. Rep. 39, Feb. 1959, 49 pp. 

On the Flow of Gases Under Free 
Molecular Conditions, by D. Roger Willis, 
Princeton Univ., Dept. Aeron. Enqnq., Rep. 
442 (AFOSR-T N-58-1093: ASTIA AD 
207594), Dec. 1958, 50 pp. 

One Dimensional Compressible Flow 
Tables for a Perfect Gas in Mach Num- 
ber Increment of 0.001 (k = 1.26), Isen- 
tropic Process, by C. Hoebich, Army 
Rocket and Guided Missile Agency, Ordn. 
Missile Labs., Div. ARGM A, TN 1H1N-1, 
Jan. 1959, 84 pp. 

Adiabatic Wall Temperature Due to 
Mass Transfer Cooling With a Combustible 
Gas, by George W. Sutton, ARS Journat, 
vol. 29, Feb. 1959, pp. 136-137. 

Energy Transfer at a Chemically React- 
ing or Slip Interface, by Sinclaire M. Scala 
and George W. Sutton, ARS JouRNAL, 
vol. 29, Feb. 1959, pp. 141-143. 

Effect of Unequal Wall Roughness on 
Flow Between Flat Plates, by H. N. 
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McManus Jr., ARS Journat, vol. 29, 
Feb. 1959, pp. 144-145. 

Scale of Separation Phenomena in 
Liquids Under Conditions of Nearly Free 
Fall, by E. T. Benedikt, ARS Journat, 
vol. 29, Feb. 1959, pp. 150-151. 

Hypersonic Aerodynamics, Part I, by 
Wallace D. Hayes, AsTrRoNAuTICs, vol. 
4, March 1959, pp. 24-25, 72, 74. 

Arc-heated Plasma for Laboratory Hy- 
personics, by T. R. Hogness, ASTRONAU- 
Tics, vol. 4, March 1959, pp. 40-42, 47. 

On Initial Conditions in Hydromag- 
netics, by G. S. S. Ludford, Cambridge 
Phil. Soc., Proc., vol. 55, Part 1, Jan. 
1959, pp. 141-143. 

An Investigation on Plasma in External 
Magnetic Fields, I: Steady State, by G. 
Schmidt, Nuovo Cimento, vol. 10, no. 1, 
Oct. 1, 1958, pp. 55-67. 

An Investigation on Plasmas in External 
Magnetic Fields, II: Varying Fields, by 
G. Schmidt, Nuovo Cimento, vol. 10, no. 
4, Nov. 16, 1958, pp. 659-674. 

Statistical-Mechanical Theory of Trans- 
port in Fluids, by Hazime Mori, Phys. Rev., 
vol. 112, Dec. 15, 1958, pp. 1829-1842. 

Estimated Performance of Radial-flow 
Exit Nozzles for Air in Chemical Equi- 
librium, by Gerald W. Englert and Fred D. 
Kochendorfer, NASA Memo 1-5-59E, 
Feb. 1959, 53 pp. 

Internal Performance of Two-dimen- 
sional Wedge Exhaust Nozzles, by William 
T. Beale and John H. Povolny, NACA 
RM E56K29b, Feb. 1957, 35 pp. (Declas- 
sified by authority of NASA Pub. An- 


nouncement, vol. 3, 2/4/59, p. 26). 


Effect of Outer-Shell Design on Per- 


_ formance Characteristics of Convergent- 


plug Exhaust Nozzles, by H. George Krull 
and William T. Be: ale, NACA RM 
E54K22, April 1955, 18 pp. (Declassified 
by authority of NASA Pub. Announce- 
ment, vol. 3, 2/4/59, p. 22). 

T-F Theory of Electron Emission in 
High-current Arcs, by T. H. Lee, J. 
Appl. Phys., vol. 30, Feb. 1959, pp. 166- 
171. 


An Application of Irreversible Thermo- 
dynamics to the Study of Diffusion, by 
Richard W. Laity, J. Phys. Chem., vol. 
63, Jan. 1959, pp. 80-82. 

Hydromagnetic Equations for Two Iso- 
topes in a Com oy! Ionized Gas, by 
Joseph Slepian, Phys. Rev., vol. 112, Dec. 
1, 1959, pp. 1441-1444. 

Plasma Oscillations with Diffusion 
Velocity Space, by A. Lenard and Ira B. 
Bernstein, Phys. Rev., vol. 112, Dec. 1, 
1958, pp. 1456-1459. 

Microwave Propagation in Hot Mag- 
neto-plasmas, by J. E. Drummond, Phys. 
Rev., vol. 112, Dee. ‘: 1958, pp. 1460-1464. 

Transverse Plasma Waves and Plasma 
Vortices, by O. Buneman, Phys. Rev., vol. 
112, Dec. 1, 1958, pp. 1504-1512. 

Microwave Emission from High-Tem- 
neg Plasmas, by David B. Beard, 

hys. Rev. Letters, vol. 2, Feb. 1, 1959, pp. 
81-82. 

Fast Ion Heating, by Dary! Reagan, 
Phys. Rev. Letters, vol. 2, Feb. 1, 1959, pp. 
82-83. 

Spontaneously Growing Transverse 
Waves in a Plasma Due to an Aniso- 
tropic Velocity Distribution, by iE rich 8. 
Weibel, Phys. Rev. Letters, vol. 2, Feb. 1, 
1959, pp. 83-84. 

The Generalized Thermal Energy in 
Boundary Layer Type Flows with Chemi- 
cal Reaction, by S. F. Shen, Zeitschrift fur 


Angewandte Mathematik und Mechanik, 
vol. 38, Nov.-Dec. 1959, pp. 431-436. (In 
English. ) 


Heat Transfer in Reacting Systems: I. 
Heat Transfer to N.O, Gas under Turbu- 
lent Pipe Flow Conditions, by Walter F. 
Krieve and David M. Mason, Calif. Inst. 
Tech., Jet Prop. Lab., Progr. Rep. 20-366, 
Nov. 1958, 13 pp. 

A Review of Two Phase Heat Transfer 
(1935-1957), by J. G Collier, Gt. Brit. 
Atomic Energy Res. Estab., AERE CE R 
2496, 1958, 93 pp. 

Generation of Supersonic Dissociated 
and Ionized Non-equilibrium Streams, |! 
D. E. Rosner and H. F. Caleote, Aero- 
Chem. Res. Labs., Inc., TM-10 (AFOS/- 
TN-58-1080; ASTIA AD 207590), Oct. 
1958, 30 pp. 

Rayleigh’s Problem in Hydromagnetic: . 
Impulse Motion of a Pole Piece, by G. S. ». 
Ludford, Maryland Univ., Inst. for Flu 
Dynamics and Appl. Math., TN BN-151 
(AFOSR-TN 58-1073; ASTIA A!) 
207242), Nov. 1958, 22 pp. 

Sublimation in a Hypersonic Environ- 
ment, by Sinclaire M. Scala, Gen. Ele. 
Co., Missile and Space Vehicle Dept., Tec’. 
Info. Series R58SD289, 1958, 45 pp. 

Gas Dynamic Mixing Processes; High 
Thrust through Jet Injection, by O. Lutz, 
Wissenschaftlichen Geselischaft fiir Lufi- 
fahrt, Jahrbuch, 1957, Braunschweig, 
Friedr. Vieweg & Sohn, 1959, pp. 374-380). 
(In German. ) 

Heat Transfer and Boundary Layer 
Formation in Sweated Materials, by |’. 
Meyer-Hartwig, Wissenschaftlichen Gesel'- 
schaft fiir Luftfahrt, Jahrbuch, 1957, 
Braunschweig, Friedr. Vieweg & Sohn, 
1958, pp. 439-468. (In German.) 

Calculation of the Effective Ionization 
Potential ia Plasmas, by H. Rother, Ani 
Phys., vol: 2, no. 5-6, 1958, pp. 326-338. 
(In German. ) 

Multiple Recording Manometer for 
Rapidly Changing Pressures in Flow 
Measurement, by Walter Wuest, Arch 


fiir Tech. Messen., no. 274, Nov. 1958, pp. 


225-228. (In German.) 

The Structure of a 

Shock in Steady Plane Motion, by G. S. 5. 
Ludford, J. Fluid Mech., vol. 5, Jan. 1959, 
pp. 67-80. 

On the Structure of Jets from Highly 
Underexpanded Nozzles into Still Air, by 
T. C. Adamson Jr. and J. A. Nicolls, 
Mich. Univ., Engng. Res. Inst., Rep. 
2397-5-F, Feb. 1958, 32 pp. 

Hydromagnetic Instabilities of a Cylin- 
drical Gas Discharge, II: Influence of 
Viscosity, by R. J. Taylor, Gt. Brit., 

Atomic Energy Res. Estab., AERE T/R 
1888, 1958, 38, pp. 

On the Surfaces of Discontinuity in the 
Magneto-fluid-dynamics, by L. G. Napoli- 
tano, Aerotecnica, vol. 38, Aug. 1958, pp. 
210-220. (In Italian.) 

Thermal Conductivity of Multicom- 
ponent Gas Mixtures, by Charles Mucken- 
fuss and C. F. Curtiss, J. Chem. Phys., 
vol. 29, Dec. 1958, pp. 1273-1277. 

Statistical Mechanical Theory of the 
Thermal Conductivity of Binary Liquid 
Solutions, by Richard J. Bearman, J. 
Chem. Phys., vol. 29, Dec. 1958, pp. 1278- 
1286. 

Transport Processes in Binary Liquid 
Systems, by R. R. Irani and A. W. Adam- 
son, J. Phys. Chem., vol. 62, Dec. 1958, pp. 
1517-1520. 

On the Structure of Jets from Highly 
Underexpanded Nozzles into Still Air, by 
T. C. Adamson Jr., J. Aero/Space Sci., 
vol. 26, Jan. 1959, pp. 16-24. 

International Conference on the Con- 
densed State of Simple Systems (The 
Liquid State), N wovo Cimento, Supplement, 
3d Quarter, no. 1, 1958, pp. 1-368. 
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Wave Theory of Heat Transfer in Film 
Boiling, by Yan Po Chang, J. of Heat 
Transfer (ASME Trans.) vol. 81, Series 
C, no. 1, Feb. 1959, pp. 1-12. 

A Boundary-Layer Treatment of Lami- 
nar Film Condensation, by E. M. Sparrow 
and J. L. Gregg, J. of Heat Transfer 
(ASME Trans.), vol. 81, series C, no. 1, 
Feb. 1959, pp. 13-18. 

Transient Heat Transfer for Laminar 
Forced Convection in the Thermal En- 
trance Region of Flat Ducts, by R. Siegel 
and E. M. Sparrow, J. of Heat Transfer 
(ASME Trans.), vol. 81, series C, no. 1, 
Feb. 1959, pp. 29-36. 

Generalized Correlation of Boiling Heat 
Transfer, by 8. Levy, J. of Heat Transfer, 
(ASME Trans.), vol. 81, series C, no. 1, 
Feb. 1959, pp. 37-42. 

Heat Transfer to a Boiling Liquid- 
Mechanism and Correlations, by Kurt 
Engelberg-Forster and R. Grief, J. of 
Heat Transfer (ASME Trans.), vol. 81, 
series C, no. 1, Feb. 1959, pp. 43-53. 

Transient Heat Conduction in Elliptical 
Plates and Cylinders, by E. T. Kirpatrick 
and W. F. Stokey, J. of Heat Transfer 
(ASME Trans.), vol. 81, series C, no. 1, 
Feb. 1959, pp. 54-60. 

Free Convection, Forced Convection, 
and Acoustic Vibrations in a Constant 
Temperature Vertical Tube, by T. W. 
Jackson, W. B. Harrison and W. C. Bote- 
ler, J. of Heat Transfer (ASME Trans.), 
vol. 81, series C, no. 1, Feb. 1959, pp. 
68-74. 

Summary-Selected References on Gas 
Turbines, Regenerative Cycles, Rotary 
Regenerators, and Associated Heat Ex- 
change and Pressure-Drop Information, 
by W. E. Hammond, J. of Heat Transfer 
(ASME Trans.), vol. 81, series C, no. 1, 
Feb. 1959, pp. 75-78. 

Heat Transfer and Friction for Fluids 
Flowing Over Surfaces at High Temper- 
atures and High Velocities, by R. G. 
Deissler and A. L. Loeffler, Jr., J. of Heat 
Transfer (ASME Trans.), vol. 81, series 
C, no. 1, Feb. 1959, pp. 89-94. 

Thermal Relaxation in Gases: A New 
Method of Analysing the Experimental 
Data, by A. C. J. Johnson, Phys. Soc., 
London, Proc., vol. 73, Feb. 1959, pp. 
273-279. 

Flow of Heat in a Composite Solid, 
by J. P. Newcomb, Brit. J. of Appl. Phys., 
vol. 10, May 1959, pp. 204-205. 

Theoretical Structure of Plasma Equa- 
tions, by M. N. Rosenbluth and N. 
Rostoker, Phys. Fluids, vol. 2, Jan.—Feb. 
1959, pp. 23-30. 

Propagation of Disturbances at High 
Frequencies in Gases, Liquids, and 
Plasmas, by H. 8. Green, Phys. Fluids, 
vol. 2, Jan.—Feb. 1959, pp. 31-39. 

Experimental Demonstration of Hydro- 
magnetic Waves in an Ionized Gas, by 
G. A. Sawyer, P. L. Scott and T. F. Strat- 
ton, Phys. Fluids, vol. 2, Jan.—Feb. 
1959, pp. 47-51. 

On the Confinement of a Plasma by 
Magnetostatic Fields, by Erich 8. Weibel, 
Phys. Fluids, vol. 2, Jan.-Feb. 1959, pp. 
52-56. 

Magnetohydrodynamic Flow in a Shock 
Tube, by M. Mitchner, Phys. Fluids, vol. 
2, Jan.-Feb. 1959, pp. 62-71. 

Some Compression Waves in Plasma, 
by Daryl Reagan, Phys. Fluids, vol. 2, 
Jan.-Feb. 1959, pp. 93-94. 

High-Current Ion Source, by Russell G. 
Meyerand Jr. and Sanborn C. Brown, 
Rev. Sci. Instr., vol. 30, Feb. 1959, pp. 
110-111. 
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Interested in Systems Engineering? 


) ) 


...and TOTAL systems 
in which the big bird and support 
equipment may rank only 
as a component. 


This difference between systems can make 
a big difference in your career 


IF YOU ARE QUALIFIED and interested in contributing to programs of 
“total” scope, it will be of value to you to investigate current oppor- 
tunities with General Electric’s DEFENSE SYSTEMS DEPT., whose 3 
work lies primarily in providing total solutions to large scale defense , 
problems of the next 5, 10 and 20 years. _ 


= 


The work here lies almost entirely in the areas of systems engi- 
neering and systems management. ; 


Inquire about these positions: 


Guidance Equation Engineers Systems Test Evaluation Engineers 

Systems Logistics Engineers Engineering Psychologists 

Electronic Systems Radar Equipment Engineers 
Management Engineers Weapons Analysis Engineers 


Operations Analysis Engineers Weapons Systems Integration 
Systems Program Engineers Engineers 
Data Processing Engineers Engineering Writers 


oo Forward your confidential resume at an early date. 


Whereas the growth potential is evident — both for DSD and the 

90 999999))) engineers who join us — the positions we fill during these early 
CL NSO months will carry significant “ground-floor” benefits. 
Write fully to Mr. E. A.Smith,Room 11-M 


DEFENSE SYSTEMS DEPARTMENT 


A Department of the Defense Electronics Division 


«GENERAL ELECTRIC 


300 South Geddes Street, Syracuse, N. Y. 
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Combustion, Fuels 
Propellants 


Burning Mechanism of Ammonium 
Perchlorate Propellants, by Martin Sum- 
merfield, G. S. Sutherland, M. J. Webb, 
H. J. Taback and K. P. Hall, ARS Pre- 
print 737-58, Nov. 1958, 21 pp. 

The Dependence of Linear Burning 
Rate upon Pressure for Ammonium Per- 
chlorate-Polyester Resin Composite Solid 
Propellant, by Maurice J. Webb, Prince- 
ton Univ., Dept. Aeron. Engng., (M.S. 
Thesis), May 1958, 10 figs., 40 pp. 

Theoretical Combustion Performance of 
Several High-energy Fuels for Ramjet 
Engines, by Leonard K. Tower, Roland 
Breitwieser and Benson E. Gammon, 
NACA Rep. 1362, 1958, 18 pp. 

Combustion of Various Highly Reac- 
tive Fuels in a 3.84 by 10-inch Mach 2 
Wind Tunnel, by H. Allen Jr. snd E. A. 
Fletcher, NASA Memo 1-15-59E, April 
1959, 28 pp. 

Analysis of Ammonium Nitrate III, 
Assay, by Eugene A. Burns, Calif. Inst. 
Tech., Jet Prop. Lab., Progr. Rep. 20-365, 
Sept. 2, 1958, 15 pp. 

Some Thermal Stress Design Data for 
Rocket Grains, by M. L. Williams, ARS 
JOURNAL, vol. 29, April 1959, pp. 260-266, 

Turbulent Flame Studies in a Two- 
_Open Burner, by J. H. 
Grover, E. N. Fales and A. C. Scurlock, 
Jou RNAL, vol. 29, April 1959, pp. 
275 982. 

Investigation of the Ignition of a Gas by 
an Electric Spark, by D. R. Lintin 
and E. R. Wooding, Brit. J. Appl. Phys., 
vol. 10, April 1959, pp. 159-165. 

The Heats of Combustion of Organic 
Compounds of Nitrogen, by F. W. Evans, 


D. M. Fairbrother and H. A. Skinner, 
Trans., Faraday Soc., vol. 55, March 
1959, pp. 399-403. 

Burning Mechanism of Ammonium 


Perchlorate Propellants, by Martin Sum- 
merfield and George Sutherland, Project 
Squid, Tech. Rep. PR-87-P (ASTIA AD 
207364), Dec. 1958, 30 pp. (Available 
only on microcards. ) 

Diffusion and Heterogeneous Reaction, 
II: The Catalytic Activity of Solids for 
Hydrogen Atom Recombinations, by Ber- 
nard J. Wood and Henry Wise, Project 
Squid, Tech. Rep. SRD-3-P (ASTIA 
AD 205533), Aug. 1958, 4 pp. (Avail- 
able only on microcards. ) 

Thermal Decomposition of Some Group 
I, II, and III Metal Alkyls, by Louis 
Rosenblum, NACA RM E56L18, Feb. 
1957, 22 pp. (Declassified from Con- 
fidential by authority of NASA Pub. 
Ann. 4, p. 23, 3/16/59). 

Heterogeneous Flash Initiation of Ther- 
mal Reactions, by L. 8. Nelson and J. L. 
Lundberg, J. Phys. Chem., vol. 63, March 
1959, pp. 433-436. 

Cooperative Test Projects of the Joint 
Army-Navy-Air Force Panel on Analytical 
Chemistry of Solid Propellants, Part 6: 
Round Robin 10: Normal Methyl Violet 
Paper Stability Test; Round Robin 11: 
Flame Photo Meter for Sodium, Potas- 
sium and Calcium in Small Arms Pro- 
pellant; Round Robin 12: Dibutylphtha- 
late in Small Arms Propellant, by R. H. 
Pierson, NOT'S 1937, Part 6, Feb. 1958, 82 
pp. 

Elastomeric-binder and Mechanical- 
property Requirements for Solid Propel- 
lants, by Thor L. Smith, Calif. Inst. Tech., 
Jet Prop. Lab., Mem. 20- 178, 1959, 16 pp. 

Analysis of Propyne-allene Mixtures, 
by C. M. Drew, S. R. Smith, A. Jensen, 
MeNesby and A. 8. Gordon, 


NAVORD Rep. 4913 (NOTS 1216), Oct. 
888 


1955. 13 pp. (Declassified from con- 
fidential by authority of NOTS ltr 75061 
B.S.: AM Reg 75-53864, 17 Feb. 1959.) 


Chemistry of Free Radicals Containing 
Oxygen, Part 2: Thermo-chemistry of 
the Hydroxyl and Hydroperoxy! Radicals, 
by P. Gray, Trans., Faraday Soc., vol. 
55, March 1959, pp. 408-417. 


Patterns of Free Radical Reactivity, by 
C. H. Bamford, A. D. Jenkins and R. 
Johnston, Trans., Faraday Soc., vol. 55, 
March 1959, pp. 418-432. 

Studies in Bomb Calorimetry, X: Con- 
ditions for the Combustion of Solid and 
Liquid Fuels, by R. A. Mott and C. 
Parker, Fuel, vol. 38, April 1959, pp. 189- 
204. 


Kinetics of Nitrogen Atom Recombi- 
nation, by John T. Herron, J. L. Franklin, 
Paul Bradt and Vernon H. Dibeler, J. 
Chem. Phys., vol. 30, April 1959, pp. 879- 
884. 

A Method of Successive Approximations 
for Computing Combustion Equilibria on 
a High Speed Digital Computer, by D. S. 
Villars, J. Phys. Chem., vol. 63, April 
1959, pp. 521-530. 

The Conservation Equations for Inde- 
pendent Coexistent Continua and for 
Multicomponent Reacting Gas Mixtures, 
by W. Nachbar, F. Williams and S. S. 
Penner, Quart. Appl. Math., vol. 17, no. 1, 
April 1959, pp. 43-54. 

A Note on Liquid and Solid Propellants 
for Rockets, by E. G. D. Andrews and 
A. W. T. Mottram, J. Roy. Aeron. Soc., 
vol. 63, April 1959, p. 265. 

Hydrolysis of Nitroguanidine and Nitro- 
aminoguanidine, by William R. McBride, 
Ronald A. Henry and Sol Skolnik, 
NAVORD Rep. 2011 (NOT 642), Jan. 
1953, 41 pp. (Declassified from Confi- 
dential by authority of NOTS ltr 7606/- 
BJS: AM 75-54866, dtd 19 Feb. 1959.) 

A Gravimetric Determination of Insol- 
uble Impurities and Ash Impurities in Pro- 
gang Cellulose Nitrate, by Aaron 
N. Fletcher, NAVORD Rep. 6444, Jan. 
1958, 13 pp. 

Cooperative Test Projects of the Joint 
Army-Navy-Air Force Panel on Analytical 
Chemistry of Solid Propellants, Part 5, 
Round Robin 9: Triacetin and Dioctyl- 
phthalate in OGK Propellant, by R. H. 
Pierson, NOT'S 1937, 61 pp. 

Results of Performance Calculations for 
N:0,-N2H, Bipropellant System, by G. 
W. Elverum Jr. and S. Martinez Jr., 


Calif. Inst. Tech., Jet Prop. Lab., Mem. 
20-174, July 1958, 56 pp. 
Combustion in Liquid-fuel Rocket 


Motors, by D. B. caf i Aeron. Quart., 
vol. 10, Feb. 1959, p. 1-27 


The Density of A. G. 
Streng and A. V. Grosse, J. Amer. Chem. 
Soc., vol. 81, Feb. 20, 1959, pp. 805-806. 

Contribution to the Study of Flame 
Propagation in a Closed Tube, by = 
Dalmai, Inst. Franc. du Petrole, Rem. 
Ann. des Combustibles Tiquides, vol. 
Dec. 1958, pp. 1683-1724. (In A og 


Theory of Detonations, III: Ignition 
Temperature Approximation, by C. F. 
Curtiss, J. O. Hirschfelder and M. P. 
Barnett, J. Chem. Phys., vol. 30, Feb. 
1959, pp. 470-492. 


High-temperature Gas Viscosities, II: 
Nitrogen, Nitric Oxide, Boron Trifluoride, 
Silicon Tetrafluoride, and Sulfur Hexa- 
fluoride, by C. P. Ellis and C. J. G. 


Raw, J. Chem. Phys., vol. 30, Feb. 
1959, pp. 574-576. 
Gaseous Detonations, XII: Rotational 


Temperatures of the Hydroxy! Free Radi- 


cals, by G. B. Kistiakowsky and F. D. 


Tabbutt, J. Chem. Phys., vol. 30, Feb. 
1959, pp. 577-581. 

The Mechanism of Deflagration of Pure 
Ammonium Perchlorate, by Raymond 
Friedman, Joseph B. Levy and Keith &. 
Rubmel, "Atlantic Res. Corp. (AFOSR 
TN-59-173; ASTIA AD 211-313), Feb. 
1959, 27 pp. 

Recent Advances in Gaseous Deto- 
nation, by R. A. Gross and A. K. Oppen- 
heim, ARS JourNat, vol. 29, March 1959, 
pp. 173-179. 

Compression Sensitivity of Monopro- 
pellants, by G. A. Mead, ARS Journat, 
vol. 29, March 1959, pp. 192-198. 

Isothermal Compressibility of Liquid 
Oxygen and RP-1, by George T. Y. Chao, 
ARS JouRNAL, vol. 29, March 1959, pp. 
199-203. 

Effect of Initial Temperature on Flash 
Back of Laminar and Turbulent Burner 
Flames, by Burton Fine, Ind. Engn), 
Chem., vol. 51, April 1959, pp. 564-561, 

Predicting Burning Velocities of Lean 
Hydrocarbon Flames, by S. A. Weil, 
Ind. Engng. Chem., vol. 51, April 195%, 
pp. 567-569. 

Smoke Limits of Bunsen Burner Ethy- 
lene-Air Flames, by Joseph Grumer and 
M. E. Harris, Ind. Engng. Chem., vol. 51, 
April 1959, pp. 570-574. 

Determination of the Unreacted Hugo- 
niot for Liquid TNT, by W. B. Garn, 
J. Chem. Phys., vol. 30, March 195%, 
pp. 819-822. 

Formation of Ozone From Atomic 
Oxygen at Low Temperatures, by R. A. 
Ruehrwein and J. S. Hashman, J. Chem. 
Phys., vol. 30, March 1959, pp. 823-826. 

Optical Method of Instantaneous _ 
urement of Flame Temperatures, by A 
Moutet, C. Veret and L. Naudaud 
Recherche Aéronautique, no. 68, Jan.-Feb 
1959, pp. 9-20. (In French.) 

Analysis in Flight of Combustion Gases, 
by J. Fructus, Recherche Aéronautique, 
no. 68, Jan. -Feb. 1959, pp. 49-53. (In 
French. ) 

A Preliminary Investigation of Field- 
induced Ion Movement in Flame Gases 
and Its ee by K. G. Payne and 
F. J. Weinberg, Proc., Royal Soc., vol. A 
250, March 24, 1959, pp. 316-336. 

The High-temperature Oxidation of 
Propane, by J. W. Falconer and J. H 
Knox, Proc., Royal Soc., vol. A 250, 
April 7, 1959, pp. 493-513. 

RF Loss Through Rocket Flame Easy 
to Figure, by R. J. Schmeizer, Space/- 
Aeron., vol. 31, April 1959, pp. 171-172, 


Materials of Construction 


Thermal Shock Resistance—a Design 
Parameter for Ceramic Radome Walls, 
by E. J. Luoma and E. B. Shand, J. Aero/- 
Space Engng. vol. 18, June 1958, pp. 
43-47. 

A Theory for the Ablation of Glassy 
Materials, by Hans A. Bethe and Mac. C 
Adams, J. Aero/Space Sci., vol. 26, June 
1959, pp. 321-328. 

Determination of Static Strength and 
Creep Buckling of Unstiffened Circular 
Cylinders Subjected to Bending at Ele- 
vated Temperatures, by Eldon E. Mathau- 
ser and Avraham Berkovits, NASA Memo 
6-14-59L, June 1959, 18 pp., diagrs. 

Telemetering Missile Data, by W. I 
Keller, Ordnance, vol. 43, May-June 1959, 
p. 1004. 

New Tungsten Know-How Widens 
Design Horizons, by Irwin Stambler 
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Space/ Aeronautics, vol. 32, July 1959, 
p. 48. 

Refractory Ceramic Coatings, by David 
J. Larin and E. E. Mueller, The Trend 
in Engineering, vol. 11, April 1959, pp. 
19-23. 

Ceramic Materials for Nuclear Reac- 
tors, by Lavada A. Moudy and James 
I. Mueller, The Trend in Engineering, vol. 
11, April 1959, pp. 28-31. 

Criteria in the Selection of Materials 
for High Temperature Structures, by W. 
H. Steurer, Symposium on High-Speed 
Aerodynamics and Structures, 3rd, San 
Diego, March 25-27, 1958 (Unclassified 
Papers), vol. 1, pp. 179-227. 

Ceramic Materials for Use in Hyper- 
sonic Aircraft, by John M. Nowak, Sym- 
po-ium on High-Speed Aerodynamics and 
Stiuctures, 38rd, San Diego, March 25-27, 
1958 (Unclassified Papers), vol. 1, pp. 
22')-252. 

Chermal Properties of Nose Cone Ma- 
terials, by I. V. Fieldhouse, J. I. Lang, 
A. N. Takata and T. E. Waterman, 
Symposium on High-Speed Aerodynamics 
anl Structures, 38rd, San Diego, March 
25 -27, 1958 (Unclassified Papers), vol. 1, 
pp. 253-276. 

High Temperature Structures, by N. 
J. Hoff, Symposium on High-Speed Aero- 
dynamics and Structures, 3rd, San Diego, 
March 25-27, 1958 (Unclassified Papers), 
vol. 1, pp. 383-432. 

A General Technique for Evaluating 
the RMAL Shock Characteristics of 
Materials Subjected to Temperature 
Gradients, by David Horwitz, Symposium 
on High-Speed Aerodynamics and Struc- 
tures, 3rd, San Diego, March 25-27, 1958 


(Unclassified Papers), vol. 3, pp. —: 


Spaceflight 


The Nature of Re-entry, by George FP. 
Solomon, AsTRONAUTICS, vol. 4, March 
1959, pp. 20-21, 98, 100. 

Re-entry Heat Transfer, by Lester 
Lees, ASTRONAUTICS, vol. 4, March 1959, 
pp. 22-23, 60, 62, 65. 

Planning a Re-entry and Recovery Test 
Program, by William R. Lucas and 
Myron E. Huston, AstRoNAuTICcs, vol. 4, 
March 1959, pp. 30-31, 88. 

Attitude Control of Space Vehicles, by 
M. Baron T. George, ASTRONAUTICS, vol. 
4, March 1959, pp. 34-35, 106-107. 

Communicating with the Hypersonic 
Vehicles, by Stephen G. Homie and 
Richard L. Phillips, AsrRoNAUTICs, vol. 4, 
March 1959, pp. 36-37, 92, 94, 96, 98. 

The Technology of Manned Return 
from Outer Space, by Hilliard W. Paige, 
J. Franklin Inst., vol. 267, Feb. 1959, pp. 
103-118. 

Theory of the Effect of Drag on the 
Orbital Inclination of an Earth Satellite, 
by John P. Vinti, J. Res., Nat. Bur. 
S'andards, vol. 62, Feb. 1959, pp. 79-88. 

Vanguard I IGY Satellite (1958 Beta), 
by R. L. Easton and M. J. Votaw, Rv. 
Sct. Instr., vol. 30, Feb. 1959, pp. 70-75. 

Conduction Cuts —s During Space- 
craft Re-entry, by Carr B. Neel, Space- 
Aeron., vol. 31, March 1959, pp. 42 2-45. 

Accurate Position Data Needed for 
Space Course Changes, by E. V. ey 
and W. E. Frye, Space-Aeron., vol. 
March 1959, pp. 48-51. 

Plasma Accelerator for Spacecraft Con- 
trol (Design Progress), by Bernard Kovit, 
Snace-Aeron., vol. 31, March 1959, pp. 


-59. 
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The Nation’s Leader in Rocket Engines 
Announces Limited Expansion in 


Some of the programs under study are: 

The maximum specific impulse obtainable (compared to the 
classical infinite volume thrust chamber) for a rocket chamber 
which consists only of an injector and an exhaust nozzle. 

Predicting the performance of a turbine operating with a gaseous 
working fluid which is still reacting and changing composition. The 
effect of boundary layer buildup, leading edge thickness, and radial 
flow on flow through supersonic turbine blades. 

New methods of translating chemical energy into thrust, and 
how they will work in a rocket engine. 

Types of fluid transport systems for rocket engines which offer 
simplicity, efficiency, and light weight. The fundamental param- 
eters which govern the satisfactory operation of bearings and seals 
in corrosive propellant service. 

Determine the separation point in an overexpended rocket noz- 
zle as a function of nozzle divergence angle, and the applications 
of this phenomenon. 

The dynamic interrelations of the various components of a 
rocket engine, and the control systems which best solve the prob- 
lem of controlling a rocket engine. 


Advanced degree in Engineering 

Advanced degree in Chemistry, Physics 

BS degree supplemented by appropriate experience 
For specific consideration of your background in Applied Research 
Write to: 


Mr. Stanley Greenfield, Chief 
Engineering Sciences Research 7 
Dept. 596 A.R. 


: 6633 Canoga Avenue, Canoga Park, Calif. 
‘= 


A DIVISION OF NORTH AMERICAN AVIATION, INC. 


: 


Participation in these combined theoretical-experimental _ 
activities presents a stimulating challenge and involves 
a project-oriented responsibility. 
 Qualific ations include: 
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ROCKET CONTROLS 
ENGINEERS 


1) For development and analytical de- 
sign of valves and controls for rocket 
engine systems. Develop specifica- 
tion requirements for components 
covering design, test and develop- 
ment. Prepare data for proposals. 
BS in ME or AE and 3-7 years’ ex- 

rience with hydraulic pneumatic 
and systems. 


2) For dynamic analysis of mechanical 
systems associated with rocket pro- 
pulsion. Mathematical descriptions 
of controlled rocket components 
with regard to systems stability. 
Interpretation of designs for con- 
trolling rocket components. Ana- 
logue simulation of systems. Su 
vision of above areas. BS, 

EE. Experience in servo systems 
analysis and techniques. 


To arrange interview call collect, 
Niagara Falls BUtler 5-7851, or send 
resume to: 


SUPERVISOR ENGINEERING 
EMPLOYMENT 
BELL AIRCRAFT CORPORATION 
BUFFALO 5, NEW YORK 


TOP SPOTS FOR A 
TECHNICAL SPECIALIST 
AND TECHNICAL LEADER 
ON A NEW SPACE TEAM 


Two excellent openings in Chance 
Vought’s newly formed Astronautics 
Division. 


Senior Specialist in Heat Transfer and Fluid 
Mechanics. 

Advanced A.E., Ph.D. preferred. To 
supervise group in fields of heat transfer, 
pressure distribution analysis, ablation 
theory and general fluid mechanics. Re- 
quires knowledge of real gas effects, 
capability for training people, supervis- 
ing company R & D, keeping abreast of 
technology. 


Physics M.S. with Minor in Meteorology 
(and meteorological experience.) 
To study physics of upper atmosphere, 
including compilation of data. Interpret 
data on atmospheric density, pressure, 
temperature and composition. Keep 
abreast of space probe and satellite 
data. Plan test programs. Recommend 
experimental work and conduct feasi- 
bility studies for weather satellite pro- 


rams. 

Write to: 


Professional Placement Office 
Dept. JP-8 


Evaluating APU Systems for Space- 
craft, by S. L. Sandelman and R. W. 
McJones, Space-Aeron., vol. 31, March 
1959, pp. 86, 88, 90, 92, 94, 97. 

What Instruments for Space Crews?, by 
George W. Hoover, Space-Aeron., vol. 31, 
March 1959, pp. 110, 112. 

Biosatellite Decompression (nomo- 
graphs), by Eugene B. Konecci, Space- 
Aeron., vol. 31, March 1959, p. 117. 

Space Communications, by F. W. Lehan 
and Leang P. Yeah, Space-Aeron., vol. 31, 
March 1959, pp. 120-121. 

Satellite Ascent Paths, by Charles L. 
Keller, Sperry Engng. Rev., vol. 11, Dee. 
1958, pp. 2-14. 

Characteristic Velocity for Changing the 
Inclination of a Circular Orbit to the 
Equator, by Leonard Rider, ARS Jour- 
NAL, vol. 29, Jan. 1959, pp. 48-49. 

A Proposed Kepler Diagram, by Robert 
L. Sohn, ARS Journat, vol. 29, Jan. 1959, 
pp. 51-54. 

Is Newton’s Law of Gravitation Really 
Universal?, by Fritz Zwicky, AsTRo- 
NAUTICS, vol. 4, Jan. 1959, pp. 18-19, 74, 
76. 

Interplanetary Probes: Three Prob- 
lems, by Krafft Ehricke, AstroNaAUTICs, 
vol. 4, Jan. 1959, pp. 20-22, 42, 44, 46. 

An Introduction to Space Navigation, by 
Manuel Fernandez, AsTRoNAUTICS, vol. 4, 
Jan. 1959, pp. 23-25, 76. 

Rocket Astronomy—Window into 
Space, by Herbert Friedman, AstTRo- 
NAUTICS, vol. 4, Jan. 1959, pp. 26-27, 69- 
70. 

Mars—a World for Exploration, by 
Clyde W. Tombaugh, Astronautics, vol. 
4, Jan. 1959, pp. 30-31, 86-93. 

Weather Reconnaissance by Satellites, 
by S. M. Greenfield and W. W. Kellogg, 
Astronautics, vol. 4, Jan. 1959, pp. 32- 
33, 77-78. 

The Earth Satellite Orbit, by Ballard B. 
Small, Astronautics, vol. 4, Jan. 1959, 
pp. 40-41, 46, 48. 

Re-Entry: Problems and Progress, by 
Coleman duP. Donaldson, Astronautics, 
vol. 4, March 1959, pp. 38-39, 108. 

Educational Opportunities in Astro- 
nautics, by Paul Sandorff, AsTRONAUTICs, 
vol. 4, April 1959, Part 2, pp. 2-6. 

Preparing Man for Space Flight, by 
Carl C. Clark and James D. Hardy, 
AsTRONAUTICS, vol. 4, Feb. 1959, pp. 
18-21, 88, 90. 

Isolation, by George E. Ruff, Astro- 
NAUTICS, vol. 4, Feb. 1959, pp. 22-23, 110- 
111. 

Aerodynamics for Space Flight, by E. P 
Williams and Carl Gazley Jr., Rand Corp., 
P-1256, Feb. 1958, 26 pp. 

Space Vehicle Environment, by C. 
Gazley Jr., W. W. Kellogg and E. H. 
Vestine, Rand Corp., P-1335, July 1958, 58 
pp. 

Internal Environment of Manned Space 
Vehicles, by S. H. Dole, Rand Corp., 
P-1309, Feb. 1958, 23 pp. 

Lunar Flight Dynamics, by R. W. Buch- 
heim and H. A. Lieske, Rand Corp., 
P-1453, Aug. 1958, 62 pp. 

Lunar Base Planning Considerations, by 
Holbrook, Rand Corp., P-1436, 
Feb. 1958, 24 pp. 


On Lunar and Planetary Experiments, 
by 8. M. Greenfield, Rand Corp., P-1535, 
Oct. 1958, 9 pp. 


Economic Considerations of Space 
Flight Ground Support Requirements, by 
M. A. Margolis and F. 8. Pardee, Rand 
Corp., P-1589, Jan. 1959, 17 pp. 


Space Flight Gound Facilities Require- 
ments Problems—Launching Facilities, by 
J. J. O'Sullivan, Rand Corp., P-1431, Feb, 
1958, 10 pp. 

Lunar Instrument Carrier—Landing 
Factors, by H. A. Lang, Rand Corp. 
RM-1725 (ASTIA AD _ 112403), jan 
1956, 36 pp. 

Lunar Instrument Carrier—Attitude 
Stabilization, by R. W. Buchheim, land 
Corp., RM- i730 (ASTIA AD 1124/ )2), 
June 1956, 22 pp. 

Recovery of a Circum-lunar Instrument 
Carrier, by Carl Gazley Jr. and David J. 
Masson, Rand Corp., P-1119, Aug. 1957, 
20 pp. 

Lunar Base Study Jury Report: Evalu- 
ation of an Experiment in Creative —- 
Conducted with College Students, 

R. D. Holbrook, H. A. Lang and J. tf 
Huntzicker, Rand Corp., RM-2174 
(ASTIA A 'D 205422), Jan. 1958, 223 pp. 

Outline of a Study of Extraterrestrial 
Base Design, by R. D. Holbrook, Rind 
Corp., RM-2161 (ASTIA AD 156013), 
April 1958, 78 pp. 

Space Flight Trajectories, Navigation 
and Maneuvers, by R. W. Buchhcim, 
Rand Corp., P-1387, May 1958, 45 pp. 

Motion of a Small Body in Earth-Moon 
Space, by R. W. Buchheim, Rand Corp., 
RM-1726 (ASTIA AD 123557), June 
1956, 79. pp. 

Circumlunar Trajectory Studies, by 
H. A. Lieske, Rand Corp., P-1441, June 
1958, 23 pp. 

The Status and Improvement of Physi- 
cal Constants Needed for Precision Tra- 
er by Samuel Herrick, Rand Co Ps 

1559, Des. 1958, 10 pp. " 


Russian Technical 
Articles* 


Solution of the Basic Problems of Ex- 
terior Ballistics, by L. Vorob’ev, 
PMM: J. Appl. Math. and Mech., vol. 
22, no. 3, 1958, pp. 481-492. 

On the Theory of Complicated Systems 
of Stabilization, by A. Iu. 
Ishlinskii, . Appl. Math. and 
Mech., vol. 22, no. 3, 1958, pp. 493-512. 

On the Stability of Motion of a Gyro- 
sc oe on Gimbals, by V. V. Rumiantsev, 

M: J. Appl. Math. and Mech., vol. 
22, no. 3, 1958, pp. 513-520. 

On a Gyroscope Mounted in Universal 
Suspension (on Gimbals), by N. 
Chetaev, PMM: J. Appl. Math. and 
Mech., vol. 22, no. 3, 1958, pp. 521-532. 

On a Class of Motions in Magneto- 
Hydromechanics, by V. N. Zhigule v, 
PMM: J. Appl. Math. and Mech., vol. 22, 
nol 3, 1958, pp. 537-539. 

On Gas Flow in Laval Nozzles, by O. S. 
Ryzhov, PMM: J. Appl. Math. and 
Mech., vol. 22, no. 3, 1958, pp. 549-553. 

Computation of a Laval Nozzle, by A. 
Sh. Dorfman, PMM: J. Appl. Math. and 
Mech., vol. 22, no. 3, 1958, p. 554. 

Analysis of Periodic Motions in Elec- 
trical Servomechanism Vibration Loop 
With Constant Disturbance, by I. N. 


*The Editors of Technical Literature 
Digest are making a systematic search for 

rtinent Russian articles, both in the orig- 
inal Russian and in translated form. For 
the balance of 1959, these will appear in 4 
separate section with this heading, mainly 
to draw the readers’ attention to the new 
listing. After that, they will appear item 
by item under the proper subject headings. 
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Krutova, Avtomatika i Telemekhanika, vol. 
20, no. 5, 1959, pp. 564-575. (In Rus- 
sian. ) 

On the Invariance Principle in Auto- 
matic Control Combined Systems, by 
Vy. I. Dunaev, Avtomatika i Telemekhanika, 
yol. 20, no. 5, 1959, pp. 591-594. (In 
Russian. ) 

Electric Drive With Powder Clutch and 
Its Application in Automatic Control 
Systems, by G. G. Konovalov and Ya. I. 
Flid, Avtomatika i Telemekhanika, vol. 
20, no. 5, 1959, pp. 657-662. (In Russian.) 

Stabilizing Current L-C Devices, by 
B. E. Kubyshin and A. N. Miljakh, 
Avtomatika 1 Telemekhanika, vol. 20, no. 
5, 1959, pp. 663-668. (In Russian.) 

On a Stability Criterion for Nonlinear 
Control Systems, by Chang Szu-Ying, 
Aviomatika i Telemekhanika, vol. 20, no. 
5, 1959, pp. 669-672. (In Russian.) 


Concerning Absolute Stability of Cer- 
tain Control Systems, by Yu S. Sobolev, 
Aviomatika i Telemekhanika, vol. 20, no. 
4, 1959, pp. 401-405. (In Russian.) 

Integral Estimate to Choose Optimum 
Parameters of Automatic Control Systems 
With Present Overshoot, by A. L. Tupit- 
syn, Avtomatika i Telemekhanika, vol. 
20, s 4, 1959, pp. 406-414. (In Rus- 
sian. 

Application of Square Integral Estimate 
for Finding Optimum Parameters of Pilot 
With Rate Feedback, by V. D. Matytsin 
and V. A. Ryapolov, Avtomatika i Tele- 
mekhanika, vol. 20, no. 4, 1959, pp. 415- 
421. (In Russian.) 

Dynamics of Electric Servomechanism 
Vibration Loop Under Free Oscillating 
Conditions, by I. N. Kritova, Avtomatika 
i T'elemekhantka, vol. 20, no. 4, 1959, pp. 
422-436. (In Russian.) 

Compensation of Continuous Automatic 
Control System by Means of Delay Ele- 
ments Filter, by Sin-min Van, Avtomatika 
i T'elemekhantka, vol. 20, no. 4, 1959, pp. 
437-446. (In Russian.) 

Investigation of Pneumatic Jet Element 
Characteristics, by L. A. Zalmanzon and 
V. A. Semikova, Avtomatika i Telemek- 
hanika, vol. 20, no. 4, 1959, pp. 447-467. 
(In Russian. ) 

Design Principles of Contactless Sys- 
tems of Remote Control With Exponential 
Converters, by V. A. II’in, Avtomatika i 
Telemekhanika, vol. 20, no. 4, 1959, pp. 
468-472. (In Russian. )] 

Analysis of Operation of Pulse Mag- 
netic Elements in Contactless Devices of 
Remote Control, by I. V. Prangishvili, 
Avtomatika « Telemekhanika, vol. 20, no 
4, 1959, pp. 473-479. (In Russian.) 

Concerning Stability of a Pulse Moni- 
tored Servosystem, by Ludvik Prouza, 
Avlomatika i Telemekhanika, vol. 20, no. 
4, 1959, pp. 518-520. (In Russian.) 


Analysis of Stability of Distributed 
Parameters Automatic Control Systems 
With Loss, by Ya. B. Kadymov, Avto- 
matika i Telemekhanika, vol. 20, no. 4, 
1959, pp. 525-527. (In Russian.) 


Bibliography on Automatic Control and 
Related Problems for 1956, Avtomatika i 
Telemekhanika, vol. 20, no. 4, 1959, pp. 
528-532. (In Russian.) 

Physical Phenomona That Occur When 
Bodies Compressed by Strong Shock 
Waves Expand in Vacuo, by Ia. B. 
Zel'dovich and Iu. P. Raizer, Soviet Phys.- 
JETP, vol. 35 (8), no. 6, June 1959, p 
980. 


The Stability of Shock Waves in Mag- 
netohydronamics, by S. I. Syrovatskii, 
Soviet Phys.-JETP, vol. 35 (8), no. 6, 
June 1959, p. 1024. 
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DEGREES 


... another plus benefit at the Los Alamos — 
Scientific Laboratory. 


Graduate programs for advanced degrees in 
nuclear, mechanical and electrical engineering, 
chemistry, physics and mathematics are. 
available to qualified employees of the 
Laboratory at the Los Alamos Graduate Center 
of the University of New Mexico. With 
tuition support from the Laboratory, <3 
employees may meet requirements for masters’ 
degrees entirely in Los Alamos evening 
classes, and for Ph.D. degrees with courses 

at Los Alamos and a minimum of one 
academic year at the University of New 
Mexico campus. Frequently, research for 
dissertation may be done in the Laboratory 
with some of the world’s finest facilities _ 
‘and equipment. 

College graduates in engineering and scien 
who want fascinating research assignments 
while continuing their advanced education | 
may write to: 


Director of 
Personnel 
Division 59 


OF THE UNIVERSITY OF CALIFORNIA 
LOS ALAMOS, NEW MEXICO 
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SOUTHWEST 


Ad 


SELF-ALIGNING BEARINGS 


ROD END 
TYPES 


PLAIN TYPES 


PATENTED U.S.A. 
World Rights Reserved 


CHARACTERISTICS 


ANALYSIS RECOMMENDED USE 
. For types operating under high 
Stainless Steel Ball and Race temperature (800-1200 degrees F.). 


Chrome Alloy Steel Ball { For types opera 
Race ultimate loads 

Bronze Race and Chrome For types operating under ser Seats 

Steel Ball with minimum friction requirements. 


Thousands in use. Backed by years of service life. Wide variety of 
Plain Types in bore sizes 3/16” to 6” Dia. Rod end types in similar 
size range with externally or internally threaded shanks. Our Engi- 
neers welcome an opportunity of studying individual requirements 
and prescribing a type or types which will serve under your demand- 
ing conditions. Southwest can design special types to fit individual 
specifications. As a result of thorough study of different operating 
conditions, various steel alloys have been used to meet specific 


SOUTHWEST PRODUCTS CO. 


705 SO. MOUNTAIN AVE., MONROVIA, CALIFORNIA 


needs. Write for Engineering Manual No. 551.Address Dept. ARS-59. 


RESEARCH ENGINEERS 


IF YOU— 


1. are experienced and interested in dynamic analysis 
and vibration engineering, 


have M.S. or Ph.D. degree or equivalent in specific 
experience, and have United States citizenship, 


desire to apply your knowledge to a wide variety of 
technically significant industrial and government 
problems, 


would enjoy a permanent staff position ensured by 
multiple sponsorship, 


would appreciate the opportunity to proceed with 
your graduate studies on a convenient basis while 
functioning as a full-time staff member, 


enjoy working in communications and close pro- 
fessional contact with men from other scientific 
and engineering fields, 


would like to join the staff of one of the largest in- 
dependent research organizations in the world, 


please contact: 


E. P. Bloch 


ARMOUR RESEARCH FOUNDATION 


of Illinois Institute of Technology 
10 West 35th Street 
Chicago 16, Illinois 
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